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"Development Of a Generalized Digital Simulator Bor Production Systems" 

In this study a computer simulator has been developed to 
facilitate the planning and control of the production systems. The 
proposed simulator has capabilities to handle various production plann- 
ing and control functions e.g., sales forecasting, inventory control, 
etc. The production system under study is simulated over a planning 
horizon to generate a set of decision rules. These decision rules 
provide guidelines for the production manager to take logical decisions. 
As such the simulator cm handle a production situation where fifteen 
•work stations, fifty components and three finished products are involved. 
This restriction has been imposed keeping in view the capabilities of 
IBM 7044 computing system. The simulator has been validated with the 
help of data collected from an electronic industry manufacturing tele- 


vision sets 



CHAPTER I 


I2TTBDDUCTI0N 

Production Management is the art of achieving the four basic 
aims of management - i) the right product, ii) at the right time, 
iii) at the right cost and iv) at the right place. The achievements 
of these goals can be made possible only by a proper coordination of 
planning and control of the production system. The primary function 
of production planning and control activity is to design a sound 
framework within which a set of inter-related decisions can be satis- 
factorily made. This necessitates the design of control mechanism 
incorporating the following Ingradients (l) : 

1. A current operation plan to specify the desired per- 
formance of the manufacturing system in terms of 
meaningful effectiveness criteria. 

2. A data acquisition system to determine, on a timely 
and accurate basis, the current status of the manu- 
facturing system. 

3« A method of comparison of the actual results with the 
desired goal. 

4. The corrective action to activate the goal seeking effort. 

The last two are the manifestations of the feed-back loop 


concept. The inherent complexities of the real system make the 
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design of control systems a difficult task. However, maay research, 
workers and production managers have attempted to design idealized 
control systems, keeping in view the needs of individual industries 
(_ 1 » 2 )• But, unfortunately very few efforts have been directed towards 
the control of real life production systems. 

Efforts of various re sear as workers xn the f ield of produc- 
tion plaining and control has helped the Industries to design integra- 
ted control systems for individual's requirements. Bit, in Indian 
industries even the basic concepts of production planning and control 
are not in vogue. The present work aims at developing a production 
planning and control system for an Indian industry. The model has been 
developed and validated for tire manufacture of Television Sets at 
J.K. Electronics. The particular production situation was selected 
because of the following advantages : 

1. Availability of accurate aid timely data. 

2. A mathematical model can be easily structured, because 
the size of the system is neither too big nor too 
small. 

3. Since the manufacture cf TY sets involves a continuous 
assembly line, the interactions of various elements 
within the system are easy to identify and evaluate. 

4. At present the production managers of TY industries in 
India are faced with an uphill task of systematic 
design of production planning sad control systems to 
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meet the heavy demand of the product In the country and 
abroad. It is ore of the fast expending industries 
which requires careful attention on the p art of the 
management. 

Traditionally tv.o basic approaches have been employed by 
the industries to design the production control systems. One wey is 
the ordinary individual control function approach in. which several 
control functions are treated individually. Another way is the case - 
study approach in which a similar real or hypothetical system is 
analyzed. Eoth the approaches suffer on account of several drawbacks 

( 3 ) : 

1. Analysis is given more emphasis compared to the design. 

2. They do not consider the dynamic nature of the produc- 
tion environments. 

3. A comprehensive understanding of the vhole production 
control system is lacking - the interactions are either 
undefined or not well defined. 

4. The cause "•effect relationship of decisions is not 
incorporated and thus lack in evolving corrective 
action. Further, the feed back concepts have not been 
given adequate representation. 

5. The role of management information systems in control 
has either not been realised or not stressed. 
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Most of tiiese drawbacks are offset in a third and compara- 
tively new technique - the computer simulation. Simulation is the 
process of conducting experiments on the model of a system in lieu 
of either direct experimentation with the system itself, or direct 
analytical solution of some problems associated with the system. Thus, 
with a set of alxematives in hand, simulation can be resorted to, in 
order to study the system’s behaviour and make pertinent decisions 
about the choice of the alternatives. However, simulation does not 
necessarily pinpoint optimum solutions, but it generates a response 
curve on which the optimum lies (4 ). 

The present work consists in the development of a production 
planning and control simulator for an Indian industry. A thorough 
search of the literature indicates that very few computer simulators 
have been developed by various research workers for pedagogic purpose 
(3, 5, >6, 17 ). More prominent amongst those extensively used is, 

PSD SHI V - Production System Simulator ('l). It has been successfully 
used by the instructors and research workers as an effective teaching 
aid in the field of production control. But, the simulator has basi- 
cally been designed to examine the knowledge of the students in 
designing production control systems. It is used as a tool to evar- 
luate the various decisions made by the user, but it does not provide 
guidelines to arrive at those decisions. It is a general purpose 
simulator which can be utilised for different simplified production 


systems 
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In the present work, the PBQSM 7 simulator has "been modi- 
fied considerably so as to represent the mathematical model of the 
real life production system under study. The model has been developed 
for the aforementioned industry with the basic objective of designing 
a control system that will provide the production manager with deci- 
sion policies which, on implementation in the long run, will result in 
minimum cost of production. The production planning and control sys- 
tem encompass the following functions : (a) forecasting, (b) Operations 
Planning, (c) Inventory Planning and Control, (d) Operations Scheduling 
and (e) Dispatching and Progress Control. The decision making mechanism 
covers the above functions and also takes into account the interactions 
amongst the various system elements e.g. different products, raw 
materials, work stations, men, etc. to provide optimum decision poli- 
cies to the decision - maker. 

In order to simplify the mathematical formulation of the 
model, appropriate assumptions have been male in the system, as and 
when necessary. The model has been tested using the data collected 
from J,K* Electronics for only one particular product out of several 
products maruf actured by the said industry. Although the production 
planning aid control simulator presented here has been primarily de- 
veloped and tested for a television industry, every care has been 
taken to make it general purpose simulator. The same simulator, with 
minor modifications, can be adapted for aiy other production situation 
without much difficulty* 



CHAP1ER II 


HEERAECRE HE VIEW 

As was stated in the previous chapter the use of computer 
simulation, as a technique for designing of production control sys- 
tem, is still in its infancy. However, in the past efforts for the 
development of production planning and control simulators have been 
made by educational institutes. In this chapter a review of the va- 
rious educational production simulators and their applications to the 
real life production systems has been made. 

2. 1 Educational and leaching Efforts 

Currently a lot of interest is being shown in the use of 
computer simulation as an aid in teaching concepts of industrial en- 
gineering, production management and operations research. Computer 
simulation has also been used successfully for studying higher-order 
relationships between major components of a firm aid between competing 
firms. One of the first efforts in pedagogical simul afcion was the 
work on Industrial Dynamics by J.W. Forrester (6). Forrester developed 
a computer language DYNAMO (DYNAM IC IMPELS) , .with the capacity to 
account for interactions amongst system elements and the economics 
of the entire industry. 

Roberts (8) used Industrial Dynamics philosophy for the 
control of the research and development projects. DINAMO has been 
used to evaluate and predict dynamic behaviour of production - inven- 
tory systems with changes in market demand and other factors (6, 9). 
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Several simulators have been developed for "management games". 
Peril aps the most comprehensive of these is the Carnegie Tech. Manage- 
ment Game (5)* 'Management games’ are concerned with the production 
system only in a gross way. They are used primarily for training 
purposes (id). 

Porter ex. al. (ll), Vollman (12), and Whiteman and love 
(l3) have used simulation as a teaching aid in Industrial Engineering. 
One of the latest efforts in this field being that of Greenlaw and 
Hottenstein (•*?) . They have integrated the vrorks of Porter, Tollman, 
Whiteman et. al. (.14, 12, 13) and other techniques of production 
management (l4, 15, 16, 47) to develop a computerized production manar- 
gement simulator known as PBOSBl (PBPHJCTIOl MANAGEMENT SIMDIATIOE). 

As compared to DIIIAMO and the man cement ganes PROSBI stands 
out on account of the following reasons : 

1. DYNAMO simulates the system with a macro-view, vhile in PBDSIM 

a micro-view is followed by considering the hourly events in the 
simulated production system. 

2. While DYNAMO is capable of simulating a complete industrial sys- 
tem or entire economy, it treats production system as a single 
component in the total simulated system. PBDSIM emphasizes 
centrally oa the production system and generates relationship 
with the environment. 

3. Unlike the various management games which are concerned with 
the production system only in a gross way, PBDSIM can. be used 



as a management gane emphasizing the production system in 
details. 


8 


Owing to the above mentioned advantages the PBOSIM model has 
been widely used by production management students and research wor- 
kers at the Pennsylvania State and other universities. Computer pro- 

* 

gram written in POKE EM has been developed for PB0SB1 model for two 
different IBM systems - the 700/7000 and 360 series. This computer 
package can be easily used by instructors having little or no know 1 - 
ledge of computer programming. 

2.1.1 PBOSIM Model 

The fundamental objective of PBOSIM is to facilitate the 
learning and the mastery of certain basic concepts aid techniques of 
production management by the participant. In keeping with the above, 
the participant is provided with a simulated environment wherein he 
can make decisions and study their effects. In accordance with the 
objective, decision problems in other functional areas of the business 
firm e.g. marketing and finance, have been included in the model only 
to the extent that they bear relevance to production management. 

The model has been designed so that the decision problems 
faced by the participant, possesses many of the sane fundamental 
characteristics as these faced by real world pro duetionm angers. 
PBOSIM is capable of simulating an enterprise, manufacturing at most 
three different products (X, Y, and Z respectively). The manufactu- 
ring operation is assumed to be comprising of two production lines 
*P0ETEM and POBTBM IV are used interchangably in this text. 
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each o+ -which consists of four identical machines as shown in Pigure 1. 
Each of the three products require processing on each of -these two 
lines. Saw materials are first fabricated into unfinished products cn 
line 1, referred to as 2', Y' and Z' . Then these are converted into 
finished products on line 2. 

The participant has to make numerous inter-related decisions 
during each period of simulation, i.e. a dcy. Two fundamental types 
of relationships exist in the model, as is illustrated graphically in 
Pigure 2 : 


*Derived from n PEDS32- r : A Production Management Simulation" by 
Greenlaw, P.S. and K.P. Hottensteinj page €. 
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1. The participant's decisions mil determine the number of units of 
output of finished goods to be produced by the fine. ±1 order to 
meet its demand schedules. 

2. She seme decisions will also determine the costs incurred during 
each period of the simulation. 

So arrive at these sets of decisions, the PBOSIM participant 
has to decide his objectives aid subobjectives. She main objective in 
the model is that of minimi zing the total costs involved in operating 
the sysxem in the long run. She term "long run" is important since in 
many cases it may be necessary for the firm to incur additional • costs 
in cue period in order to avoid even larger costs in future periods. 

The cost approach is followed becaise -the cost can be measured much more 
accurately than the profit, especially at the subsystem level. 

In order to meet the basic objective, numerous subob j ectives 
have to be established. for exanple, keeping inventory costs under con- 
trol, meeting the demand schedules, products assignment to machines for 
minimum set-up costs etc. The PBOSIK participant decides as to (i) how 
such sub-objectives ought to be formulated aid (ii) how various measures 
of performance of the system may be developed. 

The final decisions which the PHD SIM participant makes to 
satisfy previously mentioned objectives are as follows (figure 2): 

1. Quality control expenditures. 


2. Plant maintenance expenditures. 
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3. Regular raw materials order in units. 

4. Expedited raw materials order in units. 

5. Scheduling the operators on the machines. 

6. Decision for training the operator. 

7. Product scheduling on various machines. 

S. lime scheduling of various machines in hourly units. 

PROSE! model incorporates the various performance - measures 
and sets of decision rules in order to facilitate achievement of the 
above aims (i7). Host o f these decision rules are complex and inter - 
related. Various mathematical techniques are attempted in the model 
e.g. the technique of dynamic programming has been used for planning the 
production for short term. 

Greenlaw and Hottenstein ( 7 ) have developed a number of 
data - forms to be filled in by the participant to arrive at his final 
decisions. This reduces the number of computations and provides guide- 
lines to the participant in making final decisions. After feeding these 
final decisions as data in PROSES computer simulator, the participant 
receives a computer printout providing following informations relative 
to its operations, for each sinulation period : (i) cost information, 
(ii) product information, (iii) inventory information, and (iv) demand 
information. 

AL though there are many thousands of combinations of deci- 
sions which may be made in each period, the PROSES participant is 
constrained in his choice of strategies. The generation end evaluation 
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of decision strategies is restricted due to limited planning horizon. 

In the PB0S33I no del the problem of optimization with res- 
pect to the attainment of various objectives is left unsolved, that is, 
one objective nay be more full" net only at the expense of lesser 
attainment of another. Also generally no sroyle decision algo ri tin 
xs available to completely optimize the objective-. Due- xo these com- 
plexities in decision melting process, the PEDSBl participant is 
generally content to satisfice rather than optimize his objective, 

Ifize et. al. ( "3 ) extended the PROiTF-I model at Auburn 
University to develop a production system simulator, PROSIII V, -iiich 
is capable of simulating a variety of operating environments, inclu- 
ding norm anufaetu ring operations. vTnile PROSE- simulates tiie hourly 
activities of the production system, in PBDSBI V each minute event is 
simulated, mailing ix more accurate. 'The model incorporates move 
realistic situations like -.vaiting lines end accounts for the inherent; 
variability in different parameters, such as number of v,o rk stations, 
number of pioducts, etc. 

2,1.2 PROSE! 1 Ilodel 

PEL'SEI 7 is capable of simulating a production planning 
and control system consisting of the following functions : 

1. Sales forecasting. 

2. Operations planning. 

3. Inver.toiy planning and control 

4. Operations scheduling. 

3. Dispatching and progress control. 
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The model simulates a typical inventory-pro duct ion - sales 
system as shown in Figure 3. Such a system usually has the following 
characteristics : 

1. Several finished products, sold in discrete units. 

2. Periodic (weekly, monthly, etc.) demand for each product 
is a random variable and may or may not fellow a trend. 

3. Each product is composed of assemblied and parts. 

4. Atleast some of the assemblies are composed of sub- 
assemblies and purchased parts. 

5. There are some common components aid subassemblies 
among the finished products. 

6. lead time for purchased parts is a random variable. 

7. Fabrication and assembly operations are performed at 
"work stations". A work station is defined as ai area 
where specified work is performed; this may consist 

of a man and/or a machine, or several men and/or 
machines. 

8. Different assemblies require processing on some of the 
same work stations. 

9. Processing times at all work stations are essentially 
deterministic . 

The above charactersties are typical of many manufacturing 
systems. A particular system is simulated by initializing pertinent 
parameters in the computer model. 




PIGUBE 3 : TYPICAL IHT/MTOEY-PIiDDUCTIOSr-SALES STSTMl* 

Important characterstics of any production simulator are 
the size of operating system it can simulate and the computer time 
required for a typical simulation run. In this model these eharac- 
terstics depend upon the number of work stations, number of finished 
goods, and the number of stock numbers (a unique stock number is 
assigned to each finished good, assembly, pare, and raw material) , in 
the production system. PH) SIM V is limited in the size of the prob- 
lem it can handle by the size of available computer memory. The 
computer time required for the simulation depends entirely upon the 
size of the manufacturing system being simulated and the computational 
speed of the computer. 

* Derived from "PBDSIM V Administrator's Manual s Production System 
Simulator" by Mize, JJi. et. al; page- 6. 
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.Another important characterstic of this kind of simulator 
** is its flexibility. The flexibility of PB03III V is shovan by the 
following features : 

1. The operating period is variable^ usually a period 
of one week is used. 

2. Several periods of simulation can be run on one 
computer run. 

3. The time increment is variable. If a time increment 
other than one minute is used, then all processing 
and set - up times must be in multiples of the time 
increment . 

4. The lot size is variable. Items may be produced in 
lots of any integer size. 

5. Certain features, such as random demaid and random 
lead time cm be turned 51 on or off” with a parameter 
card. 

6. No more than five work stations may process any one 
stock number and no more- than five raw materials or 
sub-assemblies may merge into one stock number. This 
restriction is easily cvercone by specifying dummy 
assemblies. 

2. 1*2.1 Conceptualization of PEOSEI V 

In PSOSIM V the instructor acts as the Production Manager. 
He translates a hypothetical manufacturing _ystem into pertinent 
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parameters and the simulator is initialized. Then the users or stu- 
dents, analyzes this data and designs a preliminary control system, 
to make a set of operating decisions. One set of user’s decisions 
directs PBDSIM V through one period of simulated operation. 

Production orders for manufactured items and purchase or- 
ders are placed only at the beginning of each period. Each production 
order is immediately released to the first work station that performs 
work on that item. The new production orders follow a EIEO system. 
PBDSIM 7 divides the weekly operating period into five equal daily 
periods. It causes one -fifth of weekly demand to occur instantaneously 
at the end of each day. It also divides the weekly operating time for 
each work station by five and sets up five identical days of available 
operating time. 

In PBDSIM 7 a complete record of dynamic inventory is 
maintained. Whenever the stock on hand of any component reaches a low 
level, it can be replenished only by instructing the simulator to 
manufacture or purchase (depending upon the typt of item), a certain 
number of components. 

In PEDSIli V each work station is assigned a unique identifi- 
cation number. Each manufactured item must be processed on atleast 
one work station. This processing may or may not require the addi- 
tion of other stock numbers. The operations network for each item 
is defined by specifying the sequence of the work stations required 
to manufacture the item. The same work station may not appear more 


than once in the same network 
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BTGUEE 4 : CONCEPTUAL VIEW QP PEOSDI V SUBSYSTEM* 

Conceptually 5 a work station is considered as a Qxeueing 
Subsystem as shown in Figure 4. The exact number of " queue positions” 
for each work station is provided as an initialization p arameter. A 
queue position holds one production order, regardless of the size of 
the order. Queue positions are not used as temporary storage for 
items needed in the manufacture of production order, all such items 
are held in inventory. Simulator continuously moves orders forward 
in the queue. If there are no orders in the queue at a particular 
time, then the work station goes idle until another order arrives. 

PBDSIM V processes production orders through work sta- 
tions in "lots”, rather than in individual units. The ’’lot size” 
is an initial parameter and remains unchanged throughout the entire 
simulation. The sane lot size applies to all manufactured items. 

If a product order quantity other than the lot size is ordered, 

PROSE! V truncates the order to the next lowest multiple of the 
lot size. A batch quantity is also specified for each manufactured 

*De rived from "PEOSIK V Administrator's Manual : Production System 
Simulator" by Mize, J.H, et. aL.; page 17. 
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item in the initial data. Most items have a batch quantity of one, 
whereas small items such as screws, nuts, bolts, etc., which are 
manufactured and used in large quantities nay have any desired batch 
quantity. A ’’batch” is regarded as a ’’unit 1 ’ until it has been pro- 
cessed on the last work station in its network. All processing times 
on v,ork stations are specified for a "batch” of the particular item. 

When a production order for some stock number comes in the 
first queue position, the simulator checks if the raw material needed 
to manufacture one lot of the stock number are available in the in- 
ventory. If they are available then the production order is moved 
from queue position one into the work station. The lot remains on 
the work station for a known length of time (processing time). Simu- 
lator moves the processed lot from the work station to the "Hold Block” 
as shown in Figure 4. A hold block may be regarded as an area adjacent 
to the work station used to collect a specified number ("hold quantity") 
of units of an order before sending the units either (i) to the nest 
work station in the network, or (ii) to inventory, if This is the last 
work station on hie network. Simulator 'then determines whether enough 
of the raw material, needed are available to process another lot of the 
same stock number. If they are available, another lot is moved from 
queue position one onto the work station. This process is continued 
until hie entire order is processed or until items shortage is faced. 

Whenever a production order with insufficient quantity of 
raw materials, needed to produce one lot, is encountered the simulator 
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continues searching the queue for an order that can be processed* The 
orders with insufficient inventory quantities are moved back in the 
queue. If no order with sufficient inventory is found , the work sta- 
tion becomes idle until the Inventories are available for processing 
any of the orders in the queue* 

Purchased items are received only at the end of a day and 
are available for use at the beginning of the next d^ r . Manufactured 
items are added to stock on hand in lots imm ediately upon the comple- 
tion of processing at the final work station in the network* These 
items are available for immediate use , soon after the processing is 
completed for one lot. 

If all queue positions are full at a work station, then pro- 
cessing at the preceding work station in the network stops* It remains 
idle until the items in the hold block can be moved to the next work 
station. When a work station has completed processing an order, it Is 
forwarded to the next work station in the network, whether or not the 
n Hold Quantity” has been reached. 

Each production order is assigned a sequential ”order number” 
by the simulator. It is needed in order to keep separated two or more 
production orders for the sane stock number, placed in same or diffe- 
rent periods. In sane period more than one production order for one 
item may be placed due to the following reasons : (i) economic pro- 
duction lot sizes, (ii) priority assignment to certain orders and 
(iii) to force production of 1? at least a few” of several items in a 
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period. When an order arrives at a work station for processing, ihe 
simulator searches the queue for an order with the sane order number, 
if found it places the order at sane queue position. If no such 
match can be found, the arriving order is assigned the first available 
queue position. 

2. 1.2.2 Simulator Initialization 

PB3SH! V is initialised by providing values of various 
variables and parameters, llize et. al. ( 3 ^ have used the following 
parameters, pertinent to the manuf acturing system being simulated for 
initialization : 

1. Woric station network for each manufactured item. 

2. Parameters for demand generators (one demand generator 
for each finished product). 

3. Parameters lor lead time generators (one lead time 
generator for each purchased item). 

4. Process time and set up time for each manufactured 
item at each work station. 

5. Ill costs, parameters and constraints pertinent to 
the simulated system. 

PEDSIM V user is provided with the following data as a 
starting point for his analysis : 

1. Items 1, 4 aid 5 from the preceding paragraph. 

2. A demand history for each finished product. 

3. labour and work station rates for each work station. 
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shift change cost* overhead rate, idle time cost, etc. 

4. Dor each stock number, initial stock on hand, carrying 
cost, reorder cost, discount order quantity, regular 
price, discount price, average lead time, out-of- 
stock-cost, etc. 

The user uses the above mentioned data to design his pro- 
duction control system. Ihe following general procedures nay be used 
as a way of getting started : 

1. Begin the design of the control system. 

2. Analyse the historical sales data and determine the 
forecasting procedures. 

3. Determine the inventory system to be used aad the eco- 
nomic order quantity or time interval for each purchased 
part. 

4. Determine the economic batch size for each manufactured 
item and use this to help determine the hold quantity 
for each product at each v/orl; station, aid/or the 
production order quantity. 

5. Determine the point at which it is economical to switch 
from overtime to another shift for each work station. 

6. Analyze the demand trends over a long period of time 
(say one year) and attempt to "smooth." total produc- 
tion requirements over this period. 
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She user' s control system must lead him to the following 
decisions for each simulated period : 

1. Estimation of demands for all the finished goods. 

2. Purchase orders for raw materials. 

3. Production orders and desired sequence for manufactured 
items. 

4. lime available for each work station (regular time, 
overtime, extra shifts). 

5. In-process buffer sizes (hold quantities). 

2. 1.2.3 Simulator Ground Rules 

The following "ground rules" govern the manner in -which 
PROSIK V acts upon the user's decisions and simulates the operation 
of the manufacturing system : 

1. Simulator reads the forecasted demand values for fini- 
shed goods. It then generates an actual demand quantity 
for the wee k using demand parameters and the Monte Carlo 
Sampling (IS). To obtain the daily demand, the weekly 
demand is divided by 5 and the value is truncated. Con- 
ceptually demand occurs instantaneously at the end of 
each day. 

2. Por each purchase order, a random lead time is gene- 
rated using paraaeters (mean lead time aid its stan- 
dard deviation) and Monte Carlo Sanpling. Conceptually j 
all purchase orders for the week are placed at the 
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beginning of the week and lead time is measured from 
this time. On receipt of the order, the components cm 
be used in production only at the beginning of the next 
day and not on the same day. No carrying cost is char- 
ged for newly arrived components until beginning of the 
day following the receipt. It is possible for a pur- 
chase order to arrive during the same week in which it 
was placed. 

3. Simulator assigns a unique order number to each product- 
ion order it reads, which signifies the order as it 
progresses through its network of work stations. If all 
available queue positions are full at a particular work 
station, PBDSIII 7 prints a message to this effect and 
that production order Till not be processed. It will 
have to be placed again next week. 

4. For each period of simulation the number of hours of 
operation for each work station are specified as 
"lime Avail dole" and "Forced Overtime" values. These 
values must be specified in multiples of 5 to provide 
5 equal daily periods to the simulator. In deciding 
how long to operate each work station, the decision - 
maker attempts to estimate the total time requirements 
for each work station. The voik stations remain ope- 
rable until the "Forced Overtime" time value, even if 
it goes idle due to lack of production orders. Each 
work station must operate for at least one full shift . 
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i.e. 40 hours/week; and xt may operate upto a maximum 
of three full shifts during a period of five days. A 
"shift change cost” is charged vhen&ver there is a 
shift; change. 

5. The user may change the "Hold Quantity" of any stock 
number - work station combination during any period. 

This provides him the flexibility in how he moves 
partial orders through sequential work stations. He 
may move each lot of the item, as it is manufactured, 
either to the next work station or he c ax accumulate 
several lots. It also provides him with a realistic 
mechanism for smoothing material flow and for achie- 
ving high utilization of equipment. All "Hold Quan- 
tities" are given in lot sizes. 

6. PBDSIM Y simulates one day's production activities and 
then attempts to satisfy the generated daily danand 

for each finished goods. If stock on hand of a finished 
good is not sufficient to meet the danand, a back 
order for the item is generated. Back orders must be 
filled from production before the new daily demand is 
satisfied and before stock on hand can move above zero 
level. Backordered units and out of stock costs are 
calculated after all transactions (production, demand, 
receipts) for the day have occured. 
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7. Carrying cost for each stock number is computed daily. 

r 

This is determined by multiplying l/5 of weekly carry- 
ing cost with the stock-on-hand at the end of each 
day, after all transactions have oeeured. Daily carry- 
ing costs are then summed for the week. 

8, The user may enter ” Additional Sales” of any stock 
number during any period in the system. These sales are 
in addition to the regularly generated demand values. 
Additional sales are subtracted from stock-on-haad at 
the beginning of the week. 

Since PEOSIM Y is basically designed as a teaching aidj the 
instructor can play a major role in controlling and influencing the 
environment in vhich The student's control system must perform. He 
can interject "noise" into the system at ary time he wishes. Speci- 
fically, the following external disturbances car be interjected into 
the system. 

1. Change in parameters in the demand generator for fini- 
shed products. There are many types of changes which 
may be made (slope, period of cycle, variance, etc.). 

2. Additional sales order may be placed for any stock 
number during the period. This is useful for empha- 
sizing the importance of safety stock. 

3. Change of mean and/or standard deviation of lead time 
generator. This emphasizes the need to keep running 
statistics on certain important variables in the system. 
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4. Change in any other paraaeter or cost In the system. 

The simulator attempts to test the ability of the use’r's 
control system to respond to and adjust for external perturbations 
and other disturbing factors. 

2. 1.2.4 Cost Ho del 

The basic objective of PBOSIM 7 is to design a control 
system that will lead to decisions which, in the long iun, will 
result in lowest total cost. A number of periods are simulated, by 
declaring decisions for the different periods in a single computer 
run. Simulator accumulates all costs resulting from sequential de- 
cisions. All pertinent costs are recorded during the simulation to 
provide guidelines on where excessive costs may be occuring. These 
costs are calculated for each simulation period as described below: 

(a) DIRECT M ANTJF ACTU BINS COST 

(i) labor (not idle) = sum of (Productive Regular Time Cost + 
Productive Overtime Cost + Set Up Cost) for all work 
stations. 

(ii) Machine (not idle) = sum of (Productive Machine Cost + 

Set Up Cost) for all work stations. 

(iii) Materials Used = sum of (per Unit C:ost x Number of Units 
Used) for all items. Per unit cost is recomputed each week 
for each stock mmber end is determined by the ratio A/B, 
vhere, A is the total cost of labour, machines, setups, 
materials, carrying, and ordering since the beginning of 
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the simulation, and B is the total number of units manufac- 
tured or purchased since the beginning of the simulation. 

Thus, per unit cost is a cumulative averse. 

(iv) Total Cost in a Period = (i) + (ii) + (iii) 

(v) Total Cost to Bate = Cumulative sum of (iv) over all periods. 

(b) TOTAL MJdrjPACTUBDTG COSTS 

(vi) Total Labor = (i) -f Idle labour Cost. 

(v±i) Total Machine = (ii) + Idle Machine Costs. 

(viii) Materials Used = same as (iii) 

(ix) Oveihead = Costmt + (Factor of Variable Overhead) x (vi) 

(x) Shift Change Cost = Humber of Shift Charges x Ccst per Ghange. 

(xi) Total Cost in a Period = (vi) + (vii) + (viii) + (ix) + (x) 

(xii) Total Cost to Date = Cumulative sum of (xi) over all periods. 

(c) INVENTORY COSTS 

(xiii) Order Cost = Sum of Purchase Order Costs. 

(xiv) Carr; 7 ' Cost = Sum of Carrying 'Cost for all stock numbers. 

(xv) Out of Stock Cost = Sum of Stock Out Costs. 

(xvi) Total Cost in a Period = (xiii) + (xiv) + (xv) 

(xvii) Total Cost to Date = Cumulative sum of (xvi) over all periods. 

For each simulation period, th<- simulator summarizes all 
relevant aid accumulated costs in the rbm of n Cost Summary Report". 
This report is useful for the PR0SIL1 V instructor to evaluate the 
performance of the user's control system. 
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gost solziahi report 

(xviii) Total Plant Cost in a Period - (xi) + (xvi) 

(xix) Total Plant Cost to Date = Cumulative sum of (xviii) over 
all periods. 

(xx) Value of Materials Received in a Period = Value of all 
purchased items received in a period. 

(xxL) Current Total Value of Inventory = The monetary value 
of the current inventory for all stock numbers. 

(xxii) Penalty and External Costs in a Period = The sum of all 
penalties and other costs incurred external to the simu- 
lator. In each, period the user car impose penalties 
and other external costs. A penalty parameter for excess 
inventory viien multiplied by the difference between 
current value of inventory and maximum allowable inventory 
(parameter) gives the inventory penalty. 

(xxiii) Total Penalty aid External Costs to Date = Cumulative 
sum of (xxii) over all periods. 

(xxiv) Total Score ro Date = (xix') + (xxiii). By this figure 
the performance of the control system is judged. 

PROSE! V is not particular about the accuracy of estimation 
of cost - terms. However, the measure of cost values should be con- 
sistent for various cost — terms. Therefore, it is sufficient to have 
proportional values for cost parameters rather than actual ones. The 
values of various cost parameters should be such that they cause the 
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cost model to be responsive to !t good 1T system designs*. Deteimining a 
good balance of the costs involved in a particular problem is largely 
a matter of v cst and try 11 . 

2. 1.2*5 Simulation Evaluation 

The user receives feedback in PS03BI V in the form of a set 
of reports for eacn simulation period? identified and described as 
below : 

(a) Results of lb recast - This report compares the forecasted 
demand to ihe actual demand on a weekly and cumulative basis* 

(b) Status of Production System - This report shows ihe production 
orders waiting to be processed at each work station at the 
end of the simlation. 

(c) Idle Time - This report shews the amount and cost of idle 
time at each work station. 

t 

(d) Inventory Status - This report shows all inventory transac- 
tions (receipts, issues, on order, back ordered, carrying 
cost, stock on hand, etc.) for each stock number. 

(e) Total Kanufacturing Costs - This report shows manufacturing 
costs for la’ooui-, machines, materials, oveiiiead, and shift 
changes. 

(f) Cost Summary Heport - This report summarizes all costs for 
the period and the cumulative costs for the simulation run. 

The user uses the above information to make any necessaiy 
adjustments to the parameters of h is control system. He then makes a 
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new set of decisions end the entire process is repeated, until the 
reports show a fairly stable control systan, i.e ., how well it reacts 
to disturbaices. 

PBDSIM V has been successfully used by Mize. et. aL. ( fr) 
and others as a very effective teaching aid. They have provided uhe 
Initialization data for various problem like the Shirt Problem, 

Bicycle Problem, Wooden - Case Problem, Iraisistor Radio Problem, etc. 
(IS); to enable the user to test his own decisions for these systems. 

As compared to the earlier PEDSIM model of Greenlaw and 
Hottenstein ('7), the later model - PBDSIM Y of Mize, et. al. ('3)? 
is much more flexible and realistic. Though, PBOSILI covers a wider 
area of a production system than that of PBDSIM Y, it requires the 
decision mafcer to make a large number of inter- related complex deci- 
sions and requires a number of calculations to be made. In PBDSIM Y 
many subobjectives of PBDSIM like Quality Control, Machine Operator 
Proficiency, Plant Maintenance, etc., are overlooked. But it applies 
better mathematical techniques to the limited field. It requires the 
user to make very few decisions and calculations as compared to PBDSIM. 
These two are found to be very effective reaching aids in production 
management, due to the following advantages over the traditional approar- 
ches : 

1. The design of a control system is emphasized. 

2. They force the decision maker to consider the dynamic nature of 


the production environment 
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3. They provide a conceptual understanding of the total operations 
control system and of the interactions among system concepts. 

4. They foster an appreciation of the concepts of feedback, correc- 
tive action and integrated Information systems. 

2.2 Simulation Appli c ations To Real L i fe Production Systems 

So far, very few efforts have been made to apply computer 
simulation in real life production systems. le n-raade (2$) developed 
a partial simulation model for very short-term corrective action using 
actual operating data. Is was pointed out earlier, most of the work 
has been directed towards applying educational simulators to real life 
situations. Fey (9) applied DI2I.AM0 - a simulation language for con- 
tinuous system, to the customs.!.— prodr cer-employmoit model of a manufac- 
turer cf high-quality electronic components. Hie development of the 
model and results are extensively discussed by Forrester (6). 

The results obtained by using DYFJKO are quite useful for the 
manner to stabilize employment levels, factory backlogs and cash 
requirements. However, ix only provides him vdth the overall guidelines 
to make such decisions. On the contrary the PBOSIM simulators help the 
manager to make both short-term and long-term decisions. PBOSIM V 
is ideal for reaching and research purposes. The simulator provides 
the user with the derailed analysis cf his production system for the 
day - to - day planning under situations when he is faced with wide 
range of decision alternatives to choose from. 
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Despite the powerful capabilities of the educational sirxi— 
lators, heretofore little effort has been put to make use of this 
methodology, to solve industrial problems, both in Indie and abroad* 
This can be attributed to the following reasons s 

1. Cost of carrying out simulation experiments is generally 
prohibitive* 

2. Ho so of the modern production systems are highly complex 
in nature. Ihere is multiplicity of highly sophisticated 
transfer functions in the system which restrict the de- 
sign of integrated production planning and control 
system* Ihe term transfer function is used here in a 
broader sense to represent all decision processes, 
mathematical or otherwise* 

3* Lack of accurate and timely information • In majority 
of industries 9 the relevant data is not available in 
the required form, specially the past data and various 
standards of performance e.g. processing and setup 
times, Inventory costs, etc* 

4. Uncertainty prevails in most of the industries* Ig- 
no ranee is a virtue in such Industries. Hany a times, 
decisions made - have Immediate repercussions* This 
again is attributed to the lack of proper data and 
improper mmagement information systems* 
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5. The speed sad the size of computers put a restriction 
on the size of the problem that can be handled in 
computer simulation. Sometimes non-availability of 
high-speed and large-size computers becomes a big 
hurdle in implementing the computer simulation tech- 
niques. 

Currently, the compute r facilities available in this count ly 
is inadequate to cope up with the demand of the industries. Only a 
limited number of Indian industries are able to use simulators for 
the design of their control systems, firstly becaise of very few 
numbers of ’Computer Centres' all over the country aid secondly be- 
cause of the improper utilization of the existing facilities. But 
this difficulty, faced by the Indian industries, is expected to be 
removed very soon if the rate of growth of Computer Technology is 
kept at a more accelerated pace than at present. 

In the presented work an effort has been made to improve 
upon the cap abilities of BEDS Hi V, so that one can handle real life 
industrial problems. The model of Mize, et. al. ('3>) has been con- 
siderably modified and nev< features have been added. The effective- 
ness of the simulator can be best judged, by comparing the results 
of the simulation with those of real life production system. The 
data, from J.K. Electronics for the manufacture of TV Sets, was 
gathered to check the validity of the modified PEOSII! V simulator. 

The model has been developed keeping, the situations in mind -riaich 
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are more consistent with the Indian environment. It should at best 
satisfy the requirements of a typical Indian production manager, vho 
Hay get entangled with several production man^ement problems of 
varied nature due to limited theoretical background. This simulator 
holds a great premise in fhture, specially for Indian industries, 
due to its obvious advantages and because of its wide applicability 
in the fields concerned. 



CHAPTER III 


DEYELOPLJMT OP A MODIFIED PROSB1 V SIMULATOR 

PRO SSI Y simulator as developed by Mize et. al. ( 3 ) 
can. not be usee 5 In its present fonm for a real life production 
situation. This is so due to the following reasons : 

1* PROSE! V was primarily developed as an educational tool. The 
simulator serves more like a 11 tester 11 for the evaluation of the 
user T s decisions, then to aid the user to arrive at those deci- 
sions to control the production sy st em . 

2. It does not incorporate some of the very important features of 
the production systems e.g. s sales forecast, real estimates of 
variables and parameters for inventory control, etc. 

In order to account for some of the above-mentioned limi- 
tations, PRGSIM Y has been modified considerably in the present re- 
search. The modified simulator is referred hereafter as n MOD - 
PRDSIM V 1 *. The modifications incorporated in the simulator car be 
broadly categorized as follows : 

1* A general purpose m odel has been developed to forecast the sales 
of finished goods depending on the past demand history of the 
particular product. 

2. The inventory control model used in MOD-PRO SIF Y incorporates 
features for determination of the values of pertinent inventory 
parameters for each purchased item e.g. , part requirements of 
the finished good, economic and discount order quantities of 
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the purchased items, etc. These paraneters are important, to 
arrive at an optimum inventory decision policy. 

3. Tor effective production planning, a decision policy has been 
formulated for generating production orders for each manufac- 
tured item Airing each simulated period. 

Above mentioned modifications have been explained in de- 
tails in the following sections. To suit the Indian environment, 
MOD-PBOSIM V takes into account the six working days a week. While 
in the original model of Hize et. al. (3), it was assumed that the 
production system works for five d^s a week, as explained in 
Chapter II. 

3.1 Development of a Generalized forecaster flodel 

The simulator determines the forecasted demand for each 
finished good using a time series model, viiich Is capable of gene- 
rating demand trends having a wide variety of slopes, growth rates, 
cycles and variations. A tine series can generally be considered 
to consist of a random component, an arbitrary component such as 
catastrophic events and non-random components that include ti’end, 
cyclicity and persistence. The tern random component refers to 
the pure random component of the process. This component is gene- 
rally very difficult to assess and is usually neglected in stable 
systems. The arbitrary or catastrophic events are rare events about 
which no prediction can be made because of lack of data and their 


b 
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arbitrary behaviour. The non-random component may itself be composed 
of trend and cyclicity which is a function of time alone; ‘and per- 
sistence viiich defines the sequential dependence of the time series. 

The tine series used in the model for generating the fore- 
casted demand for each finished product is represented by the follow- 
ing mathematical relationship : 

¥ (t) = A + Bt + Ct 2 + D (Cos ( -n t/Ej) + P (Cos ( ?St/&)) 

+ h (cos ( n t/u) + dev (3.1) 

where 

w(t) = generated demand of the finished product for the time 
period "t n , in units, 

A = the beginning constant or n Y” axis intercept, 

B = coefficient of "straight line" slope, 

C = coefficient of quadratic growth term, 

D, E, E = coefficients which determine the amplitudes of first, 

second and third cyclic trends respectively, 

E, G, U = coefficients which determine the periods of first, 

second aid third cyclic trends respectively, 

DEV = a randomly chosen deviate from a normal population 

having a mean of zero and a standard deviation of AK. 

This normal distribution is truncated at three standard 
deviations below and above the mean. 
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This forecaster generates the demand f or each period ass umin g 
a linear-quadratic growth. As many as three Types of cycles can be 
imposed on a single demand trend. The arbitrary component of "the time 
series is represented by the variable DEV, the value of which, is cal- 
culated by using Monte Carlo Simulation technique (21 ). This keeps 
the usual off-control points in the frequency analysis of the demand 
pattern at the confidence level of 99.73$ by using the control limits 
of three standard deviations. Using this model a value for each 
parameter (a, B, C ..... etc.) can be calculated for each finished 
good. The demand trends for the several finished goods can be widely 
different. As explained in Chapter 171 , regression analysis is used 

for the estimation of the coefficients (a, B, C etc.), as 

sufficient past demand data are available. 

3.1.1 Regression .Analysis 

Regression analysis deals with the determination of the 
causal relationship between the exogenous and indogenous variables. 

The demand coefficients of regression equation are generally estimated 
by the method of least squares (22). To calculate the values of the 
demand parameters, the values of the frequencies of period cycles 
must be known. To estimate the values of coefficients E, G- and U; 
iterative procedure ana ’’Correlation aid Spectral Analysis” tech- 
niques axe used. The techniques are discussed in the later sections 
of this chapter. Once the periods of the cycles are identified, the 
values of other coefficients is calculated in a manner given below : 
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1. ‘The variables w(t) aid t axe redefined, tc reduce the computa- 
tional efforts, as follows : 

Z(t) = W(t) - ¥ 1 

\ (3.2) 

and h = t - n/2 J 

There 

n = number of available past demand data 
W = average of past demands over last n periods. 

2. Values of the coefficients in equation (3.l) after above substi- 
tution of Step 1, is detemined by solving the determinant having 
elements in first row as K(t), 1, h, h , Cos ( 7tt/E), Cos ( JX t/G) 
and Cos ( 7T t /Cf). These are the multiplicants of coefficients 

(A, 33, C .... etc.) in equation (3.l). The elements in other 

rows of the determinant are calculated by summing the product of 

first row elements with the first row element of the particular 

column over all the past n periods. To find the value of ith 

element of jth rov/ of the determinant, summation should be taken 

of product of ith and jth elements of first row. For example, 

- — 2 

the value of 4th element in the 6th xow is given by - j!— ^ 

n 

Cos ( TTfc/G). The order of this determinant is 7, as the number 
of degrees of freedom of the demand variable is 6. 

3. The value of this determinant is set equal to zero* The elemental 
determinants of this determinant are of the order 6. The values 
of these elemental determinants - B1, Cl, Dl, El, PI, 01, and 

F2; are calculated by employing Pivotal Condensation Method (23) * 
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3 y setting the value of the determinant equal to zero , the follow- 
ing equation is obtained: 

K(t) x B1 - Cl + h x 351 - h x El t Cos ( Ft/s) x PI 
- Cos ( 7Tt/&) x 01 + Cos ( TXt/t 7 ) x 12 = 0 (3.3) 


Substituting the values of E(t) snd h from. equation (3.2) into 
equation (3.3) and comparing the coefficients of equation (3.l), 
the values of the demand parameters axe obtained as follows ; 


3 

C 

P 


(Cl + PI x n/2 + El z (n/2)) / B1 + W 
(- PI - El x n) / El , 

El/Bi s P = - Pi/31 , 


(3.4) 


01/BI 


= - P2/B1 


4. The value of AK, standard deviation of demand, is calculated from 
the past data using the equation : 


AK = 


S-) (Er 
n 


xor m xorecasr 


y 


n - f 


vdaen. n f 


(3.5) 


#ie re 9 

f = degrees of freedom (6 in this case) 
and Error in forecast = difference between actual and forecasted 

demand . 


To use the above mentioned method a minimum of 7 periods 
demand history must be known because the degrees of freedom is 6. 
The accuracy in the estimation of the demand parameters improves 
with the increase in the amount of past data available for their 
estimation* 
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3.1.2 Deteim ination of Periodicity in the Cycles 

Tv/ice the values of coefficients E, G and U in equation 
(3.l) represent the periodicity in the three cycles wnich are pre- 
sent in the tine series, To estimate the frequencies of the cycles, 
the following procedure is followed : 

1. The optimum values of E, G ana U are estimated for the a inimum 
total absolute error in forecasting in the past n periods. The 
demand variable in equation (3.l) is assumed to be represented by 
a discrete time series, hence the periods of cycles are integers. 
This is true in many real life situations. By employing an ite- 
rative procedure the various values of E, G and U are examined 
for getting the minimum value of the total absolute forecasting 
error for past n periods. The values of coefficients E, G and U 
are varied in the steps of 1/2 in the range from 1 to n/2. This 
is done because the periodicity of cycles have integer values and 
the minimum possible number of periods in a cycle is 2. This 
method is useful specially when the value of n is not very large. 
The number of iterations increases considerably with the increase 
in the number of past demand periods. 

2. The calculated values of E, G and o are checked by correlation 
and spectral analysis. The presence of periodic cycles in the 
time series is indicated by the periodicity in the correlogram 
and spikes at appropriate frequencies in the spectral density 
function. These techniques have been used to calculate the 
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approximate va3u.es of E, 5 and. U for very large values of past 
data situations -when the iterative schemes are very difficult to 
use. Once the approximate values are known, the exact values in. 
their neighbourhood are found by iterative method as discussed 
above. The details about correlation and spectral analysis tech- 
niques are given in ippendix A. 

Depending on the type of product, the time series can be 
easily modified to take into account different types of trends and 
periodic components. A trend represents the long - term variation 
of a process. A plot of past data indicates whether a trend is pre- 
sent and if so, -whether it can be expressed as a linear, polynomial 
or exponential function of time. The actual nature of the trend may 
be difficult to identify because of the presence of cyclic and random 
components. To eliminate these random components, the method of 
moving averages (l) is used. Ihen n values of past data is available 
the maximum possible number of periodic components in a time series 
is n/2. Generally a much smaller number of significant harmonics 
are considered. The techniques of correlation and spectre.1 analysis 
are used, when the exact number of cycles for any product are not 
known. 

Since the time series of demand variable of any product 
may change with time, necessary corrective action has been taken 
after each demand period to modify the forecaster rrndel. Hiring 
each period the forecasted demand is compared with the actual demand 


to evaluate the effectiveness of the forecaster. 
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3.2 Inventory Control System 

The purpose of the inventory planning and control function 
is to determine the appropriate inventory policies with a view to 
keep the various costs associated 7.1 th the inventories at a minimum. 
The most common classification of inventories is - (i) raw materials 
and purchased parts, (ii) in-process inventories and (iii) finished 
goods. Star each purchased item the decision is to be made dout how- 
much to order at a time aid when (or how often) to place an order. 

In M0D-5EDS3H V the purchase orders are generated for the 
raw materials required for the production of finished products during 
the simulated period. The order quantity of each purchased item is 
determined from the forecasted demand of the finished goods during 
that period. The purchase orders for each raw material are accumu- 
lated in the simulator, o-ver the simulation run, till the Cumulative 
Order Quantity (COQ) of the particular item readies the value of 
its Economic Order Quantity (EOQ) . EOQ of any purchased item is the 
optimal ordering lot size of the item viaich minimizes the total in- 
ventory costs associated with the item. During the simulation period 
if the EOQ of any purchased item is exceeded by its COQ a purchase 
order equal to COQ is placed. 

In the above stated inventory policy it is assumed that - 
(i) sufficient safety stocks for us- during the variable lead times 
are maintained aid (ii) the forecasted demand during the lead time 
does not vary much i.e. , there are no unexpected fluctuations in 
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trie demand trends. These assumptions ere valid for most of the 
production systems. 



Quantity Ordered (q) 

FIGURE 5 

TOTAL COST FU2JCTI0IT FOR CLASSIC AL IKTOTTORY MODEL 

To implement this inventory policy the value of BOQ for 
each purchased item must be precisely known. The effect of order 
quantities on the various costs is represented graphically in 
Figure 5. The 'EOQ for each item is determined by using the Cla- 
ssical Inventory Model. The economic order quantity can be cal- 
culated using the following relationship : 
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EGQ = J units (3.6) 

where 

PC = procurement cost of each order. 

CC = carrying cost per unit per tine period, 
and D = maximum possible usage rate of the item per tine 
period. 

The values of parameters PG ana CC for each purchased part 
are known ffcr the production system to be simulated. STOP PROSE 7 
computes the values of D for each purchased item by multiplying the 
parts requirements for each unit of the finished goods (e) with the 
Estimated Plant Capacity (EPC). The value of EPC for each finished 
product is calculated by dividing the total available time per simu- 
lation period by the cycle time or the longest process time of the 
item. The total requirements of purchased items for each finished 
item (r) must be known precisely to the simulator. This is easily 
done by making a few calculations for a small problem. Hit, litis 
procedure tends to be too clumsy for larger problems involving a 
large number of parts and sub-assemblies. In the literature the 
problem of the determination of total parts is referred as the 
"Explosion Problem". Very many techniques based upon the matrix 
algebra have been reported in the literature. IIOD-PBOSB! V 
employs the algorithm of "Plug and Chug method" for the explosion 
problem. This is a much more simplified technique as has been 
explained in details in Appendix B. 



47 


For each purchased item the value of EOQ is calculated by 
using the equation (3*6)* In reel production systems a raw material 
supplier often offers a quantity discount fiien orders loiter than an 
order size are placed* MOD— F 30 SIM V takes this feature into account 
by assuming that only one discount price foi each material is pro- 
vided by the supplier* This discount price is applicable only for 
certain order size called n Discount Order Quantity (D0Q) u * The 
value of DOQ for each item is computed by multiplying a ^Factor of 
Discount Order Quantity (FD0Q) !! with the corresponding value of EOQ* 
The value of FDOQ is assumed to be same for all the purchased items* 

.Another import sit variable used in the inventory control 
system is the I! lead time 11 of each purcha.se order. As explained 
in Chapter II , MOD-PBOSET V generates the value of lead time using 
Monte Carlo Simulation. 

3*3 Operations PI arming 

MOD - PEOSDI V incorporates all the above mentioned fea- 
tures to structure a decision policy for raw materials inventory. 

As explained In Chapter II the in-process Inventory Is governed 
by the hold quantities at each work station. To control the fini- 
shed goods inventory a decision policy Is formulated for the 
production planning of all assemblies and sub-assemblies* This 
requires the generation of production order for all the Items 
manufactured during each of the simulated period. Weekly produc- 
tion orders are generated from the iorecasted demand for each. 
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finished product. This assumes that the production orders placed 
in a week satisfy lie weekly demand of the finished products. lUr- 
ther 5 it Is assumed that the forecasted demand is a goGd estimate of 
the actual demand and the planned production can be achieved in most 
of the simulation periods. 

Hos-g of the assumptions made in 'ICD-PBOSEi 7 hold good 
for real production systems. The simulator reduces the load of the 
production manner considerably and provides him with a decision 
policy which he can follow/ on a long term basis. 
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COMPUTER PROGRAM DETAILS 

Two separate packages of computer programs written in 
FORERAN IV have been developed for M0B-PRDSE1 V. These two are 
complementary to each other in use and solve the mathematical 
models proposed In earlier chapters. One ox the programs is for the 
forecaster model. It estimates the various demand coefficients for 
the simulator. The other program is the simulator itself. It uses 
the output of the earlier program as the initialization data. The 
two programs are kept independent of each other for use because this 
simplifies the handling of the data and reduces the memory storage 
requirements of the computer. The details of data requirements of 
data collection are given in Chapter V. The FORTRAN listings of these 
programs have been included in the Appendix 0. Listings contain a 
generous number of COMMENT Cards to assist the user in following the 
detailed logic. A list of various parameters and variables used in 
the program is also included in the listing of each program. 

4.1 Forecaster Program 

This program has been developed for the generalized fore- 
caster model. It estimates the values of forecaster parameters for 
the finished good, as was explained in the previous chapter. The 
package is used independently for each of the final oroducts. It 
provides op timum values of the parameters in the fom of printed re- 
ports and punched cards. The punched cards axe used directly as 
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simulator initialization data input, The size of the problem that 
the program can handle is defined by the number of available past 
demand data. The program has been designed for a maximum of 100 
past demand periods, but this restriction can be relaxed by modifying 
the DIMENSION statements. 

for 25 past demand data, the iterative procedure requires 
127 seconds of computer time on IBS 7044 system. Bit, vhen the fre- 
quency of each demand cycle is precisely known the computation time is 
reduced considerably, 'or example, the sane problem took only 6 
seconds vihen the demand cycle was precisely known. The program con- 
tains one MAIN program aid one subroutine FORCS, which is called from 
the MAIN program for each iteration. The gross logic of these programs 
is explained in Figure 6 aid Figure 7. 

4.2 Simulator Prog ran 

The computer package consists cf approximately 1200 statements 
in FORTRAN IY. A complete report of the results is printed for each 
week of simulation. Two sets of punched outputs are obtained during 
the simulation run. These punched cards are used as input data for 
further simulation runs. The simulator consists of a MAIN program 
and seven subroutines. The functions of each of them are explained 
in the following sections. The general logic incorporated into the 
simulator is given in Figures 8 to 12. 
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READ : 

Number of past demand periods ( n) aid 


the demand history. 

IS SIIN ; 

Number of degrees of freedom (ill). 

INITIALIZE ; 

Minimum total absolute error 


(MINAS! 1 = 0); Number of periods in 


cycles - (NIT , III and Ll) = 2. 


NIT = NIT +1 
= Mil +1 
£ NIT ■ 
= EL +1 
^ ML 



CALL FORCS 
(hit, LEI, ll) 


ffiUTBiE f 
JDTJfflKHL' 


MDTA3D = IA3TCT 
ENOPI = ra 
MMOPT = HM 
LLOPT = Ii 


kutabe .= o 



Iterations 
' to calculate 
the optimum 
• number of 
periods in 
3 cycles 

(mopt, maps 

and HOP!) 

: for mlninHam 
■ value of the 
total absolute 
error (lAHPOT) 












FIGURE 6 


PLOW DIJGBil! OP VJJS PROGRJ1I OP THE FORECASTER, 
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INITIALIZE 

: Elements of Matrix 'A' of order 7. 

CALCULATE 

: “ ' i ~f of Forecaster Parameters ( E , G 


aid U). Sum of past demands and 


Average of past demands. 


J + 1 



J = 


is an index. 


CALCVIiAEE : "values of Oeneats of 
lower triangular ' A' 
Matrix. 

ST I : New values with past 

values. 




SET : Values of elements of upper triangular 

' A* Matrix by symmetry. 

CALCULATE : Values of 7 Determinants by eliminating 
fist row aad each column, in turn, of 
’A 1 matrix using Pivotal Condensation 
Method and ’’clues of other Forecaster 
Parameters (a, 3, C, D, I a_d H). 

INITIALISE : Total Forecasted Demand, Total Actual 
and 'Absolute Errors in Forecasting. 


I is an index. 



(Figure continue d) 
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IIGTJRE 7 


PLOT,’ VIjaiUZl TOE 


SuBROUT HIE EOECS OR THE EOEECi-^SIEE 
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(l;gure continued) 







puitch 


: Production status values for work 
stations; stock on hand end stock 
on order for each, stock number. 

WEI IE : "SIHJLiTriOH COMPLETED” • 

BEAD : Actual unit cost sod actual in- 

ventory cost. 

WHITE : Table of summary of results 

(SDISET) . 

BE/D : Plot characters (PIT ) for 5 

plots. 

CALCULATE t Scaled values of SUMEP for 5 
plots = T (0 < T < 50). 


WEITE : Values of 'T' on tine scale for 

5 plots. 



PICT HE 8 


PLOW DIALED! POE MAIK PBQGH/JI OP THE SIMULATOR 
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CAIOULAIB 

s longest process tine (iPT) from 
all process times. 

WRITE 

: IPT. 

OiliCULiiCE 

: Estimated plant capacity (E3?c). 

EBAI) 

: Factor of discount order 
qu entity (FDO Q ) . 

WRITE 

: EPG and FBOQ. 



FIGURE 9 

FLOW Dli&RAM FOR SUBROUTINE MAPA OF THE SIMULATOR 



HEAD : 

Forecasts (actual sales data). 

G-E'sHRAiTE z 

Actual demands for this period. 

: COMPARE : 

Forecasts to actual. 

WRITE : 

i . 

Results of forecasts. 


CALCUIiAIE : Purchase orders required for 
raw materials in this period. 

’ GENERATE : Random lead time for each pur- 
chase order. 

CALCULATE J Bay of arrival of each purchase 
order. 

WRITE : Purchase order numherj stock num- 
ber and order size for each item. 



FIGURE 10 


FliOW DIGRAM FOR SUBROUTINE FORCST OF TEE SIMULATOR 




STAKE MASTER CLOCK 



(Figure continued) 
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4 * — — — 

Simulates one TIITC of ope ret ion for Work Station I. I 

Depending upon the set of conditions at tine 1, simulator will 
perform one or more of the following functions on Work Station 
(W.o.) I : 

- Continue processing the lot now on the W.5. 

- Move the lot to the 'Held Block’, if processing is complete. 

- Move the 'Hold Quantity’ in the ’Hold Block 1 to the first availa- 
ble ’Queue’ position at the following W.S.; or , move the ’Hold 
Quantity’ to inventory (vhich increases stock on hard) if W.S.I is 
last on the network*. 

- If ’Queue Position’ 1 is empty, the ’Queue' is indexed $ i.e. , all 
orders in the ’Queue’ are moved up one position. 

- If there are no orders in the 'Queue' , Hie machine is idled. 

- If an order is in ’Queue Position’ 1, stock on hand of all items 
added at W.S.I is checked to be sure, sufficient quantities exist 
to produce one lot. 

- If insufficient materials exist, the ’Queue 1 is searched for an 
order that can be produced. This order is moved to 'Queue Posi- 
tion' 1. 

- Whenever there is an order in 'Queue Position’ 1 for vhich sufficient 
materials exist, one lot is moved to the W.S. Set-up time, if any, 

is calculated. Processing time is determined. This determines the 
time vhen the lot will be finished. 

- Various costs are accumulated. 

- Inventory levels are updated. 

- Idle time (if any) is recorded. 

- Humber of set-ups are recorded. 

- — " * ~ ' " " ~ ~ 

t 



(figure continued) 
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)J5u is Ihe number of 
jy.ork stations in the 
I problem 



Fill = maximum time available; 

no more processing 
■ after this time on day 

I JR5. 


One dgy of simulated production has been completed. 

- Subtract 1/6 of weekly demand from stock on hand of 
each finished good. 

- Compute carrying cost for each stock number for day 
JR5. 

- Compute out-of-stock cost for day JR5. 

- Determine if any purchase orders arrived on day JR5; 

if so, add to stock on hand. ! 


Z' N _ ./JR5 h 



RESUM 
\ / 


FIGURE 11 

FLOU DIAGRAM FOR SUBBOUTIUE PROD IE OF THE SIMULATOR 
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1 CjiliCULASE 

? 

I 

WRITE 

TRITE 


WRITE 


Cumulative -unit cost for each 
manufactured item. 


Repor- "htveietcry status". 

Cost Reports - "DIRECT 

TiATUFACiuriuc- coses, total 

L' ART PACT TJBIU 0 00 SI F 5 
ZUTSLII CRY COSTS." 

"COST S'JlILC-v REPORT" . 



] RECAP is a subroutine 
j that collects infoi*- 
| nation for evaluation 
I in KAQT PRDGRALI 

j (sumry) . 


FIGURE 12 


FLOW DIAGRAM ID R SUBROUT HIE IITVEU1 OR THE SIMULATOR 
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4.2.1 LLAIIT Program 

This program consists of three parts. The first part 
is used to create the data set for the initialization of the simu- 
lator. It also calls the subroutine MPA to generate the inventory 
parameters. The second part calls the subroutines PCRCST, PBDKDIfl 
and ETVETT during each period of simulation. The third part of the 
program prints tables and plots charts for the results of the simu- 
lation run. It is used after the simulation has been carried out 
for the desired number of periods. The relationship between perti- 
nent variables is plotted for each of the finished products. The 
plot scales can be suitably modified to meet the individuals require- 
ments. 

4.2.2 Subroutine Programs 

Mainly three subroutines govern the simulation of the 
production system, for each simulated period. These are identified 
as IORCST ? PROEM and IRVHJT. Subroutine IORCST computes the fore- 
casted demand for finished products. It also generates the lead 
time for each purchase order. Subroutine PROBE! generates the 
production orders for the manufactured parts. It uses a ’’Master 
Clock” to simulate the production activities at each work station 
for each time increment. It also keeps the timely ana cumulative 
record of the production times, the inventory status and the 
various costs involved in the production system, "or each simu- 
lated period the subroutine WRITES prints intermediate and end-of- 
the-week status of the production system. 
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Subroutine HTYEIT calculates sad prints the various 
weekly and cumulative costs. It calls the subroutine BECAP to 
store all important results obtained in the simulated period. These 
pertinent results are used at the end of the simulation run in the 
main program to print tables and to plot the graphs. Subroutine 
MAPA generates Ihe inventory paraaetars for each purchase item, 
iflhich. are required in the simulator. As explained in Appendix B, 
the punched output of the ''Explosion Program for unit requirements 
is used as the input data for MAPA. Subroutine RAIffiTJM generates 
two random numbers for the given seed values, having normal dis- 
tribution and values betv/een zero and one. These random numbers axe 
used in Monte Carlo sampling for generating demand aid lead time 
distributions in the subroutine EOBCSI. 

4.3 Capabilities of MOP-PBQSIM Y 

MOD-PBOSIM Ths been basically developed for the 
IBM 7044 system. The size of the available computer memory puts a 
constraint on the capability of the simulator* MCr-PBOSIM Y is 
capable of handling a manufacturing system consisting of upto 15 
work stations, 60 stock numbers and 3 finished goods. This requires 
about 100,000 bits of computer memory. To accommodate larger pro- 
duction systems the procedure of "cut and try” is employed, as ex- 
plained in Chapter Y. 

The computer time required for each simulated period 
depends very much on the computational speed of th< computer used. 
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It is also a function of the size of production system, the dura- 
tion of the simulation lun ana the time increment used for simula- 
tion. The IS! 7044 system requires about 70 seconds to simulate 
the 15 x 60 x 3 problem (l5 work stations, 6C stock numbers aid 3 
finidaed goods) and about 16 seconds for 7 z 20 x 1 problem. For 
both the problems the simulation run ras 25 periods (weeks) with 
the time increment of one minute. Che computer time can be reduced 
considerably by using a larger time increment. Maximum possible 
value of time increment is the largest common divisor of all lot pro- 
cessing times and the set-up times. 

In addition to the printed reports of simulation results, 
M0F-PB0SIM V also provides a deck of punched cards as output. One 
set of punched cards consists of the initial Inventory status of 
the production system. Other set of punched output consists of 
stock-on-hand aid stock-on-order for each purchased item for the 
last simulation period. These o utputs can be used directly either 
in other computer programs or for further simulation run. In 
MOIKPEOSIM 7 the computer run is administrated and the tape/disk 
systems are used m a way similar to that done by Mize et. al. (3)» 
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DAI A EE^riHS"3 A17D DAI A COLLECTIOR 

In order to simulate a production sister on ItlCD-PBOSHI 7, 
it is necessary to structure a model for the production system under 
stud;/-. An extensive study of the production system is made by ana- 
lyzing the various operations mvol“3d sod the pertinent parameters. 

Since MOD-EBOSIII 7 is a general purpose simulator, many 
different types of production problems can be simulated on it. She 
major difficulty in problem formulation is faced ?hile structuring 
the system's model suitable to the requirements of the simulator. 

She following restrictions are placed by the simulator on tire problem 
representing a production system : 

1. MOD-P5DSE5 7 is restricted in its application to large complex 
systems due to the computer limitations. Eence the problem is 
structured using as few work stations and stock numbers as 
possible, without seriously affecting the characteristics of Idle 
production system under study. In case of large-size systems 
this is achieved by selecting an independent sub-system . for 
viiich relevant; data are available. If possible the purchased 
items are grouped together according to their functional use. 

All this results in a compact system - enabling easy analysis. 

2. Any item carried in the inventory is assigned a unique stock 
number. [Finished products are assigned numbers 1, 2, 3, — n, vhtere n 
is idle total number of finished products. All other items can 
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oe numbered 1 n any desired orders starting with n -r 1. ITo dis- 
continuity In numbering is permissible. 

3. Items requiring processing must have an operations network. The 
network may consist of one or more work stations* Multiple pro- 
cessing of any item at the sane work station is not permissible 
to avoid the looping* 

4-* There is a limitation on the number of components In any sub- 
assembly* In IIOD-PHDSIK V this parameter is specified as J5G* 

The maximum value attributable to 140 is 5. In the operations 
network of a processed item, the full amount of aiy component re- 
quired must enter at one aid only one work station* 

5. I? lliminy assemblies 11 are incorporated to account for situations 
where an assembly requires more than J4C components. 

6* Similarly 5 work stations processing more than J40 components 
are sub-divided into two or more work stations* 

7* Sub-assemblies are assigned the same stock number as the final 
assemblies. 

The following sequence of steps are suggested for structuring 
the model of a production system, keeping in view the above mentioned 
restrictions. This also aids in determining the data necessary for 
simulation. The procedure delineated herein Is in context to the 
TV industry under study. 
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TABLE 1 

IlTTETl'OHY DESCPJPTIOi: 
SET 

1. If Stage - finished Pioanct 

2. Printed Circuits 


3. 

Hesistors 

35 

different 

values. 

4. 

Ccndensors 

26 

different 

values. 

5. 

Transistors - 

6 

different 

values. 

6. 

Diodes - 

3 

different 

values. 

7. 

Peaking Coils - 

3 

different 

values. 

8. 

Coils VIP 

6 

different 

values. 

9. 

Coils SIP 




10. 

Trap Coils 




11. 

Disc rimin ato r 




12. 

Detector Can 




13. 

Test pin 





14. Wire Jumpers - 3 different values. 

15. Solder-Consumable Item 

16. Potentiometer. 

17. Pirst Dummy Assembly. 

18. Second Dummy Assembly. 

19. Third Dummy Assembly. 

20. Iburth Dummy Assembly. 
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low, sz.'zio:: 

LIZOUT 


_ bw .. 


W.S. 1 
— 

! 17 ; 


is) 


’;.s. 


@ 

©L 

W.S. 3 



F 


Y * b* b 


W.S. 7 


* 

© 


DESCRIPTION 


Solder 35 Ho sis tors 


Solder 25 Cordensors 8 3 Wire Jumpers 
and the Potentiometer. 


Solder 6 -Transistors, 3 Modes, 
3 Peaking Coils and 7 Test Pins 


Solder 2 SIP Coils, 3 Trap Coils 
Discriminator and Detector Can- 


Solder 5 PIP Coils 


Testing Sound Section and Repairs 
if required. 


Testing TIP Section, Tuning, 
iligrnent ana necessary repairs 


FIGURE 14 
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1* For production systems, vice re the number of to ik stations and stock 
numbers run into hundreds , analysis is carried out only for a par- 
ticular product out of the several products bsir r manufacture^* 
Thus, in the present case, a particular sub-assembly called n IF 
STAGE 11 is selected for detailed analysis. This particular pro- 
duct was selected as the management considered it to be a bottle- 
neck in the final assembly of T T ~ sets* In addition, data for 
l? IE STAGE** was easily available* 

2* All inventory items were listed end unique stock numbers were 
assigned* Preferably, the stock numbers are to assigned that 
larger stock numbers go info the assembly of smaller stock num- 
bers. This avoids renumbering of stock numbers for !f Explosion 
Program 11 as explained In Appendix B. Items list for the predic- 
tion of U IE STAG-E ! ‘ is given in Table 1* The stock numbers are 
reduced by grouping different items together as shown in Table 1 
i # e. resistors, condensers, transistors etc. of different values* 
Grouping can be done only for items requiring similar operations, 
seme operation tine, ordering costs and purchase orders* 

3* The sequence in which the components are used in the production 
system has to be indicated* Specific work areas are established 
on the basis of machines or operators working in that work area. 
Each work area is assigned a work station number- The nature of 
work done md the components used at each work station are indi- 
cated in Figures 13a snd 14* 
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1. For production systems, where the number of voik stations and stock 
numbers run into hundreds , analysis is carried out only for a par- 
ticular product out of the several products terry manufacture^. 
Thus, in the present c&se, a particular sub- assembly called n H? 
STAGE” is selected fer detailed analysis. This particular pro- 
duct was selected as trie maiagement considered it to be a bottle- 
neck in the final assembly of TV sets. In addition, data for 

n IE SC/GE" was easily available, 

2. All inventory items vrere listed aid unique stock numbers were 
assigned. Preferably, the stock numbers are so assigned that 
larger stock numbers go into the assembly of smaller stock num- 
bers. This avoids renumbering of stock numbers for "Explosion 
Program" as explained in Appendix B. Items list for the pro dic- 
tion of "IP ST ICE" is given in Table 1. The stock numbers are 
reduced by grouping different items together as shown in Table 1 
i.e. resistors, condensors, transistors etc. of different values. 
Grouping can. be dcr.e only for items requiring similar operations, 
seme operation time, ordering costs and purchase orders. 

3» The sequence in which tie components are used in the production 
system has to be indicated. Specific work areas are established 
on the basis of machines or operators working in that wcik area. 
Each wozk area is assigned a work station, number. The nature of 
work done aid the components used at each work station are indi- 
cated in Figures 13a aid 14. 
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4. Ill operations are examined to find the maximum possible number 
of earponants (J4C>) required for an assembly. Hinny assemblies 
are created viien ever the number of components required exceeds 
the maximum value of J40. The production of "IF STiGE" requires 
the assembly of 15 items. Therefore, the dummy assemblies labe- 
lled 17, 18, 19 and 20 had to be evolved. 

5. Additional work stations were incorporated on similar grounds. 

The values of processing time, set-up time, in-process inventory, 
man rate and machine rate, of the original work station are divi- 
ded between newly created work stations. For "IF STAGE" the origi- 
nal work station 4 processing 6 items has been divided . into work 
stations 4 and 5 (Figure 13a). 

6. Assembly digram is drawn for each manufactured item (including 
the dummy assemblies) as shown in Figure 13a* For assembly of 
"IF ST/GE", although solder is a consumable item required at five 
different work stations (Figure 14), it is allowed to enter the 
system at work station 1 only but is carried to other work stations 
as and vhen necessary (Figure 13a). 

7. Figure 14 depicts toe work station layout indicating tire operar- 
tions network for each manufactured item and tie dummy assemblies. 
For "IF STAGE" the processed Item at work stations 5, 6 and 7 is 
referred to as stock number 1 although, it has not been processed 
completely at these work stations. 

8. For convenience, work station orders are tabulated as shown in 


Table 2 
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TAELE 2 

TOE K ST ASIDE 0E3EE 
Paraaeter PST (i, J, 4) IC = 1 


WS 


I 

J 

1 

2 

3 4 

5 

6 

7 

1 





1 

2 

3 

17 


1 






18 



1 





19 




1 




20 




i 





TABLE 3 a 

COMPONENT HEQUIHSIEBTS 





Parameter 

EMR 

( I, 

J, K 

) 

10 

= 2 





J 

5 

! 

I 



II 


1 


III 

”1 

(, 

17 



V 


K 

[ 

fHf 

2 j 

P 


2 

l 

jj 

1 

rri 
2 y 

3 1 

1 

1 n I 

l J 

3 

1 1 
1 i 

2 I 

JLi 

3 

1 

20 

1 


8 

6 


5 










17 

2 

1 

1 

15 

1 


1 

3 

35 

1 







18 

17 

1 

2 

4 

26 


2 

16 

1 

2 

14 

3 

2 




19 

18 

1 

3 

5 

6 


3 

6 

3 

3 

7 

3 

3 

13 

*7 

4 

3 

20 

19 

1 

4 

12 

1 


4 

9 

2 

4 

10 

3 

4 

11 

1 

4 


- - 
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9. Table 3a gives raw material requirements as arrived at from 
Figure 13a. For each manufactured item, the table contains as 
many as J40 sets of data. Each set contains the info mat ion m 
the following order for each component : 

(i) Stock number of the component. 

(ii) Quantity of the component required for cne unit or 
batch of the manufactured product. 

(iii) Work station ax which the component enters the pro- 
cess. 

10. Set-up times are determined for manufactured items at each work 
station. As the time increment used in the simulation xs cne 
minute, set-up times are assigned integer values (Table 4). 

11. Processing times per unit or batch (also integers; are determined 
for each product at each work station (Table 5). 

12. For each manufactured item the in-process inventory is determi- 
ned in lots, corresponding to the respective work station. For 
dummy work stations this value is zero (Table 6). 

13. The labour aid machine rates for aLl work stations are tabula- 
ted for each minute of operation (Table 7). Machine costs are 
ascertained for working and idle conditions. For IF Si-GE 
they are not different as the power is not disconnected vhen 


not soldering. 
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TABLE 5 

PBOCESS TEIl/ONIT (OE BATCH) 117 TIIITUTES 

Parameter PST (i, l) IC - 1 
Work Station 


I 

~jT 
L- 1 

1 : 

_1 2 

I T~ 

! 5 i 4 

t~t 
x 5 

“1 7" 

I 6 

t~ 7 

I 7 

1 




0 

-> 

1C 

15 

17 

IS 






18 


17 





19 



13 




20 



10 





TABLE 6 

HOLD QUANTITY H7 LOTS 


Parameter PST (l s J, 3)? It = 1 
Work Station 


I 

ji 1 

n~ 

i 

~~r 

— i — 

3 J 4 

~r 

i 

i 

5 i 


H 

7 

1 






1 

i 


1 

17 

1 









18 



i 







19 




1 







20 


0 
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TABLE 7 

1TA5 £13 ) MACHETES RATES ET 3s ./MINUTE 
Parameter OHR (i, j), IC ■- 3 


j' 15 

1 

l 5 

6 

I I 

Fen 


Machine Idle 

1 

0,C15 

G.0C25 

0.0025 

2 

0.015 

0.0025 

0.0025 

3 

0.015 

0.0025 

0.0025 

4 

0.008 

C • 00 *1 

0.0013 

5 

0.007 

0.0012 

0.0012 

6 

0.020 

0.0025 

0.0025 

7 

0.030 

0.0025 

0.0025 



78 


14. Tlie inventory information for "I? ST AGE" is Tabulated in Table 8. 
The unit costs of manufactured products are computed by summing 
up the unit costs of purchased items repair-.,-! to produce one 
unit of the naauf actured product. Unit costs of purchased com- 
ponents are set ecu . 1 to toe regular purchase price. Carrying 
costs are assumed proportional ho the unit costs. The discount 
prices are provided bp the raw material supplier. Generally 
these values are 5 to 20 percent lower than the regular prices 
( 3 ). For "IF STAGS" this value is 10 percent. The values of 
standard deviation of lead rime are determined from past supply 
data. Generally the value of standard deviation does not exceed 
one third of the average value. Backorder or out of stock costs 
for finidied products are generally 2 to 5 times the costs of 
the units ( 3 ). The value of initial stock on hand depends en- 
tirely on the management policies. This value is estimated by 
multiplying the average demand with the maximum lead time. Stock 
on hand is kept too high for items needing lesser investment. 

15. The values of other parameters ore determined in a manner out- 
lined below an entered :n Table 9 £ 

(a) The maximum vak's of stock numbers, work stations, J40 

and the fin ish ed products are known from Table 1; Figures 
13 a and 14. The values of producer 1 orders queue 
length before a work station depends on tht particular 
production system. 
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TABEE 8 

iythttohy iyicey^ioit 
Parameter 5TY3 (i, j), IC = 7 


jT 

JL 

1 ■ 

nr 

Li 

I " 
l 3 

■1 

■i 

mm 

H ■ 

Q 

o 

T 

\ 

s 

14 1 

19 


2 

10 

0.26 


30C 67.44 





200 

1 


1 

1 

0.012 

60.0 

1500 

6.0C 6.0 

5.4 

S.G 

2.0 





1 

1 

0.00016 

24.0 

52500 

0.08 0.03 

0.072 

8.0 

2.0 





1 

1 

0.00054 

32.0 

39000 

0.27 0.27 

0.2* 

8.0 

2.0 





1 

1 

0.008 

60.0 

9000 

4. CO 4.0C 



2.0 





1 

1 

0.0025 

60.0 

4500 

1.25 1.25 

1.125 

8.0 

2.0 




- 

1 

1 

0.0005 

31.0 

4500 

0.25 0.25 

C.225 

8.0 

2.0 




* 

1 

1 

0.0035 

70.0 

9000 

1.75 1.75 

1.575 

8.0 

2.0 




9 

1 

1 

0.0035 

70.0 

3000 

1.75 1.75 

1.575 

8.0 

2.0 




10 

1 

1 

0.002 

58.0 

4500 

1.00 1.00 

0.9 

8.0 

2.0 




11 

1 

1 

0.006 

71.4 

1500 

3.00 3.00 

2.7 

8.C 

2.0 




12 

1 

1 

0.0012 

60.0 

1500 

0.60 0.60 

0.54 

8.0 

2.0 




13 

1 

1 

0.00004 

8.0 

10500 

0.02 0.C2 

0.018 

8.0 

2.0 




14 

1 

1 

0.00002 

5.0 

4500 

0.01 0.01 

0.009 

8.0 

2.0 




15 

1 

1 

0.0033 

68.0 

1500 

1.65 1.65 

1.485 

8.0 

2.0 




16 

1 

1 

0.0014 

56.0 

1500 

0.70 0.70 

0.63 

8.C 

2.0 




17 

2 

10 

0.0209 



10.45 






1 

18 

2 

10 

0.0364 



18.20 






1 

19 

2 

10 

0.092 



46.84 






1 

20 

2 

10 

0.112 



56.94 






1 


('Table continued) 











TABLE 8 (Contd. ) 


EZBIMAEIOIT OB COLUI'DTS ; 


Column 


2 

3 

*4 

5 

6 

*10 

*11 

*12 

*13 

**14 

***-19 


Description 

Code indie sting whetner the item is purchased (l) or 
or manufactured (2). 

lot size. 

Carrying cost per unit per we ek in rupees. 

Reorder cost per order in rupees. 

Initial stock on hand in number of items. 

Unit cost in rupees. 

Regular price per unit in rupees. 

Discount price per unit in rupees. 

Average lead time in number of weeks. 

Standard deviation of lead time in number of weeks. 
Out of stock cost per unit per week in rupees. 

Batch quantity. 


* Applies only to purchased items. 

** Applies only to finished goods. 

*** Applies only to manufactured items. 



si 


b. 

c. 

d. 

e. 

f. 

g* 

h. 

i. 

i* 

k. 

l . 

m. 


n. 


0. 


P- 

q.* 
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PRODUCTION OP "IP ST/JIB" 

Costs, Parameters, and Constraints. 


Number of Stock Numbers (150) 

Number of Work Stations (J50) 

Maximum number of components permitted per SB (J40) 
Queue size at all work stations (lQ) 

Number ox finished goods (m) 

Number of allowed outstanding purchase orders (NOO) 
Number of sub-periods (iSP) 

Number of expected runs (iBOW) 

Starting time period (lip) 

Shift change oost per work station (sc) 

fixed overhead per week (CH) Rs. 

Variable overhead (VOH) 0.4 times total labour 

Initial number of shift minutes (NMS) 

Random number generator seed values in generating 
random demand values (lx) 76500267 and random lead 
time values (ip) 76 500267 . 

Maximum allowable inventory (PITVTY) Rs. 

Penalty factor of excess inventory (VTIIV?) 


20 

7 

5 

20 orders 

1 

200 orders 

6 

25 

1 

3,000 

cost. 

480 


1,50,000 

0.25 


Actual average unit production cost of final product 


(ACUNCO) Rs. 80 

Actual average weekly inventory costs (aCUTVM) Rs. 1,200 

Number of minutes in one time increment (TING) 1 

Factor of Discount Order Quantity (fdOQ) 1.2 
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(b) The number of allowable outstanding purchase orders (lIOO) 
depends on the purchase policy of the management, lead 
time of the items and the suppliers policy. If the numcer 
of purchased components is high, a high value of 1100 is 
chosen to avoid shortages. 

(c) Humber of sub-periods represents the number of to iking 
days per week. 

(d) Expected number of simulation periods is the planning hori- 
zon for which demands have been estimated. Care is taken 
to select this value such that capacity of the production 
system is not grossly exceeded during this time period* 

(e) Simulation periods are numbered similar to the demand 
values. Hence, the value of starting time period is speci- 
fied as one. 

(f) Shift change cost includes the cost of changing personnels 
and an arbitrary penalty for haphazard change. This value 
is not known for ’IE STAGE ' , since at present the industry 
under striy works for only one shift during each day. 

(g) High values of seeds arc chosen to generate random numbers* 
This prevents the random numbers from approaching the zero 
values. The random numbers of zero values create computa- 
tional difficulties in calculating normal .eviate values 
for demand aid lead time. 
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(ii) Time increment used in the simulator is the highest common 
factor of the various processing and set-up times* 

In addition to the above mentioned data values s others shorn 
in Table 9 are determined for the particular production system under 
study* These data values are necessary for the initialization of the 
simulator* Punched output of the supplementary computer programs - 
?T Emplosion Program T? and n Forecaster ProgreP. n ? are also used as input 
data for the simulator. After initialization^ MOh-PBDSIIl ¥ requires 
demands as input data during the simulation run* These values are 
specified in the second column of Table 15* 



GHAPTEH VI 


Essjncs aud Discussions 


The working of the alOD—PBDSILI V simulator has been tested, 
with the help of the data obtained from J.K. Electronics. The follow- 
ing parameters are accounted in the simulator in order to develop the 
model for the production system under study : 

1. The period of simulation is kept as are wed: comprising of six 

wo riding days. This is done due to the fact that a weekly planning 
period is observed in J.K. Electronics. Veekiy demands are gene- 
rated and production planning is done accordingly. 

2. The entire simulation run consists of 25 simulated periods. The 
planning horizon of 25 weeks is choosen so as to coniorm with the 
actual conditions prevalent for the ! IE STAIE' product at 

J.K. Electronics. 'IE 3TJGE' is a newly introduced product for the 
assembly of TV sets. Hence ? the past demands are available for 
only 25 weeks. This planning horizon is considered sufficient to 
analyse the system's behaviour as events for e"very minute are 
simul ated . 

The forecasted demand, of 'IE ST HE 1 is determined from 
the past production data. It is assumed that the production data 
provide good approximations to the demand aid sales figures. The 
forecasted demand is estimated by the forecaster model. As shown 
in Eigure 15, a sudden drop in ihe production of 'IP STJ&E' occurs 
from first week to the fifth week. This drop in production is 
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attributable to the catastrophic events of strike and labour troubles 
during these periods® These periods have also been included in the 
past demand data in order to test that suck events can be taken care 
of by the proposed simulator. 

6 • 1 To re caster 

The various decisions made In a production system are 
either directly or indirectly based on the forecasts, Therefc re 5 in- 
accuracies and errors in forecasts result in non-optima! decisions. 
Using the line ar— quadratic— cyclic forecaster^ the forecasted demands 
match closely with the actual production values for past data as is 
evident from Figure 15* Assuming that the past demand distribution 
has a bearing on the future demands 3 the proposed forecaster can be 
used for future plaining of the production. The development of fore- 
caster requires the estimation of the frequency of cycles as has 
already been discussed in Chapter III. The ” Correlation aid Spectral 
Analysis” is detailed in Apt end i:: A. As mentioned earlier, these 
techniques have been usel tor the dete:minatior of frequencies. 

6.1.1 Use of Correlation and Spectral Analysis 

The past demand data of t IF STASE* axe tested by calcu- 
lating the different values of correlation coefficient (r^) and 
smooth spectral density function (g^) Tor various values of mazimum 
lag number (m) - the value of m varied in the range c ? 3 to 22. This 
is done in order to select a proper value maximum lag number (m) by 
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following ’'vrindow closing procedure” ss explained in Appendix A. The 
values of r^. and G^_ are aaslysed for different values of n. The 
significant results are observed when n readies the values of 6, 9 and 
19. The values of correlation coefficient (r, ) and smooth spectral 

K 

density function (& ) are tabulated (Tables 10 to 12) for the signi- 
ficant values of m. 

Correlogram. of past demand data does not indicate any 
significant periodicity within 95?» confidence interval as is observed 
in figure 16 . But, non-randomness of "die time-series is indicated by 
the non-zero values of r^. The smooth power spectra for the maximum 
lag value of 6 indicates "the demand as a simple Maikov process, as 
shown in figure _17. Ihen maximum lag has a value of 9, a spike in 
the power spectra indicates the periodicity of 3. Again a periodicity 
of 3 is obtained by using maximum lag value of 19. However, the 
spikes in power spectra for periodicity of 3 are found outside the 
lower control limit of 95^ confidence level. The only significant 
spike observed in the power spectra vithin control limits represents 
the periodicity of 9 for fie maximum leg value of 19 (figure 17). 

Correia b. on analysis indicates insufficient data values 
to identify periodic components. But, the presence of an eutocorre- 
lated time series is observed. Spectral analysis identifies the 
presence of only one cycle of 9 periods with 95/i confidence level in 
the past don and distribution. To confirm the periodicity of 9 and 
to identify other significant cycles, the iterative procedure as 
discussed in Chester III is followed. 
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lisle 


10 


CXSEBEIiiiEIOII /ITD SPECHL'i ilUiYSIS 
Maximum Leg (m) = 6 


S.Ho . 

Lag 

k 

Correlation 

coefficient 

*k 

Smooth Spectral 
Densi+v function 

n 

Frequency 

f 

1 

0 

1.0000 

130799.60 

C.OCOO 

2 

1 

0.9748 

56707.63 

0.0833 

3 

2 

0.9204 

1082.54 

0. 1666 

4 

3 

0.8503 

407.71 

0.2500 

5 

4 

0.7834 

233.45 

0.3333 

6 

5 

0.7699 

200.94 

0.4166 

7 

6 

0.7764 

174.84 

0.5000 


OOEBEDOGrE/M : N = 25 Mem = p -0.0417 


Standard Deviation = o~~ = 0.1998 

For 95% Confidence Level - 


Upper Limit (UL) = 0.3500 
Lower Limit (ll) = - 0.4333 
POWER SPECTRl : V = 7.666 r ( 

Mean = G = 27086.69 

Por 95% Confidence Level - 

2 2 
Z 2.5%,V = 2 *° X 97 


* 




i 7. 1 


Upper Limit (UL) = 60400*0 

Lower Limit (ll) = 966 0*0 
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TEHEE 1 1 

coBHEircio:: ;:;d spsczfji xi/xysis 

Llaximus Lag (m) = 9 


S.Ho « 

lag 

k 

Correlation 

Coefficient 

KI 

Smooth incctral 
Dcncix;; function 

f . 

Frequency 

f 

i 

0 

1.0CC0 

194998 c 39 

0.0000 

2 

1 

O.S74S 

84444.60 

0.0556 

3 

2 

0.9204 

1783.22 

0.1111 

4 

3 

0.8503 

1011.73 

0.1667 

5 

4 

0.7834 

467.61 

0.2222 

6 

5 

0.7699 

258.76 

0.2778 

7 

6 

0.7764 

224.78 

0.3333 

8 

7 

0.7925 

229.10 

0.3889 

9 

8 

0.8324 

174. 5C 

0.4444 

10 

9 

0.8290 

171.82 

0.5000 

POWER 

SPECTRE 

V = 4 

.8889 ~ f - 

28376.45 


2 

*2. 

5 $S,V = 0,79 

X '//~5;t. 12 * 5 



UL = 72500,0 

II = 4-580.0 
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SC 


OOEBEL.JTIQH J11D SPSC1RJX ILiLYSIS 
ICaximum Lag (m) = 19 


S.Ho. 
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k 

Correlation 

Coefficient 

■v* 

~k 

Smooth Spectral 
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G k 

Frequency 

f 

1 

0 

1.0000 

405489.84 

0.0000 

2 

1 

0.9748 
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0.0263 

3 

2 

0.9204 

2991.94 

0.0526 

4 

3 

0.8503 

1776.90 

0.0789 

5 

4 

0.7834 

1811.03 

0.1053 

6 

5 

0.7699 

1601,07 

0.1316 

7 

6 

0.7764 

1199.98 

0.1579 

8 

7 

0.7925 

576.02 

0.1842 

9 

8 

0.8324 

520.68 

0.2105 

10 

9 

0.8290 

321.72 

0.2368 

11 

10 

0.8112 

304.76 

0.2632 

12 

11 

0.7828 

157.14 

0.2895 

13 

12 

0.7455 

219.65 

0.3158 

14 

13 

0.6839 

223.46 

0.3421 

15 

14 

0.5915 

284.36 

0.3684 

16 

15 

0.4390 

238.89 

0.3947 

17 

16 

0.2016 

175.26 

0.4211 

18 

17 

-0.2220 

150.47 

0.4474 

19 

18 

-0.5559 

156.96 

0.4737 

20 

19 

-0.8605 

203.01 

0.5000 

POWER 

8 

8 

1-3 

£ 

#• 
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= 29794.59 



2 


2 



ON 

! C\J 

K 

V “ °*° 5 

Z ^7.55&,V~ 7,1 



UL = 108200.0 


LL 


758.0 





0 ^ n ^ „ 

JL. ^ W' U 


— S*“*j*> 



FIG. 17. SMOOTH POWER SPECTRA OF PAST DEMAND DATA 

(95% confidence level } csates) 
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6.1.2 Use of Iterative Procedure 

By using the iterative procedure the periodicities of the 
three cycles in time-series are found to be 9, 11 snd 14- periods 
(Table 13} • The periodicity of 9 is also confirmed by the spectral 
analysis. The most significant cycle has tine periodicity of 11 ana 
an amplitude coefficient of 11.09. The iee-t significant cycle has 
periodicity of 14. degression analysis also generates tne estimates 
of other demand parameters for minimum absolute error in forecasting. 

The demand coefficients axe tabulated in Table 13. 

The total absolute error in forecasting for the 25 past 
weeks is 65. The total actual demand for those periods is 3505. Thus 
the forecaster gives an error value of 1 .85% over 25 periods. The 
value of standard deviation of demand is 3.1° for average demaad 
of 140 uni ts of the finished product. In percentage this gives a 
value of 2.33. As shown in Figure 15, the forecasted values tally 
pretty well with the actual demands. This comparison has been carried 
out for 99.73$ confidence level in forecasted demands. The forecasted 
and actual demand values for ’IP STJC-E’ are tabulated in Table 1 j 
and the absolute difference per period for them is plotted in Figure 18. 

6 . 2 Utility of the Inventory Policy 

The inventory policy followed currently for the produc- 
tion of ’IF STiGE' is to keep initial stocks to leet the demand 
during the entire planning horizon i.e. 25 weeks. Tie purchase or- 
ders for raw materials are placed at regular time intervals of 16 weeks 



TABLE 


E0230ASTE2 PLRJ3EEP.S 

1=1 Eor ”11? Stage” 


IG 

Variable 

Value 

11 

A (l) 

97.644529 

12 

3 (I) 

1.56 3497 

13 

0 (l) 

0.1C46SS 

14 

D (i) 

7.833^83 

15 

E (i) 

4.5 

16 

S' (i) 

11.085214 

17 

G (I) 

5.5 

18 

H (I) 

6.752066 

19 

u (I) 

7.0 

20 

IK (I) 

3.162278 




to meet the forecasted demands. This results in large investments in 
initial stocks. The inventory carrying costs are also increased sig- 
nificantly. 

The proposed inventory policy generates the pertinent in- 
ventory parameters for the inventory model, The parts requirements 
of finished products are calculated by employing n Explosion Model* 1 * 

The estimated values of the inventory parameters for each purchased 
item are listed in Table 14. Initial stocks are maintained for the 
Increase in demand considering the variable lead time for each pur- 
chased item. Sufficient raw materials are maintained at the time 
period zero to meet the average demand during the maximum possible 
lead time of each purchased part, for *IE STAGE* this inventory po- 
licy results in reduction of initial stocks from about Rs. 7 S 00 ? C00 to 
Rs. 1 9 21 9 342, This represents the seven fold improvement over the 
inventory policy currently used for *IP STJGE 1 . 

The simulation results in significant improvement of the 
total inventory costs over the actual inventory costs (average). This 
improvement in inventory cost increases from 16,84$ in first week to 
34.17$ in the 25th week of simulation as shown in Eigure 19. These 
imp rovements in inventory costs are obtained mainly due to the reduc- 
tions in inventory carrying costs, 

6.3 Effectiveness of M0D-PB3SIM V 

The main performance measure of MGE-PBOSIT Y is the !! Unit 
Production Cost (UPC)" of the finished products. In calculating Hae 



TABLE 1^, 
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HTVEBTOBT PARAMETERS FOR PURCHASED ITEMS 


a. Longest Process Time = 18 Minutes 

b. Estimated Plant Capacity = 160 pieces per 

c. Factor of Discount Order Quantity = 1.2 

week 



Renumbered 

Actual 

Unit 

Maximum 

EOQ 

DOQ 

Stock Number 

Stock 

Requirement 

Weekly 


IT A3 


Number 

of SIT 1 (IR) 

Usage 


(i 9 9/ 


(I) 


Rate (R) 



17 

2 

1 

160 

1267 

1518 

18 

3 

35 

5600 

40988 

49185 

19 

4 

26 

4160 

22204 

26645 

20 

5 

6 

960 

3795 

4554 

6 

6 

3 

480 

4800 

5760 

7 

7 

3 

480 

7715 

9258 

8 

8 

6 

960 

6197 

7436 

9 

9 

2 

320 

3578 

4293 

10 

10 

3 

480 

5276 

6332 

11 

11 

1 

160 

1951 

2342 

12 

12 

1 

160 

4000 

4800 

13 

13 

r*7 

l 

1120 

21166 

25399 

14 

14 

3 

480 

15492 

18590 

15 

15 

1 

160 

2568 

3081 

16 

!6 

1 

160 

3578 

4293 

5* 

17 

1 




4* 

18 

1 





3* 

19 

2* 

20 


i 

1 


^Manufactured Item 


IT 
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values of UFO, all relevant oasts (Eefer O^ter II 
included, except the oveiheais. Jlllocation of overt 
particular product out of the several products beir* 
the npst difficult task faced iy the industries. 


for deupils) are 


oads to a.17 cr±e 
msr^facturei ? is 


The value of UPC/week is calculated iron ie cumulative 
average of the costs over past periods. Figure 20 represents the im- 
provement in UPC over the actual value (averse). Phis improvement 
is due to the reduction in inventory costs. In the first week of 
simulation the improvement in UPC is only 6.142 as compared to the 
13.86% achieved in the second week. This is due to the fact that; 
cost of initial stocks and other overhead costs are contributed to 
the value of UPC in the first period of simulation. 


The UPC of ’IF STijCE' increases from second period of 
simulation onwards to the last simulated period (Table 15). This 
also causes the drop in the percentage improvement as shown in Figure 
20. This is due to the fact that the costs of purchased items are 
also included in calculating the value of UPC for each period. The 
parts for viiich purchase orders are placed in the initial periods 
start arriving only during the later periods. ’ 


The effectiveness of MOU-PBOSIU Y is also .-judged by 
analysing the values of percentage utilization of production capar- 
city and the total plant operating costs during each of the simulated 
periods. The values of these vaiiables-hnve been pl< tted 
with respect to time in Figure 21 and_22. The charac oerstics of these 





21 ABIE 15 

RESULTS OF SIMUIu'JTIOB HJ3 

(i) Actual Averse Unit Cost = Ks. 80.00 

(ii) Actual Average Weekly Inventory Cost = Hs. 1,200. CO 

(iii) Calculated Weekly Out of Stock Cost = ?_s. 0.00 

(iv) Total Weekly Penalty = Rs. 0.00 


1C1 


io. iActuallFore- jlversge Inventory ][Total potal {Total {Weekly faumulativ 

of {Weekly Jcasted {Unit $ Carrying {inventory *?roductionflnvcntory|piart fpiaat 
Week {Demand {Weekly iCo st ][ Cost {Cost {Tice Rvalue {Cost {Cost 

{ {Pemand{Bupees { Rupees {Rupees ^'finutes {purees I Rupees {Rupees 

1 2 3 4 5 6 7 8 9 10 


1 

123 

121 

75.09 

264.48 

997.88 

2 

111 

111 

68.91 

250.23 

983.63 

3 

101 

99 

69.01 

237.08 

970.48 

4 

90 

90 

69.08 

223.97 

957.37 

5 

81 

86 

69.14 

212.25 

945.65 

5 

94 

90 

69.18 

199.48 

932.88 

7 

107 

101 

69.21 

185.60 

919.00 

1, 3 

115 

115 

69.22 

169.67 

903.07 

9 

134 

129 

70.00 

164.63 

898.03 

iO 

135 

138 

70.86 

159.65 

893.05 

\ 11 

136 

142 

70.90 

152.83 

886.23 


137 

141 

71.01 

145.92 

879.32 

\ ;3 

138 

139 

71.12 

139.00 

872.40 

; .4 

140 

138 

71.20 

131.32 

864.72 

< 15 

■ 16 

144 

140 

71.27 

124.55 

857.95 

147 

145 

71.38 

120.11 

853.51 

17 

150 

152 

71.49 

115.79 

849.19 

: 18 

156 

159 

71.59 

111.96 

845.36 

19 

162 

166 

71.68 

107.33 

840.73 

i 20 

168 

171 

71.78 

102.00 

835.40 

21 

175 

176 

71.87 

95.59 

828.99 

t 22 

181 

181 

71.95 

87.73 

821.15 

t 23 

187 

187 

72.02 

78.62 

812.02 

c 24 

193 

195 

72.09 

68.25 

801.65 

1 . 25 

200 

204 

72.16 

56.55 

789.95 


I 


11080 

114769.09 

12578.73 

12578.73 

10150 

106183.99 

11887.42 

24466*15 

8280 

99522.47 

10526 .80 

34992.95 

8280 

93472.71 

10513.69 

45506.64 

7360 

87681.68 

9827.56 

55334.20 

8280 

81623.74 

10489.20 

65823.40 

9200 

74817.73 

11149.71 

76973.12 

10120 

6" 057. 28 

11808.19 

88781.31 

11020 

77737.14 

12478.44 

101259.75 

11760 

72108.58 

15292.81 

116552.56 

12500 

67375.15 

15400.85 

131953.40 

12880 

63157.27 

15353.27 

147306.68 

12340 

59590.75 

14625.04 

161931.71 

11960 

55809.29 

-14474.69 

176406 .40 

12500 

52044.97 

15099.99 

191506.39 

12880 

49589.78 

15165.89 206672.28 

13240 

47 165 .88 

15530.15 222302.43 

13800 

45463.14 

15801.39 

238103.82 

13980 

43376.26 

15947.95 254051.77 

14720 

40936-38 

16558.59 270610.36 

14720 

38131.16 

16389.61 

286993.96 

14720 

34922.18 

16341.65 

303335.61 

1 r7 20 

31434.04 

16304.67 319640.28 

14720 

27745.40 

16269.17 

335909.45 

14720 

23928.00 

16232.80 

352142.25 
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closely resemble those of the demand distribution (Figure 15). This 
indicates that the demand distribution may be treated as reo resent sr* 
tive of the system performance. 

In calculating the values of plant opening costs for 
each period the approximate values o: o” -he'--; c-rts hove also been 
included. The plant operating costs cculd mi be compared with the 
actual values. This is due to the dilfirulties faced in estimating 
actual plant costs, since the allocation of overheads to a particular 
product could not be dete min ed. 

MOB-PBDSIM 7 simulates the production system for each 
minute. For 'IF STJS&E' the simulation is carried out over 25 weeks* 
This is quite a considerable period for the system to stabilise. The 
external perturbations settle down sufficiently over the .simulation 
period. The initial disturbances in the demand of 'IF SliGE' due to 
str ik es etc. (Figure 15) causes disturbances n other variables also 
(Figures 20 to 22). Bit these fluctuations in the variables settle 
down considerably by the end of the sinu-atm. 

In addition to the effectiveness measures mentioned above, 
M0D-PBQSD.1 V also provides the user with detailed reports. The printed 
reports for each period of simulation provides the production manager 
with the necessary information regarding the production s; item e.g. , 
production status, inventory position, ouista-rdin^ pu’iiase orders 
etc. This facilitates him in making both short -te-m ’’d long-term 


decisions 



CHAPTER VII 


CONCLUSIONS AND RECOIIX'DillCITS EOS 

IUKDHER UOBK 


7.1 Conclusions 

MOD-KDS3H V incorporates various inportant Unctions cf 
the production systems. The Amotions considered ar- forecasting, 
inventory control and production planning. Tne simulator generates 
a set of decision policies. These decision rules provide the man^e- 
ment with certain guidelines for rational planning. In developing 
the model an attempt is made to achieve an "adaptive" control systan- 
one that monitors its own performance and at leasx partially updates 
its own parameters. Every care has been taken to incorporate the 
realistic features of a production system in KOD-PBDSBI 7. 

The following inferences are drawn from this study : 

1. The simulator uses cost minimization as the optimization criteria. 
The various sub-objectives are satisfying customer demands, opti- 
mum plant utilisation etc. 

2. The model structured for MOD-PEOSUI V is very general - it 
permits application to varied manufacturing situations. The 
model is validated for the sub— system of a typical Indian indus- 
try manufacturing television sets. 

3. The formulation and validation of the model for large complex 
production system is prohibitive due to the res' -action imposed 
by the Computer for vhich the simulator has bee .eveloped. 



Effectiveness of the simulator is best judged by carrying out a 
comparative analysis. The results obtained from the simulation 
run are compared with -the actual performance of the production 
system. 


Demand forecasting is an important feature of liQD-PHOSEI V. A 
line ar-quadratic-cyclic forecaster is used in the proposed model 
with capabilities of handling upto three cycles in the time - 
series. Dor the 25 weeks of past demands for the TV industry 
under study only 1.85$ of forecasting error is found. Thus the 
proposed forecaster is sufficiently accurate and flexible for 
application to different types of demand distributions. 


. The inventory model developed in the simulator provides signifi 
cant improvement over the inventory costs of the production system 
under study. By carrying out the simulation the investment in 
initial stocks are reduced by as much as seven tines. The reduc- 
tion in inventory costs varies from 16.84$ to 34.17/'-’ Over the 


simulation run of 25 weeks. 


7. 


8 . 


The simulator provides an overall view of the operation of prod 
tion system «iti respect to time, me production mareger gets 
useful information for each simulated period e.g., production co- 
sts of final products, plant utilisation, total plant operating 
costs, production and inventory status etc. 


The simulation is carried over a planning honzo 
ability of me control system to respond ad; 


to test the 
st to external 



perturbations. The model as tested for the production system 
under study is found quite effective to take care of catastrophic 
events like strikes and labour troubles. The results of simula- 
tion indicate the presence of a responsive control system vfcich 
results in significant improvement over "the actual system per- 
formance. The system is found to be sufficiently stable over 
the planning horizon of 25 weeks, since the simulation period is 
a minute. 


7 » 2 Ec co m mendations for Further fork 


There are several channels for extension in the preset- c 
work. In attempt may be made to relax some of ihe assumptions -edc in 
developing the simulator. More realistic simulation would achie/ed 
by incorporating the following features of production systems on 
MOD-PBDSIM V : 


1 . The forecaster developed in the present model is Capable o so 
ling linear, quadratic and cyclic demands. Spending on me 
demand distillation of a particular product, me .r ^' cr esa 
be rrodificd to account for different trends (polynomial, logarith- 
mic, exponential etc). *cn simulation is carried oux over the 
entire life cycle of a product, then growth curves lihe Pearl 
and Oomperta (24) can be used to advantage for forecasting. Phis 
is generally done for products with short Ufa e‘' ole - 
2 . In the present „* demand is ronsidered only as a time dependent 
variable, since it's relationship with other vast abies is unhno*.. 


ICc 


In real life situations tie demand of a product may have a 
correlation with several variables - seasonal sales, demands 
of other products etc. The forecaster car. be codified to incor- 
porate the multiple correlation of demand, if the inter-- relation- 
ship and inter-dependence of variables is precisely known. 

3. Processing times at work stations are not essentially determinis- 
tic. a random variable processing tine car. be considered for 
each of the operations. This requires the information regarding 
processing time distribution. 

4. The quality (percent defective) of the incoming raw materials and 
of manufactured items in the production system affects the values 
of other variables. The quality can be considered as a random 
variable. This requires complete information and data connected 
with quality of all items o.g., parameters of quality distribu- 
tion, quality control costs etc. The relationship of quality of 
items with other variables can be incorporated in the model. 

5» In real life production systems the machines experience breakdowns. 
The breakdowns occur at random intervals aid the repair time is a 
random variable. These features can be incorporated in the simu- 
lator if various parameters and distributions for plant mainte- 
nance are precisely known. 

6. The proficiency of all machine operators is not ?aae, as is assumed 
in the simulator. To improve on it, requires a thorough evaluar- 
tion of each of the operators. After analysing ne various human 
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and. other factors, planning can be done for the operator’s train- 
ing, hiring, firing or replacement for cost minimisation criteria 
( 7 ). 

7. In the proposed model the assembly lint is assumed to be balanced. 
An assembly line balancing algorithm can be accommodated in the 
simulator for any particular manufacturing system. 

Each of the above mentioned features require an indepen- 
dent research in the related field. Incorporation of these features 
in the simulator, necessitates a thorough comprehension of the inhe- 
rent complexities and inter-relationships prevalent in ihe system. 
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appendix a 


COEKELJffilON MB SPECTRAL M ALYSIS 

The parameters of any time series model cm be estimated 
to greater accuracy if sufficient past data of ihe system are known. 
The stationary of the model is also identified by using the data. 
Correlation and spectral aialysis are two of the important tools in 
the identification oi the periodic and persistence components of the 
time series* 

Persistence refers to the linkage between successive members 

of a time series* The sequential dependence m a time series can be 

tested, by autocorrelation analysis* The autocorrelation function of 

a time series describes the general dependence of the values of the 

data at one time on the values at another time* The sample value of 
tli 

the h * serial correlation coefficient (r^) between and x^ t ^ can 
be calculated by - 

1 Jfch 1 N-h IT-fc 

N-h x i X i+h ” 7 v2 if X i 21 X i+h 

_ __ irJ QmQ. i=1 


r. = ■ 
h 

s. • S - ,v 
i l+h 


h = 0 S 1 ? 2, 

where 


N 

= total number of data values. 

X t 

= given time series at time t, 

h 

= lag number s 

m 

= maximum lag number 5 
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s ± is the standard deviation T/shich can be calculated by 



and s ±+h is .liven #ien x_.^ is cubs tit?;. -jed fcr x^ in 'die above 
equation. 

The correlation coefficient (r , ) "hida are functions of h 

can be represented graphically. Tire g^ayhicaL representation is 

called ‘Oorrelogrcj’. The value of r. as given by equation (A-l) does 

n 

not depend o i t since the process is assumed to be stationary. Though 

r, is defined only for discrete values of h. The plotted points are 
h 

generally joined by straight lines. The important properties of the 
autocorrelation function are summarized beiov : 


1 . 

2 . 

3. 

4. 


r Q - 1, 




S' 1 


for h = 0 
for all h. 


r h " r -h 


i,e. s is symmetrical about h — C. 

Tor a pure random process = 0, for h £ 0 S but _or a 
finite sample in fie rsodom process, fie computed values of 
, may be - fferent from zero Out to sampling errors. Sene- 
rrtLly for smell ample sites, to statistical testing can be 
done for the values of h > 1/10. A series •.* be considered 
random if r h is not Sign! f icurtl;- different from zero, fcr a 
circular random series, r h is nearly noemall distributed uith 
a mean — 1 / (N — 1 ) and a variance of (5 - - 1 ) / ( !r ~ 1( 
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So the confidence limits for probability level (i - o( ) are 
given by 


- 1 
B-1 


" t (</2) 


/( I - 


.H - 1 


< ^< 



^(^ 2 ) 



U-3) 

vbere is ibe standardised noxmal variate for confidence 

level (l - ). If the computed value of r^ lies with in this 

interval? is not significantly different fiom zero and the 
process can be considered as random vith (l - -\ ) confidence 
level* 


5# If cyclicity is present? then the cor relog ram will also be pe- 
riodic with the same frequency* 

6# For a simple Markov process? the correlcgram is a monotonously 

decreasing curve ? for r^ 0. lor negative r^ a strong high 
frequency oscillation with period unity is indicated and the 
cor re log ram has a decreasing but nonvanishing amplitude* 

7* For a simple model with equal weights, the corns log ran decrea- 
ses with h and vanishes for m >. xl * 

Autocorrelation function (rg establishes the influence of 
values at present time over the values at fliture time in ihe given 
data. A sine wave, or any other deterministic data, mil have an 
autocorrelation function vhich persists over all time displacements, 
as opposed to random data which diminishes xo zero for large time 
displacements. An autocorrelation measurement dee-iy P™^ es a 
tool for detecting deterministic data vfcidi ^ t: Qa ^ ed 311 a 
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random background. The details about a time series am generally 
better interpreted from the power spectral density function than from 
autocorrelation function. 


The power spectral density function of random data describes 
the general frequency composition of the data m terms of the spectral 
density of its mean square value. The power spectrum density function 
is usually estimated as the Fourier transformation of the autocorrela- 
tion function or Hie auto covariance function 3^ (k^" sato covariance) 
given by - 

\ = i^T 1- (l i - K+k - * i+k > + *i • s i + k 

i=1 

k = 0 , 1- 2, s m (ii-4) 


where 


x 


= sample mean value of 3U 
k = lag number 


if the transformed record is stationary with x - 0 ? then transxcimdd 


auto covariance function is given by 

I-k 


N 


S x. ♦ x. _ 

— * i l+k 


then - 1 fC 3^ ^ 1 


(A-5) 


i-1 


For sampled data from a transformed record ^ vihida is 
stationary with x = 0 ? a raw estimate of a true power spectral 

density function is defined for an arbitrary frequency f in 
range 0 f f “ 

m-1 ^ - Kk_f ' 2 /Tim. f 

f 

k=1 c 


G (f) 
x 


_2 

m 


, a Tik f n , £ Cos 

1+2 21 \ C0S ( ~T~ } f c 

o C c 


( A-6 ) 



values, 


Hence it is necessary to smoother the raw spectra. 
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here are several procedures available to obtain final 
practical smooth estimates Hiich have more suitable statistical - - 
ri ability properties. The frequency smoothing is done with a pro- 
cedure called "Hanning", for #iieh the smooth spectra aro given by 
( 2 ) - 

G = 0.54 G + 0.46 T 

o oi 

G m = °* 46 Vi + °- 43 \ 

G h = 0,23 G h-1 + 0,54 G h + °* 23 Vi’ h = ■*’ 2 ’ » 

U-9) 

A typical plot of power spectral density G_^ (f) versus 
frequency f is called power spectra. In order to estimate a sui- 
table value of maximum lag number m a "window closing procedure" is 
followed (3). The estimate is said to have high fidelity when for 
all frequencies, the bias (difference between the estimate end the 
true value) is small. For high fidelity, the bandwidth must be of 
the same order as the width of the narrowest important detail of 
the spectrum. This requires a small bandwidth sid a large el. 

The estimate is said to have high stability if the variance 
of the estimator is small. High stability requires a large value 
for the degrees of freedom and hence a small value of m. Generaii;. 
a compromise between these two requirements may be reeded. A-- m 
decreases, Hie raw spectra is appro ached and spurious peaks increases 



where 
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f c " CUt “° ff or ffyquist frequency = ; 1 

J 2 At 

m = maximum value of leg number k. 

St = tlme 3nteiYal over time horizon I between r samples 
of discrete data. 

A 

R q - value of transformed autocovarisnce function (% ) 

£. 

for zero lag number. 

A A 

R m ~ value of R^ for maximum lag number. 

Since the values of G^ (f) are calculated only at the m + 1 


special discrete frequencies where - 


h = 0, 1, 2, , m 


Ur-7) 


At these discrete frequency points, tie value of power spectral density 


function is given by 


G (f) 
x 


2 A a 

— E + 2 p R 
mo k 


p R Gls ( — ) + (-if 2 
t~L k m 1 1 r 


(A-8) 

Tile index h is the harmonic number. This will provide m/2 independent 
spectral estimates of G (f) which may be considered to be an average 

x 2 i 

of the spectra ove '• a bandwidth 3 = around frequency f. 

c e m 

Spectral estimates at points which are less than apart are 
correlated. (1 ) 


has a negative exponential distribution and so the 
sample value of G^ may be widely different from the true population 

*The references have been given at the end e . the j., oendix. 
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By comparison of the remits wilt theoKticsl models ar.d erpeete 
behaviour of the physics process, a suitable value of E may be 
chosen. 


The general characterstics of the power spectra for severs! 
generating processes are summarized below - 

1. A constant value time— series has no power spectrum. 

2. A sloping line has non— zero value of power density function at 
zero frequency. 

3. A simple harmonic tame series has non zero values at appropriate 
frequencies in the power spectra. 


4. For a pure random process, the spectral density, in the frequency 
range 0 < f ^ f c , is constant aid is equal to the average of 

tile (m + l) computed values of the raw spectrum. The sample 
spectral estimates are distribute d about the population spectrum 
according to Chi-square Distribution. The number of degrees of 
freedom for spectral calculation y= (2N/m) - (2/3). From tables 
of (x*V -'jj ) distribution, the confidence limits for a given level 
of confidence (l - ) can be obtained. The series is pure 



m + 1 


average of spect.um density 
function for pur • random process. 
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5. For a simple Markov process the spectral density decreases 
continuously with the frequency for f 0. 

6. For the general weighted moving average process the spectra i*= 

a sinusoidal distribution about the variance of the time series. 


The principal application for a power spectral density 
j^inctio ^measurement of physical data is, to establish the frequency 
composition of the data which, in turn, bears important relationsnips 
to the basic characterstics of the physical system involved. In the 
identification of the components of 'the time seiies, correlograr. Sid 
spectral analyses are to be used complementary to each other. ¥nxle 
correlogram may indicate the general nature of the components, adja- 
cent values are not independent; vhile in spectra, they are nearly 
independent. Moreover details may be more clear in spectra tna n in 


corxelogram. 

a, presence of periodic cycles m the tine series will be 
indicated by the periodicity in the correlogras Sid spikes et appro- 
priate frequencies in the spectral density function. Once the 
frequencies are identified other nefficients of idne series can be 

estimated by using regression ®alysis. *rom the past data *» 

n +• so -Pnnr+ion sad smoothed power spectral den- 
values of autocorrelation function Ji 

sity function axe calculated, for different vines of maximm 1« 
number m. Kter plotting hie power spectra for 11 *.» values 

v, i of m c® be decided vhi* Si'’, s sufficient 
of m, a suitable value of m a 

■ 0 _ mhe values of n may vaxy f* 0111 

details to analyze the time se 



Ui. 

N/ 20 to a maximum value of H/2. The nln+ , 

' ine pio t of correlogram remains 

same for all values of m and it ohnnu v „ „ ^ „ 

aua xx should be preferably plotted for 

maximum value of m. 

A computer program is available for autocorrelation, cross- 
correlation and spectral analysis ( 4 ). The program has been written 
in FOBTEAtT IV and the listing is included in .Appendix G. The computer 
program input requires number of samples IT, number of variables, past 
history data of time series aid different values of a. The output 
results contain values of autocorrelation function, cross correlation 
function, raw and smooth spectra, frequency, upper aid lower confidence 
limits for correlogram for 95$ confidence level, degrees of freedom 
and mean value of smooth spectral density function T? (f), for each 
value of m. Sufficient number of COMMENT Cards have been sprinkled 
to aid the user in understanding the general flow of control within 
the program. 
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APPEtDIX B 


EXPLOoIOlT PBDCEHJRE EOR PARTS DETERMINATION 

Explosion procedure is aa analytical tool for the conver- 
sion of sales forecasts to the estimate of production quantities. 

It is a systematic procedure based upon matrix algebra. It provides 
a good framework for performing the planning function by parts de- 
termination. This technique is specially usefhl vfcen large numbers 
of products and parts are involved. In this procedure a Gozinto 
analysis is carried out for a final assembly by drawing explosion or 
Gozinto chart showing the relationships of components, materials and 

t \* 

sub-assemblies (1). This directed network indicates the order in 
which parts go into other parts to make the final product. Eurther 
it indicates the number of parts of each tjpe needed to make one of 
the other type. Erom this the "lext Assembly Quantity Matrix (ll) n 
is developed. The I matrix shows the bill of material relation^iips 
for all parts, finished or unfinished. In I matrix a row of zeros 
indicates a final assembly, vhile a column of zeros indicates a 
raw material or a detail part. 

The goal in the explosion procedure is to calculate the 
"Total Requirement factor" or "Matrix I", showing the total quantity 
of each part necessary for all sub-assemblies aid assemblies. The 
diagonal elements of T matrix have values equal to un. ty. The 
relationship between T and H matrices is represented 
* The' references have been given at the end cu Ap f 
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T = f I - N] 

viie re 


(B-i) 


I - identity matrix of order equal to the total number 

of items (raw materials, arb-assemblies and finished 
products. 1 ) 


An algorithm is available to compute the inverse of (i - II) matrix 

( 2 ). 

Even after using the (i - S’) ^ algorithm the number of cal- 
culations required in G-ozinto analysis are quite large. In the 
"Explosion rind Netting" procedure (3) for parts determination, the 
parts are renumbered so that large stock numbers always go into the 
assembly of small stock numbers. By doing this the IT matrix becomes 
lower triangular matrix, for which inverse is easy to calculate. The 
calculations required in explosion and netting procedure are simplified 
by employing "Curado and Hassell' s" algorithm (4). The latest effort 
in this field is the "Plug aid Chug" technique (5), vihich further 
simplifies the Curado and Hassell's algorithm. 

In Plug and Chug method the computations are done backwards. 
Prom relationship (B-i) it is clear that since IT is a lower triaigular 
matrix, T is also a lower triangular matrix with unit diagonal ele- 
ments. The total requirement of the ith assembly or sub-assembly is 
given by the recurrence relationship * 


t 


k=i+1 


(B-2) 
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T = [ I - N] (B-l) 

vhe re 

I = identity matrix of order equal to the total numbei 

of items (raw materials, sub-assemblies and finished 
products.) 


An algorithm is available to compute the inverse of (i - N) matrix 

( 2 ). 


— 1 

Even after using the (i - F) algorithm the number of cal- 
culations required in Gozinto analysis are quite large. In the 
’’Explosion and Netting” procedure ( 3 ) for parts determination, the 
parts arc renumbered so that large stock numbers always go into the 
assembly of small stock numbers. By doing this the N matrix becomes 
lower triangular matrix, for which inverse is easy to calculate. The 
calculations required in explosion aid netting procedure are simplified 
by employing "Curado and Hassell' s” algorithm (4). The latest effort 
in this field is the "Plug aid Chug” technique (5), which further 
simplifies the Curado and Hassell's algorithm. 


In Plug find Chug method the computations are done backwards. 
Prom relationship (B~l) it is clear that since IT is a lower triangular 
matrix, T is also a lower triangular matrix with unit diagonal ele- 
ments. The total requirement of the ith assembly or sub assembly is 
given by the re cur re nee relationship s 


m 

q - Y. \± \ + h 

k=i+1 


(B-2) 
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where 

= ith column vector of matrix I. 

m — total number of parts (raw materials, sub-assemblies 
and finished products). 

~ ©lament in the kth row and ith column of S' matrix. 

If = column of the identity matrix with value of ith 
element as unity. 

The column vectors of matrix T are represented as — 


T = 




(B-3) 


The value of vector T is known, because I is a lower 
triangular matrix. The values of other vectors of T matrix are com- 
puted by moving backwards from vector m to the first vector and by 
using the recurrence relationship as indicated in equation (B-2). 

The computer program listing for Plug and Chug method has 
been included in .Appendix D. The program, written in BOETBAB IV is 
capable of handling upto 50 different items for each finished product. 
This restriction can be easily relaxed for bigger problems by modify- 
ing the Dr.Vk'iCJOU statement in the program. Before using Ihe program 
the item numbers ore renumbered so as to make the input data matrix 
B, a lower triangular matrix, as was explained earlier. The progran 
provides the results in the form of printed report and punched cards. 
This punched deck output represents the elemr-rts of total require- 
ments matrix T for the f inal product. These punched cards are useful 
as input data for other computer programs. 
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APPENDIX G 


COMPUTER PROGRAM LISTING 



I' "i 





' i ' . 1 ■ i I ; { ) * J > J ( > > I A - ( ■ ) » ( » ) 1 ■ 1 

" ' ‘ f 

< ’ u , ' ; ) 

i ■ i ( J i 1 ( J ) » J •” . > f' ) 

■ ( » ) 

I t/////// « »*1 10. JF P\S-, *i • F <- ' ! < F . 


) f I 


»?..///; 



*|Ml = NM+i 

ALLOWING LOOP CALCULATES HE OPTIMUM SID 

: 3R MINIUM TOTAL A3S0LUTE ERROR I M F0R = ! 

30 23 NM=2» M 

33 2 J MM= 2 , M 

IF ( NN. £3- MM ) 33 TO 20 

30 i 9 Ll=2»M 

IF C LL - = 3 - MM ) 33 TO 19 

IF (LL.E3.MM) 33 TO 19 


3F PERIODS 
ASTIM3. 


CALL F3RCS SUBROUTINE TO CALCJUTE TOTAL A8S3 JT 
FORECAST! MG. 


I M EACH CYCuE 


ERROR IM 


;all for: s ( n, id, jo, ierr, i eras, a 

i TOT, H,IMM» F, G,MM,lL, J,NM,MM i) 


»I3T3T,JDT3T,IERT3T,AA,8B 


D,E,I A3T 


STORE MIMIMUM VALJt: OF TOTAL ERROR Is| FORECAST! M3 AMO MO. OF 
PERIODS. 


IF (MINA3 T.EQ.O) 33 TO 10 
IF (MIMABT-IA3T0T) 23,1 ,13 
I) -IIMABT=I ABTOT 
MM3PT=NM 
MM3PT=MM 
_L0PT=LL 
19 3 0MTIMUE 
23 3 0MTIMUE 

; *** CALL FJRCS SJBR3JTINE TO 3ENERATE F0RECASTIM3 PARAMETERS FOR 
; MI IM IM JM TOTAL ABSOLUTE ERROR IM FORECASTI MG. 

3 ALL FORCS( N, ID, JO, IERR, IERA3, A, I3T3T , JDTOT, IERTOT, A A, BB, C , D, E, MI M 
DBT,H,NM3PT,F,S,MM3PT,LL0PT,U,MM, NM1) 

3RIMT 112, MM0PT,MM0PT,LL3PT 

112 =ORMAT { / / / / 8 X » 5 J ( 1H*)////8X»*M3. OF PERIODS PER CYCLE ( NM ) = *, 13 
1////8X,*M0. OF PERIODS PER CYCLE (MM) = *, 1 3/// / 3X» *N3. OF PERIODS 
2 PER CYCLE (LL) = *, I 3 ////8X, 80( LH*) //// ) 

3 RI MT 1J2, AA,BB,C, D,c, F,S,H, J 

132 =3RMAT ( / / 3X » A1 H**** THIS IS THE OE M AMD FORECASTER ****////5X, S3 
HD(T )=A + 3*T + C*T**2+D*COS ( PI *T/E ) +F*C3S( PI*T/3) +H*C3$ ( PI *T/ J ) +YAR// 
2//8X, *A * *,Ei 3.8, 5X»*B = *, 

3Ei9 . 8 » 5X , *C = *, EL8.9, 5X,*3 * *, El 8. 8/8X, * = = *, E13. 3, 5X,*F = *, Ei 
43.8,5X,*G = * » £1 8 • 8 , 5X »*H = *, E18. 8/8X, * J= *,E18.8//) 

5 RIMT 103,(1,10(1) , JD( I),IERR(I),IERAB(I),I = i,M) 

103 =3RMAT ( //3X, A3'H***THI S IS DEMAMO FOR PRESEMT * M' PERIODS *****/// 
1/3X,*PERI 30 ACTUA_ CALCULATED ERROR IN ABSOLUTE ERR3R*/15X»*DEMAM3 
2 DEMAMO FORECAST IN FORECAST*/// / ( 4X» 5 1 8/ /) ) 

>RIMT 1J4.IDT0T, JOTOT, IERTOT, MINA3T 
104. : 3RMAT ( / / BX, *T3TAL ACTUAL DEMAMD = * , I B////3X, *T0T AL FORECASTED 
DEMAND = *, 1 8////3X, *TOTAL ERROR IM FORECASTI M3 = *, I8////8X,*T3TA 
2. ABSOLJTE ERROR IM FO RECASTI MG = *,I8//) 

********* ************************************************************** 


THIS CALCULATES THE MINIMUM VALUE OF STAMDARD DEVI ATI DM . 

*********** ********** ************************************************* 


IERTOT=IERR( 1 ) 


b) 
\l 3 
121 


I 

I 

I 


£*TQT*IEiRT 3 T + IiRR(I) 
SJM 2=IERT0T**2 


^SI3 = FLDM( (I *1 • R2T-ISUM^)/{ ) 

; I SMA = S3R T ( PS 1 3 ) 

[F (I.Lc.5) G3 T3 S3 
;)SIG=FL3AT( I ER2T/( I-MM )) 

SlS v lAL = S3^r ( QSIG) 

SIS'1A = A W II MitSIG^Af S I 3 ^ A1 ) 

33MTINUE 

3 S.I MT I23»SIGMA 

= 0*MAT ( 3 Xt *AK = *»Ei8.B) 

JJ'ICH 12i»AA»B3tG»0»E» F» G» H» U» SIGMA 
=3RMAT ( 3X» FIS. S ) 


*##**#*** fr************************************************************* 


rHIS IS =3RECAST FOR MEXT M PERIODS 


* ********************************************************************** 
\)2 = N!*2 

33 50 I*M*M2 


53 jou )=AA+-B3*AI+C*AI**2+3*C0S( (22. *AI)/(7. *=))*. 5* F*C0S 1(22. *AI)/(7 
I.*G) ) +H*33S ( ( 22 . *AI ) / ( 7. #U) ) 

135 =3^AT 1 (//'3xlt3H^**rHIS~IS FORECAST FOR FJTURE • M‘ PERIODS**** 
I7///8X»* t> cRI3D F3RECAST*//( ^X t 2I8//) ) 

STOP 

:MD 


SJBROUTI N1E M «S«<t,ID.JO.«E»*.I«M.*-»”T.J3nr,IER7DT.»».»B.O. 

13,E* IABTOT» H, i4M» F* Gt^ |LL»U* UM» NM1) 


*****************************^^ 

. . , Or -A- <**■<*■ & 


*** THIS IS SUBROUTINE F3RCS *** jrtBS0UTrERR0^I^P 3;iEC&STIMS * 
**» PURPOSE - TO DETERMINE MINIMUM ABS3LJTE tKKJK 


rf UNru3C ~ itj * * ' * * '*** 

>I'ENS13* ID.U,),J D aOD),URMnDI.Ml J ,n).BIU,l J ),XE«W.L«) 



1 1 'I. '' ' 

5 = F j N / 2 • 

3 = FLUaT ("3 ':) 
j =F L OA T ( L L ) / :: . 
i PT :.i I = 

no j= ■» j 

l IUT')T = I <)T‘JT+Il){ J ) 
QinT = 1 D F3 T 
0^= '* 

3 A V ' > — f"> T jl /I '• I 
3D / L = . » i”': 

DO 7 0 = -i ''P'i- 

7 A(L, ! i) = . 


^ t D •_ F I N , U A 1 '! 1 1 T RY AK ADD H TO DFCRtASh HO, OF CALCULATIONS 

DAVG= A VCR A CL AuTUAL DEMaND IN N P£RI OPS 

iK= !)-[■ A VC 

T = J = P,.R1s)D 

D=T-N/2 

DO : J= . » <i 

3 = 1 b t J ) 

ak=d-cavg 

T = J 

D = T-lH/1. 

:g=cus( (/2.-.*T)/(v.*d ) 

COfc =Ci ! S ( ( , 1 ’ ,*r)/(7.*G) ) 
cfj:=CuS((o - .*rj/(7.*u} ) 

A )*',( , 1 )+AK 

A ( A ? ;> ) = A ( / » 3 ) +! I 
AC,4) = A(^.,4)+H*H 
A (i f . )»/.( „ )+H*AK 

=A< . , 4 ) + H**3 
A { 4 , i ) = A ( *t » 1, ) +• H # H * A K 
n K ) =A( V , A- ) +H**4 
A ( *r » b ) = 7, ( “■ , S )+H*ll*CO 
m:,5)=ai .,5 )+co 
A(2o)=A( j,5)+H*C3 
A ( 5 » C ) = A { i, 5)+C(J*CO 
A ( £ , fc ) = ; , ( » 6 ) + C J 2 

A ( „ » 6 ) = A ( , 6)+H*Ll72 

A ( 4 , t> ) = A ( 4 » 6 ) +H#H* COL 
A(5,6)=A( 5,6) +C<)*C02 
A(6,*.) = A( j j t ) +aK*Cl)2 
4(6,6)=A(6,6)+CL2*L02 
A ( 2 , 7 ) = A ( 1,7) +CP. a 
A ( 3 , 7 ) = A ( 3 , 7 ) +H *CO 3 
A ( 4 » I ) — a ( , 7 ) +H£H*C03 

A (5, 7 )=A( b, 7) +Ui*CUjs 
A (C»7 ) =4( 6» 7 ) +CG2#CQ3 
A { 7 » 7 ) = A ( 7,7) +CU3*CG3 
A ( 7 , _ ) = A ( 7,1 )+AK*C03 
2 A (5, a )=A( 5*1 ) +AK*C(J 
A ( ?. , 2)=DM 
A ( .3 , 3 ) = A ( 2,4) 


I =K 

3 I = 1+1 

A ( I , K ) =/■• ( K t I ) 

IF ( I —Nil 1 ) 3 , ' , 

* CO.iTli'L’ 

************************************£ 

DLT 1 - KM I. i Ai ;T OF i iRD '.: K \K ,Y PIVOTAL CONOcNSATIOM 

************* *****************: * 4 :*****:**:**^^^^^^ 


DP ;;3 Ii=i*Nf! 

DO 6 J = - » i PI 
DU b K = . » 'MM 

J , = J + .*. 

<*. = K 

I F ( K.Gl-. i I ) K .. = <* + .*. 

6 3 ( J » K ) =a ( J ,. , K„. ) 

<=*' 

L = « 

Di,T r R M = _ . 

M = 0 

5 DU .. I =K j MM 

11 IF UBS(o( L,L ) ) .LF.l i. Lti-25 ) GO TO 3 
Gil Td id 

£ M= M .+ „ 

DO 9 KM = L » !M 
LM=L+M 

3!KM,L )=B(KM,LM) 

B ( K M, LM ) - u ( KM » L ) 

9 DDT _RM = -i)(. r.-RM 
SO TO i, 

12 UTItJ = B{ 1 ,L )/M L,L ) 

DU . J=K,NM 

1J BU,J)=R( I,J)-o(L,J)*RATiO 
IF ( K-N M ) ..5,10,1:0 

15 L = K 
K =K+ * 

SO Tb 5 

2D DO c 5 L = i , ' JM 

25 DDT .RM = QEr_RM*H(L,L) 

SD TG ( 1,32, 13,34, 36 »37)»II 
3D B x = j c TERM 

31 C_ = i)ETEK(i 

32 D_ = 0_TE iM 

33 B_ = L)tTERM 

2^ f„ = d:temm 

26 D:. = ;j1TERM 
37 Fl = t)ETERM 
35 CUNTINU:- 

F=D,,/1 . 

U=(C ,.4Ui^F + C , .*F**2)/B1+DAVG 
B6=(-D1~2.*ua*F )/Bi 


METHOD 



D=-F x/B - 

F=0* /3i 
H=-F2/iU 
JuTiJT = C 

i;R'for=o 

I A&l 0T= 

DU 11 - I » < 

41 = 1 

jUf 1 ) = A;' + HR*4 1 + L*AI**;:+D^C0S< ( 22.*AI ) /{ 7 .*e ) )+.5+F*CQS( (22.*A1) / < 7 
1, #G ) ) + H u j ( ( .^• s > c AI)/(7. -£U ) ) 
i:rr( i )-jt'( I )-i D(i ) 
i;rab( i ) = i abs { i t RR(in 
1 ABT UT = 1 A 6 TU I + 1 ' RA ri { I ) 

J DT 0T= J 0 T "J T + J D ( 1 ) 

43 I ^RTuT = IE RT (JT + I LRRU) 

RirTURH 

END 


-it , ! * If.* it.. I J- ; ; 


r' L r j I ' 


' '.' ' ’• '£ ;• i} $ > ■!= .£ ^ ^ ,|t j|c .j; ;{•.;; s';# £ # £ £ ^ J: : 


*, L ; f C. . v,S F T.' !: 1 I 


j f & P. ‘j !j i J C i I ’ J “< U U A , m T I 


jAL'i ) l ' i 'J j 3‘k. PL 


JO ' , J J I J J 


;«> « ! p p ": ,'.•*< n s s - 

i 1 ( • : i V . i ! i ■ ** { £R I ‘ L t r L’-IC H A i •_ u P A K T uX MJL ■ 

t T .. 1 l/. 1 J ” i "l ;f. 1 LY TC "VC Cni (JJIT uF CCCC* 


!'■■) = ! J. J ) --- -.,(: t J ) 3 ■ J.u, C 


! ■ J ) = , i 


'■■ ■> 1 • L :.L - 

,Ul,< , M . Lif C |‘TS I BY P',;<T J FIj-\ b U 6 - b L c v r . L Y 


( i - •! 3 


. J- 


( { ’ ; J ) — IF J • ‘3 1 * 

A! : i ‘ 


-i= 3 


i'i i > ‘ i j)»:( Ct b! 
/ / 


.j » i j 


(,,.]}= j 


j J ) » J =. * 1 ) j i - _ * • ' 3 


f 


)•,*--!( "i , *i } 


4 


T ( 


i - f 
“ * ' s 
^;(Ji 

I ; 


} '>' i t f J ) 4 I l I f ' ^ ) 


] = J 

, r 

(/ /// 


i- . . i / 

r • • ' // ' 

f . , !/// 

r . , ! i ) 

i < i ) 

• !//// 


• - ! " ) 

( ■ ! - j J ) j J = _ j ! ) i i = _ » ') 

) 

t ) ( h * ) f •* ' • J i - L ' 1 Z . Of- ! J ’ '< T ’ = 
'■ ' 1 r <-T 1(0-17 ,.,)*///) 

,i j ' iJ :■ / / ) ! 

t •j ) 


if 


* , 1 _////:,.,.. 0 ( . i * ) » * 


j '*• ) 


1 1 l j 


nTp.i, 


OF P (i i< T £ li; iKIAr. J */// ) 


ooonoooooooooooooc-iooooor>ooc~>oo<nc-)oooooo<r>ooooo 


************ 


SIMULATOR PROGRAM 
♦ ****************** 


C *** THIS IS MAIN PROGRAM ******************************************** 
C ** PRODUCTION SYSTEM SIMULATOR - FOR A TV INDUSTRY 

** DESCRIPTION OF INPUT PARAMETERS - 
150 - NUMBER OF STOCK NUMBERS. 

J50 - NUMBER OF WORK STATIONS. 

J40 - MAXIMUM NUMBER OF COMPONENTS PERMITTED PER SN. 

LQ - QUEUE LENGTH (SAME FOR ALL WORK STATIONS}. 

NT - NUMBER OF FINISHED GOODS. 

NOO - NUMBER OF ALLOWED OUTSTANDING PURCHASE ORDERS. 

ISP - NUMBER OF SUBPERIODS. 

I ROW - NUMBER OF EXPECTED COMPUTER RUNS. 

IP - RANDOM NUMBER GENERATOR SEED VALUE FOR ESTIMATING LEAD 
TIME 

IX - RANDOM NUMBER GENERATOR SEED VALUE FOR ESTIMATING DEMAND. 
PST(ItJfK) - AN ARRAY CONTAINING PROCESS TIME, SET UP TIME, HOLD 
QUANTITY AND WORK STATION ORDER. 

I = STOCK NUMBER. 

J = WORK STATION NUMBER. 

K = A CODE THAT INDICATES TYPE OF INFORMATION. 

IF K = 1 - PROCESS TIME IN MINUTES. 

K = 2 - SET UP TIME IN MINUTES. 

K = 3 - HOLD QUANTITY IN LOTS. 

K = 4 - WORK STATION OROER (THIS NUMBER INDICATES THE 
POS IT I ON OF WORK STATION J IN THE NETWORK 
SEQUENCE OF STOCK NUMBER I). 

RMR ( I , J , K ) - AN ARRAY CONTAINING COMPONENT REQUIREMENTS DURING 
MANUFACTURING - STOCK NUMBER OF COMPONENT, QUANTIT 
OF COMPONENT, AND WORK STATION AT WHICH THE ; 

COMPONENT ENTERS THE MANUFACTURING PROCESS. 

I = STOCK NUMBER OF ITEM BEING MANUFACTURED. 

J = A NUMBER FROM 1 TO J40. THE VALUE OF J( 1 , 2 , 3 , ) 

INDICATES THE SEQUENCE IN WHICH A COMPONENT ENTERS TP 

process. ; 

K = A CODE THAT INDICATES THE TYPE OF INFORMATION. 

IF K = 1 - STOCK NUMBER OF COMPONENT ENTERING THE PROCESS 
K = 2 - QUANTITY OF COMPONENT (IN UNITS) REQUIRED TO , 
PRODUCE ONE UNIT OF STOCK NUMBER I. 

K = 3 - WORK STATION NUMBER AT WHICH THE COMPONENT 
ENTERS THE MANUFACTURING PROCESS. 

OHR ( I , J ) - AN ARRAY CONTAINING MONETARY RATES IN RUPEES PER 
MINUTE FOR MEN AND MACHINES. 

I = WORK STATION NUMBER. 

J = A CODE THAT INDICATES THE TYPE OF INFORMATION. ; 

IF J = 1 - MAN RATE ON FIRST SHIFT. 


• J = 2 - MAN RATE ON SECOND SHIFT. 

J = 3 - MAN RATE ON THIRD SHIFT. 

J = 4 - MAN RATE ON OVERTIME. 

J = 5 - MACHINE RATE WORKING. 

J = 6 - MACHINE RATE IDLE. 

VOH - VARIABLE OVERHEAD AS A FRACTION OF TOTAL LABOR COST. 

UH - FIXED OVERHEAD Ifl RUPEES. 

SC - SHIFT CHANGE COST. 

FTAB ( I , J ) - AN ARRAY CONTAINING ALL INVENTORY OATA FOR EACH 
STOCK NUMBER IN THE SYSTEM. 

I = STOCK NUMBER. 

J = A CODE THAT INDICATES THE TYPE OF INFORMATION 

IF J = 1 - A CODE = 1 - STOCK NUMBER I IS A PURCHASED ITEM 

2 - STOCK NUMBER I IS A MANUFACTURED 

ITEM. 

J = 2 - LOT SIZE IN UNITS. 

J = 3 - CARRYING COST PER UNIT PER WEEK. 

❖J = 4 - REORDER COST PER ORDER. 

J = 5 - STOCK ON HAND IN UNITS. 

J = 6 - COST PER UNIT. 

J = 7 - ISSUES THIS WEEK. 

J = H - RECEIPTS THIS WEEK. 

= 9 - DISCOUNT ORDER QUANTITY IN UNITS. 

*J = 1J - REGULAR PRICE PER UNIT. 

*j = 11 - DISCOUNT PRICE PER UNIT. 

*J = .2 - AVERAGE LEAD TIME IN WEEKS. 

s 13 - STANDARD DEVIATION OF LEAD TIME IN WEEKS. 

S 14 - OUT UF STOCK COST PER UNIT PER WEEK. 


& J 
J ; 
J 

** J 

J 

^APPLIES 

^❖APPLIES 


= 15 ■ 
= 16 
* 17 
= 18 
= 19 
ONLY 
ONLY 


- STOCK ON ORDER IN UNITS. 

- CARRYING COST THIS WEEK. 

- NUMBER OF UNITS BACKORDERED. 

- OUT OF STOCK COST THIS WEEK. 

- BATCH SIZE IN UNITS. 

TO PURCHASED ITEMS. 

TO MANUFACTURED ITEMS. 


US - NUMBER UF MINUTES PER SHIFT PER SUBPERIODS. 
a t t \ oit \ err) n ( T ) • P ( I ) « F { I ) « G ( I ) t H ( I ) t U ( I ) t &K ( I J 

ArSaYs’cONTAIN PARAMETERS FOR DEMAND GENERATOR FOR 
FINISHED GOODS (PUNCHED OUTPUT OF FORECASTER PROGRAM). 
T = STOCK NUMBER OF THE FINISHED GOOD* 

PINVTY - MAXIMUM AL LOWABLE INVENTORY IN RUPEES. 

^ an* arr ay*c gntain^ng^the^number^qf^shi ft-minutes p j=R 

SUBPERWD ThIt NS I NAS OPERATING AT THE END OF THE 
LAST PERIOD OF SIMULATION. 


I = WORK STATION NUMBER. 

NWS - NUMBER OF WEEKS OF SIMULATION. 

TINC - TIME INCREMENT IN MINUTES. 

GPNO - USERS SIMULATION RUN SERIAL NUMBER. 
ttd t M T T 1 A I TIME PERIOD IN MINUTES* 

Is I suitCH^TQ PUNCH INITIAL INVENTORY STATUS. 

IIS SW T ... T t rnST OF FINISHED PRODUCT. 

ACUNCO - ACTUAL UNI rQST PER PERIOD (AVERAGE). 

ACINVN - ACTUAL INVENTORY COS! KtK 


BL - A BLANK. 
DOT - A PERIOD. 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


J • 
IF 


PLT ' Rf L SP>-rn3ELY E ?HF T ?p?Tr» UNCHED IN t0Lums 3 thru 7 ARE 

asIerJsk anS’ T the Uim xl™ 6 LETTER U ’ THE LETTER *' AN 

** DESCRIPTION OF VARIABLES - 

VOI - INITIAL TOTAL INVENTORY VALUE IN RUPEES 
I F TP - LAST TIME PERIOD IN MINUTES. 

SUMRY { I » J ) — SUMMARY OF FINAL RESULTS 
I = NUMBER OF SIMULATION PERIOD. 

loJIr IND1CAT£S THE TYPE 0F INFORMATION. 

FORECAST FOR FINAL PRODUCT. 

ACTUAL DEMAND FOR FINAL PRODUCT. 

CALCULATED UNIT COST OF FINAL PRODUCT. 

ACTUAL UNIT COST OF FINAL PRODUCT. 

CALCULATED INVENTORY COST. 

ACTUAL INVENTORY COST. 

WEEKLY PRODUCTION TIME. 

INVENTORY ORDERING COST. 

CURRENT INVENTORY CARRYING COST. 

' CUMULATIVE OUT OF STOCK COST. 

* WEEKLY TOTAL PLANT COST. 

- CUMULATIVE PLANT COST. 

- CURRENT INVENTORY VALUE. 

■ CUMULATIVE PENALTY ASSESSMENT. 

CUMULATIVE COST OF OPERATION. 

TOTAL NUMBER OF SHIFT-MINUTES PER DAY OF ALL WORK 
STATIONS. 

TOTAL NUMBER OF SHIFT-MINUTES AVAILABLE EACH WEEK 
THE SYSTEM. 


CO DE 

„ n 

- JL - 
= 2 - 
= 3 - 
= 4 - 
= 5 - 
= 6 - 
: 7 _ 

= 8 - 
= 9 - 
= 10 ■ 
= 11 - 
= 12 • 
= 13 - 
= 14- 


J = 15 - 


IOTNMS 

TOTNMS 


TO 


** OTHER ROUTINES USED - 
**MAPA#* **FORCST** 


**PRODTN** 


**I NVENT** 


THE FOLLOWING SET OF INSTRUCTIONS REPRESENT THE FORTRAN SOURCE PRO 
WHICH IS USED FOR INITIALIZING THE DATA SETS. ITS REFERENCE NAME I 
PROSGO 

DIMENSION NMS (15) » A C 3 ) ,B(3),C(3),D(3),E(3),F(3),3(3),H(3),U<3), 

1AK ( 3 ) » PST (59,15,4) ,RMR ( 59, 5, 3 ) , OHR{ 15 ,6) ,HD(3) , FTAB(59, 19) 

i , # *** * * r * A t s r t i A Hi \/ S C n \ CTADO/KO) . DOnCO f 1AA . 1 \ . D f 1 R . 



,TA,US 

. .... ^ . . . . t t ’ai « 7 i ' ' i ■ » ’”»* • « 

INTEGER PST, RMR 
INTEGER RUN, TINC, ORDER, GPNO 
INTEGER P , BL , DOT , PLT 

DATA $UMRY/75»)*0./ cr nu AIn) r ., M ,<- r TnT TnM 

DATA 0,HRT,OC,SPNLT,CMU,CMC,VOH,PINVTY,SC,OH,AIDL,CUM,SSC,TDI ,TDM, 


CD LTi 



1VINVP, DUM , ITP ,1 X,TINC, US, IP,IFTP, IDL, ITTA,N , 1 50, J5D , J40, KTP/17*0 . » 

213*0/ 

data PST, RMR, Q» ORDER, HOLD, MW, NMS, FTAB ,OHR, APPLY, FTAB8 , RMC 1, OHRP, 
10HRM , A , B , C , D , fc, F,G,H,U,AK,SER,ABE/2515*3. 780*0./ 

THIS INPUTS THE NUMBER OF STOCK NUMBERS, THE NUMBER OF WORK 
STATIONS, THE MAXI MUM NUMBER OF COMPONENTS PER STOCK NUMBER, 

THE MAXIMUM QUEUE LENGTH, THE NUMBER OF DIFFERENT FINISHED GOODS, 
THE MAXIMUM PERMITTED NUMBER OF OUTSTANDING PURCHASE ORDERS, 

THE NUMBER UF WORKDAYS PER WEEK, THE NUMBER OF PLANNED WEEKS OF 
SIMULATION AND TWO SEE C NUMBERS FOR RANDOM NUMBER GENERATION 
RE AD ( 5, 15 H) I 50 , J5 0, J4C, LQ, NT, NOQ, I SP , IRQW, IP, I X 
153 FORMAT (0 15,211 J) 

READ INITIAL TIME PERIOD 
READ(5,124) ITP, IIS 
124 FORMAT (215) 

iftp=itp-i 

THIS READS ALL THE VALUES FOR THE VARIOUS DATA ARRAYS AND 
CONSTANTS USED IN OPERATING THE SIMULATOR. IC IS A CODE WHICH IS 
USED TO DETERMINE WHICH TYPE OF DATA HAS BEEN PUNCHED INTO A DATA 
CARD. I, J AND M RESPECTIVELY REPRESENT THE ROW, COLUMN AND PLANE 
SUBSCRIPTS OF THE MEMORY CELL INTO WHICH A DATA VALUE IS TO BE 
STORED (IF THE ITEM IS NOT PART OF AN ARRAY, THESE ARE LEFT BLANK) 
CM IS THE DATA VALUE WHICH IS BEING STORED. 

RE AD ( 5 , 12 7 ) IC,I,J,M,CH 
FURMAT ( 4 12, F 16. 6) 

IF ( IC.LE.O) GO TO 108 

(IC.LQ.i) PST ( I , J, M )=CH 
( IC.EQ.2) RMR( I , J, M ) = CH 
(IC.t'Q.B) OHR( I , J ) =CH 
( IC.EU.4) VQH=CH 
t)H=CH 
.’fQ .6) SC=CH 
, F Q . 7 ) FTAB ( I , J ) =CH 
i EQ . B ) US=CH 


106 

127 

107 


IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

CC 


(IC. 1:0.5) 


( IC 
{ IC 
( IC 
: CH 


I ND = IC 

IF (IND.EQ.ll) 
(INC.EQ.12 ) 
( IND.EU. 13 ) 
( IND.EQ.14) 
( IND.Eti.15) 
( IND.EQ.16) 
{ IND.EQ. 17 ) 
( IND.EQ. 18) 
( IND.EQ. 19) 
( IND.EQ.2U) 
( IND.EQ. 21 ) 
( IND.EQ. 22) 
( IC.EQ.15 
( IC.EQ.19 


IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

GO 


A ( I )=CC 
B ( I )=CC 
C { I )=CC 
D ( I )=CC 
E(I )=cc 

F ( I )=CC 
G( I )=CC 
H ( I ) -CC 
UC I )=CC 
AK( I )=CC 
PINVTY=CC 
VI N VP=CC 
AND.CC.LE.O. 

ano.cc.le.o. 


) E ( I ) — 1 
)U( I )=1 


( IC.EQ 
TO 106 
WRITE (6,22) 

WRITE ( 6, 23 ) P I NVTY 
WRITE (6,24) VINVP 


7.AND.CC.LE.0.)G(I)=1 



WRITE 
WRITE 
WRITE 

22 FORMAT 

23 FORMAT 

24 FORMAT 
IR Y* ) 

25 FORMAT 

26 FORMAT 

27 FORMAT 
If 8 CONTINUE 


( 6 . 25 ) 

( 6 . 26 ) 
( 6 , 27 ) 

( * ,'* 


OH 

VOH 

SC 


r4.;X,*START IMG PARAMETERS*/41x,19(lH*) 

/-tny.JkTMwCMtnnu , r,, , ^ . f A J 


(*0*,30X,*INVENTGRY LIMIT =*, F11.2) 
(*J*, 3 ;X, INVENTORY PENALTY =*, F ?.3,* 

( * , * »3v>X , *F IXED OVERHEAD =*,F10.2) 

( *0* » 30X » *VARI A 8L t OVERHEAD =*,F7.3,* 
(* )*»30X,*SHIFT CHANGE COST =*,F8.2) 


TIMES EXCESS INVENTO 


TIMES DIRECT LABOR* ) 


C *** CALL MAPA PROGRAM - TU GENERATE INVENTORY PARAMETERS. 

CALL MAPA (PST, FTAB, 1 5 G, J50 , US , I SP ) 

C THIS INPUTS THE NUMBER OF MINUTES PER DAY EACH WORK STATION WAS 

C OPERATING DURING THE * PREVIOUS WEEK*. EVEN THOUGH, IN REALITY, 

C THERE WAS NO WEEK PREV IOUS TO THE FIRST ONE TO BE SIMULATED BY THE 

C USER, A VALUE USUALLY 480) IS NECESSARY SO THAT SHIFT CHANGE COSTS 

C CAN BE CALCULATED SHOULD THE USER DECIDE TO OPERATE TWO OR THREE 
C SHIFTS ON HIS INITIAL RUN 

115 READ (5,162) ( NMS ( J ) , J =1 , J50 ) 

162 FORMAT (5113) 

C THE NEXT FOUR WRITE STATEMENTS GENERATE THE HEADINGS FOR AND PRINT 

C THE INITIAL STATUS (CARRYING COST, REORDER COST, QUANTITY ON HAND, 

C UNIT COST , DI SCOUNT ORDER QUANTITY, REGULAR ORDER PRICE, DISCOUNT 
C ORDER PRICE , AVERAGE LEAD TIME AND BACKORDER COST) OF EACH ITEM IN 
C INVENTORY. 

WRITE (6,500) 

500 FORMAT ( 1 HO , 5CX , *1 NI T I AL INVENTORY STATUS*) 

WRITE (6,5-H) 

501 FORMAT (1H , 7X»*ST0CK*,5X»*CARRYING*» 6X,*REORDER*,7X,*STOCK*,9X, 

1* UN I T*,6X,*DI SCOUNT*, 6 X,*REGULAR*,5X,*DI SCOUNT*, 6X,*AVERAGE*,5X, 

2* BACKORDER* ) 

WRITE (6,502) 

502 FURMAT (1H ,6X»*NUMBER*»7X,*C0ST*,9X» *COST* » 8X , *0N HAND*, 8X, *COST* 

1, 2X , *ORDER QUANTITY ORDER PRICE*, 2X, *ORDER PRICE*, 3X,*LEAD TIME*,< 

2, 6X , *COST *, / ) 


DO iOU 1 = i , 15.) 


I F ( FTAB ( 1,3) .LE.D. ) GO TO 100 

WRITE (6,504) I , FTAB ( I ,3 ) , FTAB ( I ,4) , FTAB ( 1 , 5 ) , FTAB( I , 6 ) , FTAB ( I , 9 ) , 
1FTAB(1, 10), FTAB (I, 11), FTAB ( 1 , 12 ) , FTAB ( I , 14) 

504 FORMAT (1H , 1 11 , FI 4. 5 , F13. 2 , Fi3 .0 , F 13.4, F13 .0 , 2F13.4, F13 . 2, F13 .4 ) 

C THIS SERVES AS A SWITCH TO PERMIT OR NOT PERMIT THE PUNCHING OF 

C THE INITIAL INVENTORY STATUS. SUCH PUNCHED OUTPUT MAY BE USEFUL. 

I F ( IIS.LE.O) GO TO IOC 

WRITE (7, 32 9) I , FTAB ( I ,3 ) , FTAB ( I ,4) , FTAB ( 1 , 5 ) , FTAB ( 1 , 6 ) , FTAB ( I ,9 ) , 
1FTAB ( I , ID ) , FT AB ( 1 , 11 ) , FT AB ( 1 , 12 ) , FT AB ( 1 , 14 ) 

329 FORMAT { I 3 , 2X , F8. 5 , F8 . 2, 2X » F8 .0, IX, F8 .4, F8. 0, 2F8.4, F8.2 ,F8.4 ) 

100 CONTINUE 
1001 V0I=0 

C THIS CALCULATES THE INITIAL INVENTORY VALUE (TOTAL) . 


DO 117 M= 1 , 1 50 

117 V0I=V0I + FTAB(M, 5)*FTAB (M,6) 

WRITE (6,156) VOI 

156 FORMAT ( 1H0 »42X ,*CURRENT TOTAL VALUE OF INVENTORY*, 13X, F10. 2 ) 
WRITE (6,295) 



nn,_D*h L u)1!r G SFT 0F 1 INSTRUCTIONS REPRESENT THE FORTRAN SOURCE 
^ HIEF [ USED F0R the PRODUCTION SIMULATIONS. ITS 
RfcFcRENCE NAME I S PROS IM. 


123 


41 


C *** 


READ (5,123) NWS ,T I NC , GPNO 
FORMAT (315) 

WRITE ( 5 , 4 1 ) IP, IX 

FORMAT (1H0,5X,*TH£ VALUES OF SEED FOR THIS SIMULATION RUN ARE AS 
IF OL LOWS* / / /iCX»*SEED FCR GENERATING LEAD TIME = *, 1 13//1GX,*SE£D 
2F0K GENERATING ACTUAL CEMAND =*,113//) 

DO 105 KKK = 1 » NWS 
CALL FORECAST PROGRAM 

CALL FORCST ( RMCl , NMS, APPLY, FTAB8 , FTAB,PST, RMR, OHR, Q, HOLD, MW , OHRP , 
13HRM, ALT, COST, ORDER, A, e, C, D, £, F, G, H, U , AK , SER , ABE , I5U, J 5 G , J 40 , W » NOO 
2, LQ , NT , I SP, FCST, ERR) 

CALL PRODUCTION PROGRAM 

CALL PRODTNM I5U»J5u»J40,LU,N00,NT) 


CALL INVENTORY PROGRAM. 


CALL I NVE NT { RMC a »NMS » A FPLY, FTAB8, FTAB »PST ,RMR ,OHR,Q, HOLD, MW, OHRP , 
13HRM, ALT, COST, ORDER, A, B,C,D, E,F,G,H,U, AK, SER, ABE, 153, J5C, J4G,W, 
2N0U,LU,NT, ISP.SUMRY, IROW,FCST) 

105 CONTINUE 

WRITE (6,296) 

FORMAT ( 1H ;,*SIMULATI0N COMPLETE*) 

DO lill M=1 , I 5u 
I F ( FT AB( M, 1 ) -LE.D. ) GC TO 1111 
WRITE (7,171) GPN0»M»FTAB(M,5) ,FTAB(M,15) 

FORMAT (215, 2 FI ''.0) 

CONTINUE 


295 


171 

1111 


THE FOLLOWING SET OF INSTRUCTIONS REPRESENT THE FORTRAN SOURCE 
PROGRAM WHICH IS USED FOR LISTING THE USER SUMMARY TABLES. 


I T 


REFERENCE NAME IS LSTTBL. 


THIS 

COST 

read 


53 


60 


READS ACTUAL UNIT CUST OF FINAL PRODUCT AND ACTUAL INVENTORY 
(AVERAGE ) . 

(5,5 0) ACUNCO , AC I NVN 

FORMAT (2F10.C) 

DO 60 1=1, NWS 
S UMRY ( I , 4 ) = ACUNCO 
SUMRY (1,6 ) =ACI NVN - 
CONTINUE 

PRINT SUMMARY TABLE. 

WRITE { 6 * 42 ) 

FnRMAT(lHl,4X,*FCST*,2X,*ACT*t3X,*CALC*»4X,*ACT*,5X»*CALC*,6X»*ACT 

1*,4X,*T0TAL*,5X,*INV*,5X,*INV*,4X,*0UT 0F*,4X, *PL ANT* »6X,* PLANT ) 

W R I T - ( 6 » 4 3 ) 

43 F OR MAT ( 1H * 4 X * *UNI T*» 1 X* *UNIT* * 3X* *UN IT*t 3X* *UNIT** 5Xf *1 NVhl** 5X t , 

1*1 N VN* , 4X , *PRODN*» 3X,*CRDER*,3X,*CARRY*,4X,*ST0CK*,6X,*C0ST*,7X, 
2*COST*,3X,*INV£NTORY*,3X,*TDTAL*) : 

WRTTFf6«44) 

44 F ORMAT ( 1H ,5X,*0NE*,2X ,*UN£*, 3X ,* C0ST *t3X,*C0S T *, 5X, *C0ST*,5X, 

1*COST*,5X,*TIME*,4X,*CCST*,4X,*COST*, 5X,*C0ST*,7X,*PER*,8X,*CUM , 


42 


c 


c 

c 


c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


45 FORMAT (1H0, I3,2F5.0,2F7. J»2F9.2, F9.0, 3F8.3,3F11.2,F8.2, F12.2 ) 

I OTNM S = 0 

TOTAL MIMUTES AVAILABLE IN THE SYSTEM PER DAY. 

DO o'. I=l»J5o 
IOTNM$=IOTNMS+NMS( I ) 

80 CONTINUE 

TOTAL MINUTES AVAILABLE PER WEEK. 

T0TNMS=I3TNMS*ISP 
READ PLOT CHARACTERS. 

READ (5,4,) BL» DOT , PLT 
40 FORMAT (7A1) 

DO 73 K=I » NWS 
DO 76 J=l,5 
76 T ( K» J ) = ;. 

ABSOLUTE DIFFERENCES BETWEEN THE FORECASTED AND ACTUAL VALUES FO 
FINAL PRODUCT IS ENLARGED ON PLOT. 

T (K, l ) = (ABS(SUMRY(K f U— SUMRY(K,2) )*6. ) + .5 

PERCENTAGE IMPROVEMENT IN UNIT COST OF FINAL PRODUCT IS ENLARGED 
T(K,2)=ABS(SUMRY(K f 3)-SUMRY(K,4)}*250./SUMRY{K,4)+.5 
PERCENTAGE IMPROVEMENT IN TOTAL INVENTORY COST OVER ACTUAL 
INVENTORY COST. 

T(K,3)=ABS($UMRY(K,5)-SUMRY(K,6) )*10j./SUMRY( K,6) +.5 

MINUTES OF PRODUCTION AND OPERATING COST ARE SCALED TO FIT INTO 

PLOT AREA. 

T(K,4)=$UMRY(K,7)* 50. /T0TNMS+.5 
T ( K , 5 ) = SUMRY ( K, 11) /b'Jti .+ .5 
73 CONTINUE 

THESE LOOPS DETERMINE THE POSITION FOR AND INSERT PLOT CHARACTER 


INTO THE PLOT ARRAY IN MEMORY. 

DO 77 J J= i » 5 

THE NEXT SEVEN STATEMENTS CLEAR THE PLOT ARRAY AND THEN INSERT 
DOTS FOR THE X AND Y AXIS. 

DO 70 K=1 » 50 
DO 7u J=l» 100 
73 P < K , J ) =BL 
DO 71 K=1 t 50 
DO 71 J = 3 » 100 » 8 

71 P(K,J)=DOT 

DO 72 J=l, 100 

72 P ( 50 , J ) =DOT 
DO 78 M=1,NW$ 

L=M*2 


K=5i-T(M, JJ ) 

IF(K.LE.0.0R.K.GT.50)GC TO 78 
P ( K, L ) =PL T{ J J ) 

L P=L + 1 

P (K,LP ) = PLT( JJ) 

78 CONTINUE 

WRITE (6,46) ( I, 1=1 » 5u , 2) 

45 FORMAT ( 1H1, 1 .*X ,*GRAP H IC AL REPRESENTATION OF 
11 5X ,25 14/ ) 

DO 79 K = 1 i 50 
J K = 5L-K 


RESULTS*/11X,35(1H 



ooooooooc^c'jct. 


47 MJKiiM l un , X A, 13 I , iU , A1 J 

79 CONTINUE 

WRITE (6,48) ( J,J = I,50,2 ) 

48 F ORM AT ( 1H » 14X, 251 4//40X , *N0. OF WEEKS * ) 

30 TO (*-0 » 40 i , 4U2 ,403 ,4 4 ) , J J 

400 WRITE( 6,14} ) 

141 FORMAT ( iHO , 17X, *PLQT CHARACTER F REPRESENT ABSOLUTE DIFFERENCE 
IBP. TWEEN FORECAST AND ACTUAL FOR FINAL PRODUCT*/ I 8 X,* VERTICAL SCAL 
2 IS 6 DOTS = 1 UNIT OR ONE INCH = 1 UNIT*//5GX»*FIG. - 11 *) 

SO TO 77 

401 WRITE (6,144) 

144 FORMAT ( 1 HO , 17X, *PLOT OF PERCENTAGE IMPROVEMENT IN UNIT COST OF 
IF INAL PRODUCT*/ 18X,*VERTICAL SCALE IS ONE DOT = 0.4 PERCENT OR ON 
2 INCH = 2.4 PERCENT*) 

SO TO 77 

402 WRITE ( 6 , j.45) 

145 FORMAT ( 1 H'V , 17X ,*PL0T CF PERCENTAGE IMPROVEMENT IN TOTAL INVENTOR 
1 COST OVER ACTUAL INV. C0ST*/18X , *VERTI CAL SCALE IS ONE DOT = l.C 
2 PCRCENT OR ONE INCH = 6 PERCENT*) 

SO TO 77 

403 WRITE (6,142) 

142 FORMAT ( 1H0 , 1 7X ,*P LOT CF PERCENTAGE UTILISATION OF PRODUCTION CAP 
1CITY*/18X ,* VERTICAL SCALd IS ONE DOT = 2 PERCENT OR ONE INCH = II 
2 P URGENT* ) 

GO TO 77 

404 WRITE (6, 14.; ) 

143 FORMAT (1H0 , 17X, *PLOT OF PLANT OPERATING COST. VERTICAL SCALE IS ( 
IE INCH =RS360v OR ONE COT =RS6G0*) 

77 CONTINUE 

STOP 

END 


C 

r *** THIS IS SUBROUTINE MAP A **************************************** 
*** PURPOSE - TO CALCULATE INVENTORY PARAMETERS. 

*** DESCRIPTION OF INPUT PARAMETERS - 

FDOQ - FACTOR OF DISCOUNT ORDER QUANTITY. 

IR - UNIT REQUIREMENT FOR MANUFACTURING FINAL ITEM (REARRANGE, 
- PUNCHED OUTPLT OF EXPLOSION PROGRAM). 

*** DESCRIPTION OF VARIABLES - 

LPT - LONGEST PROCESS TIME IN MINUTES. 

EPC - ESTIMATED PLANT CAPACITY IN PIECES PER WEEK. 

R - WEEKLY USAGE RATE. 

C EOQ - ECONOMIC ORDER QUANTITY IN NUMBER OF PIECES. 



<~> o o o 


c 


c 


DIMENSIUN FSIl 24, b, 4) , FTAB { 29,19) 

INTEGER PST, EPC » US 

lpt=;* 

CALCULATE LONGEST PROCESS TIME IN ASSEMBLY LINE. 

DO 20 1=1,150 
DO 20 J=1 , J50 

IF { L PT-P ST ( I , J , 1) ) ID, 2 ,20 
10 LPT = PST ( I , J , 1 ) 

20 CONTINUE 

WRITE (6,100) LPT 

CALCULATE ESTIMATED PLANT CAPACITY IN PIECES PER WEEK. 

EPC=US*ISP/LPT 

READ (5,200) FDOQ 

WRITE ( 6 , 110 ) EPC, FDOQ 

WRITE (6,120) 

CALCULATE DISCOUNT ORDER QUANTITY FOR EACH PURCHASED ITEM (FTABII, 

DO 30 1=2,150 

READ (5,213) IR 

I F ( FTAB ( 1,1) .EQ. 2.) GC TO 3u 

R=IR*£PC 

SE0Q=2.*R*FTAB( I , A ) /FT AB ( 1 , 3 ) 


30 

100 

1 ID 


EOQ = SQRT( SEOU) 

FTAB ( I , 9 ) =FDOQ*GOQ 

WRITE (6,130) I ,IR,R,ECQ,FTAB(I,9) 

FORMAT (ilU,5(/).33X, (LONGEST PROCESS TIME 14 , 5X , *»JNUTES*/ /// > 
FORMATdH ),20X, ^ESTIMATED PLANT CAPACITY =*,I8,5X,*PdCES PER WEE 
1*////20X, ^FACTOR OF DOQ =*, F10.4/// ) 

120 FORMAT ( 1H0 , 7X , * STOCK* , 7X , *UN IT*,7X,*WEEKLY*,7X, *eC GNOMIC 

1, #D I SCOUNT 0RDER*/6X,*NUMBER*,3X,*REQUIR£MENT*»2X»*USAGfc 

2* QUANTITY*, 9X,*QUANTIT Y*// ) 

L30 FORMAT (/7X, I5,3X,I8,F13.i,6X,F13.1,6X,F13.1) 

200 FORMAT ( F10. 4 ) 

2110 FORMAT ( 15 ) 

RETURN 
END 


ORDER*, 5) 
RATE*»7X. 


c 

c *## 
c *** 
c 


THIS IS 
PURPOSE 


iqd mil T MF FDR CST ^W*************************^****** 4 
TO GENERATE DEMAND OF FINAL PRODUCT AND TO COMPARE IT 

COMPARE IT WITH ACTUAL DEMAND. 


C 


*** 


ESCRIPTIOM OF INPUT P^^ETER - 

FCST(K) ~ USERS FORECAST FOR THIS 
ACTUAL SALES DATA FROM 
K = STOCK NUMBER OF FINISHED 


PERIOD FOR FINAL PRODUCT - 
PAST DEMAND HISTORY. 
PRODUCT. 



ooooooooo 


SUBR0UTIN1E MAPA ( PS T, FT AB , 1 50 , J50 , US, I SP ) 

DIMENSION PST ( 29, 8,4) , FTAB ( 29,19) 

INTEGER PST, EPC , US 
LPT = D 

C CALCULATE LONGEST PROCESS TIME IN ASSEMBLY LINE. 

DO 20 1=1,150 
DO 20 J=1 , J50 

IF ( LPT-P5T ( I , J , 1) ) 10,2 ,20 
ID L PT=PST ( I , J , 1 ) 

20 CONTINUE 

WRITE (6,100) LPT 

C CALCULATE ESTIMATED PLANT CAPACITY IN PIECES PER WEEK. 

EPC=US*ISP/LPT 
READ (5,200) FDOQ 
WRITE (6,110) EPC, FDOQ 
WRITE (6,120) 

CALCULATE DISCOUNT ORDER QUANTITY FOR EACH PURCHASED ITEM (FTAB 

DO SO 1=2,150 

READ (5,21.:) IR 

I F ( FTAB( 1,15 .EQ.2. ) GC TO 3U 

R=IR*EPC 

SE0Q=2.*R*FTAB( I , 4 ) / FT AB ( I , 3 ) 

EQQ = SCRT( SEOQ ) 

FTAB ( I ,9) =FDOQ*EOQ 

WRITE (6,130) I,IR,R,ECQ,FTAB(I,9) 

30 CONTINUE 

10D FORMAT (\Hj.,5(/),3DX,*L0NGEST PROCESS TIME = *, 1 6 , 5X , *M INUTES*/ //, 
1 ID FORMAT ( IH ), 2 OX,* ESTIMATED PLANT CAPACITY =* , I 8 , 5X , *P IECES PER Wl 
1*////2QX, ^FACTOR OF DOQ =*,F10.4///5 
12D FORMAT ( 1H0 , 7X ,* STOCK*, 7X , *UN IT*,7X,*WEEKLY*,7X,*£C0NQMIC ORDER* 
1, *D I SCQUNT 0RDER*/6X,*NUMBER*,3X,*REQUIREMENT*,2X,*USAGE RATE* 
2* QUANTITY *, 9X , *QUANT I T Y*// ) 

130 FORMAT (/7X,I5,2X,I8,F13.i,6X,F13.1,6X,F13.1) 

200 FORMAT ( F10. 4 ) 

21D FORMAT ( 15 ) 

RETURN 

END 


*** THIS IS SUBROUTINE FORCST ************************************** 
*** PURPOSE - TO GENERATE DEMAND OF FINAL PRODUCT AND TO COMPARE IT 
COMPARE IT WITH ACTUAL DEMAND. 

*** DESCRIPTION OF INPUT PARAMETER - 

FCST(K) - USERS FORECAST FOR THIS PERIOD FOR FINAL PRODUCT - 
ACTUAL SALES DATA FROM PAST DEMAND HISTORY. I 

K = STOCK NUMBER OF FINISHED PRODUCT. 


oooooooooooooooooooonooo 


*** DESCRIPTION OF VARIABLES - 

W ( K ) - CALCULATED DEMAND VALUE OF ITEM K. 

ERR { K ) - ERROR IN FORECASTING FOR ITEM K. 

SER(K) - SUM OF ERROR IN FORECASTING (CUMULATIVE) . 

ABE i K ) - SUM OF ABSOLUTE ERROR IN FORECASTING (CUMULATIVE). 
SALES - OUTSIDE SALES (IF ANY). 

NVAS - STOCK NUMBER OF RAW MATERIAL FIRST IN NETWORK QUEUE. 
PNU - STOCK NUMBER CF PURCHASED ITEMS FOR WHICH PURCHASE ORC 
IS PLACED. 

SIZO - QUANTITY OF ITEM PNO ORDERED IN THIS PERIOD. 

ORDER{ J,K ) - ORDER MATRIX FOR PURCHASE ORDERS. 

J = ORDER NUMBER. 

K = A CODE THAT INDICATES TYPE 6F INFORMATION. 

IF K = 1 - STOCK NUMBER OF ITEM ORDERED. 

K = 2 - QUANTITY ORDtREO FOR THE ITEM. 

K = 3 - TIME PERIOD THIS ORDER WILL ARRIVE. 

ALT ( J ) - TIME PERIOD WHEN ORDER WILL ARRIVE. 

J = ORDER NUMBER. 

JOC - CUMULATIVE COST OF ORDERING IN CURRENT PERIOD. 

*** ROUTINE USED - 
**RANNUM** 

SUBROUTINE FORCST ( RMCi *NrfS, APPLY, FT AB8 , FT AS, P ST , RMR, OHR , Q, HOLD, 
13HRP,UHRM, ALT, COST, ORDER, A, B,C,D,E,F, G, H, U, AK , SER , ABE , 1 50, J5Q , J 
2W , NOO , LQ , NT , I SP » FCST, E RR ) 

COMMON /CLCJR/ 0 , I TP , I X, T I NC , HRT , OC , S PNLT ,US, CMU, CMC , KTP , I P , VDH 
IP INVTY, SC, OH, AIDL, CUM, SSC , TD I , TDM, I FT P , I DL , I TTA, V I NVP , N , DUM 
DIMENSION RMC* { 29 ) ,NMS( 8),APPLY( 29) ,FTAB8( 29),FTAB< 29,19) 
1PST ( 29, 8, 4 ), RMR { 29, 5,3 ) ,OHR{ 8 , 6 ) ,Q( 8,23 ,3 ) ,HOLD( 8,4 

2MW( 8,4), OHRPt 29),0HRM( 29 ) , A LT ( 200 ) , COST ( 29 ) , ORDER ( 200 , 3 ) , 
3A( 3) ,B( 3) ,C ( 3 ) , D ( 3 ),E( 3) ,F{ 3),G( 3) ,H( 3),U( 3) ,AK( 3) , 
4SER( 3) , ABE ( 3),W( 3) , FCST( 3 ) , ERR { 3 ) 

INTEGER TA, US , PN, QT »TM ,PuU , HQ , T , HOLD, SN, Q , SUT , TEMPI , TEMP2 
INTEGER PST, RMR 
INTEGER T INC, ORDER 

DATA TEM?2,TA,PN,QT,TM ,PjU,HQ,T ,SN, SUT, TEMPI, DUMB / l 1*0, 0 ./ 

DATA NP, S ALES , PNO, SIZO /C', 3 * 3 . / 

DO 50 ICC 2=1, NT 
40 = CST ( I CC2 ) =0 
53 ERR ( ICC2) = ' 

C THE NEXT 6 STATEMENTS WRITE HEADER FOR FORECAST DATA 

WRITE (6,135) 

WRITE (6, 136) 

WRITE (6,137) 

117 DO 118 K=i , NT 
C READ USERS FORECASTS 

63 READ (5,138) FCST(K) 

CALL RANNUMK I X, RRAN ,RA AN ) 

NP=K 

3EV=(-2.*AL0G(RRAN) )**C. 5*C0S ( 6. 283*RAAN ) *AK( NP ) 

V = ITP 

C CALCULATE DEMAND VALUE FOR CURRENT PERIOD. 

W ( K ) -AC NP)+B(NP) *V+C (NP ) *V**2+D ( NP ) *COS ( ( 3. 141 6*V ) / E ( NP ) ) +F ( N! 



o o 


1CUS C U. i.H , *.*V }/G('JP3 ) + M ;P)*COSC ( 3. 1416* V ) /U( NP ) ) +D£V 
I w = W ( K } + . j 
W i K ) = I h 

r RR { K 3 = FC G M K 3 - W ( K ) 

A bl: ( K 3 = AS _ ( K 3 + A C S ( t R R { K 3 3 
SjR(KJ=SZFt(KJ+'f<R(K» 

W ^ I TF U>, 9) K, ITP,a( K) , FCST { K 3 , ERR ( K 3 , S £R { K 3 , ABE < K 3 

IF { W ( K ) . L . 0 . 3 W ( K J = D 

118 CONTI )Lr 
X YZ= . 

119 cr^'Tir-Uc 

C PUI fj p=u 2 r NO OUT SIC: ■», L7-S. 

M P= 

IF ( r!P . LF • J Of TO 12 

IF (SfiLcS.OT. FT^.fcC mP, 5 >) SALES= FT AB { i\ P t 5 3 
FT A- ( JP, 7 ) = FT --'it NP, F ) + SAL:_S 
FiOr,( i F , b 3 = F T ;P, 53-SALES 

XVZ = . . 

WRIT” {fij .■ : ) , PjSALiS 

o; to i;. 3 

121 CluTUlr 

WRIT1 ( S, - ) 

YY=0 

3t\=?AT 7,0 CRT 2RS TO REET CALCULATED DEMAND IN THIS 

P R I Of: . 

K TP = KTPu 
I K = U 
I J=- 

:< ; 4 if ( i k . ; o . j 4 ) : j= i j+: 

IF ( I J.OT .15 3 SO TO 2 09 

IF UK.b'.JAu) ' K=<: 

I t- = i K + 1 

N V / S=RMR( I J, I K, .. ) 

IF { (FT; L { . ' V ' j ? _ ) . - U . ? . 3 . OR . ( NVA S . L E . 0 3 3 GO TO 204 

Pi Ll-i "V A S 

SIZij = FLlb-{Rr^UJ, IK, 2 })*W( 1 3 
C P <- R T liL'-T; !, IS CHECK-C T j INSURE IT IS A VALID NUMBER. 

DU / 5 1 1 = . j I 

IF ( I F I /. ( ? • 1 3 .. U.I I ) GC ; 0 205 
2,5 COOT Lit 

WRITE ( 6 j US) P • 0 
Y Y = 

I F ( IJ.OT.Taw) GO TO 2C9 

GO TO 2 4 

150 WRIT! { Gj - 3 P .*0 

300 F OR!-) AT ( 1H‘>, 5/» *QU4NT I TY ORDERED FOR SN *,F5.0 f * LESS THAN ZcRt 

1 ORDER IG.'TRl-D*) 

yy=: 

I F ( IJ.GT.i5' 3 GO TO 209 
GO TC 20 m- 

206 I F (FTA'it II, 13..U.2.) GO TO 2- 0 
I F (SIZO.LT. ..) GO TO 15 

17 1 2 7 J= 1 ? 2 . 

IF (CRCcRt J, i J.'b.O) GC "sO 207 
C , L L It A' ."IJ' 4 ( I P j RRAN » R A AN 3 

3=(-2.*A‘ OG { IRAN 3 3 . 5*C0S ( 6 . 233*R AAM 3 
C C A L u U L A T f ” I u (L-A!J Tlf^E) WHEN THE ORDER J WILL ARRIVE 

I T I I \ - £_T I 7 ‘ > \ -lCT A D / . T 1 1 I *C jl C 1 fit. T ( V Tit \ 


UNITS 


1RP •>’ ( J l = SI 2. : 

0 P 0 R ( J , -• ) = A L 7 ( • i ) 

FV C ( I It - : ) = F7‘.'_ U I , iS )+iIZCJ 
<K= IF I X ( p ! ■) 

3C = 'jC+F7Ar( KK, ) 

] PA J. = PU I 
I $IZ*=SIZ J 

n 'Li C 6 , ■ ' » ) l PiJQ„ , i $ IZ ^ » J 

>0: F OR '' lA 1 ( i H 7 , „ ,',*JROER PL ACZO FOR S'. *, 15,* FOR*, 112,* 

iORD-R N'J *,I51 
IF (IJ.GF.I5 ) Go T'J L‘,9 
GO TU 20*+ 

207 CONTI NU’- 

WRITL: (6,2-}+} P.fU 

YY = j, 

If ( I J.0T.I50 ) GO T 3 2C9 
SO T 0 2 •'+ 

£ 8 . 1 1 I f . ( 6 , 5 ) PIG 

Y Y= . 

IF ( i J.bT.lBu) GO TO 2C9 


W* J i t 4 4 ^ 

I F <YY. . } .CITE (6,2 *'6) 
3 T P = 1 T P + 


in 

ccm 1U 

. 




R ’TOR j 




122 

Fl-R.IAl 

! 

* . 

i»> 

123 

F OR > ia r 

c 

f 

>. ) 

1 :: 4 

F0R5AT 

( I 

A ) 


125 

FOR 1AT 

( I 

) 


1-6 

FOR te f 

( i 

.. , 

. 1 T ' -1 ) 

i 27 

FORMAT 

(4 

I ■ 

, F„ . ) 

1£B 

F.-pvlA T 

(2 

i / 

-A . 

,F . ) 

i29 

Fur: .at 

( I 

f 

- » F j • 5 , 

130 

FiJR" 1 a : 

( 5 

* 

) 

* 31 

F OR MAT 

(i 


» - A, 4 0 


132 FUR UT (-0 , /. »*TI Mrr* »£ . ) 

A 35 FORMAT (9,., I 5, 6 .,4( „0F7. ./: OX) ) 
j.34 FliPOO-.T C* ;, I ,, :C,4(‘.0F7.. /iOX) ) 
1)5 FORMA f (-.0,5 X , *R iSULTS OF FOF 


FORMAT (;. I ,, j,( f !+{ ‘.GF7. . /iOX) ) 

FORMAS’ ( _ 0 . ,5 ;'.,*R:GULTS OF FOREC AST*/51 X , 19 { 1H* ) ) 

FOR V- T (.Mi f .. 0 , *S TUCK * , o'A ,*T I ME* , l&X, *4CTUAL* , 6 X , *FORECAST* , 


I2 0X,*GUF. ’ i F * , 0 ,* 5 UM OF aBSOLUTt*) 
i.27 FORMaT Co *MUyB ; - R* , 9X , *PER I OD*, 8X, *V4L UR*, SX, * VALUE *, 9X, 


). * ; R . \ U R * , o / , * RRoR* , RX , *£ RKOR* ) 
123 F-R :*■'•;• T ( F . . ) 

*39 format ( „io, j) 7 ,1'h. i .o»,I5,5X,f. 


i39 FoR. -1 A I („ HO, ~u 7 , 14 , 102,15, 5X,F10.2,3X»FI0.2»6X,F3.2,3X,F10.2»3X, 

*F A (j . — ) 

143 FORMAT CM , 4 JY , *A DD I T 10 iAL SAL<=S*) 

141 FORMAT ( I . , F 4 

142 FORMAT ( i, Hi 1 , i \>A , *S T CCK NUMBER*, I5,5X,*SALtS OF *,F10.0,* UNITS*) 
aA 3 FORMAT Cii , - X,*NO ADDITIONAL SALES*) 

22s) FORMAT (111;, 5 X , *P UKCH AS ORDERS AND ERR0RS*/51X, 26 ( iH* ) ) 

2 format cf. ■ . ) 

224 FORM 4 r CHU, v;,*UMADL5 TL FILL ORDER FOR STOCK NUMBER*, F6. 2 , 5X, * 
-0RU..R BANKS FULL*) 

223 FUR MAT (iH , C\,*STUCK NUMBER* ,F 6 . 2 , 5X,*DDES NOT EXIST*) 

2 >5 FORMAT ( iH , NX , *5T0CK NUMBER* , F6.2 , * IS HOT A PURCHASED PART*) 
236 FORMAT ( 1 HD, iC'., *ALL REQUESTS PROCESSED*) 
t N D 



suB^nunn*. p r< 


TN( i 5 » J5 » Ja i » LO » 30 » NT ) 


c ### 

THIS IS SiF-RUO 

I ’• i. PR LOT 

c ❖ 

PURPOSE - T 0 S 

i :‘ULAT. C 

C 

'iitR K 

;TUIU' . 

C 

t ** 

D ..SCRIPT! 0 J ijr 

VAR i A EL IS 

c 

U { I » J » K ) - 

•laird: CCS 

c 


i- TURK ST 

c 

I - ,!~K/ 

j ! h i i 1 j ! L > ■ 

c 

J = P’-J' 

POSITION 

c 

h = \ CD i 

. 7-i A >’ I NO 

c 

IF K = .. 

- st.jCk nu 

c 

K = _ 

- S’JJ.'TITY 

c 

j j. 

- UR. R N't 

c 

( I ? J ) - • 

S\R S T A" I.J 

c 

I = ' r 'RK 

Si d If'., L. 

c 

j = CD.. 

.. TH a t ! hi 

c 

IF J = . 

- STOCK NO 

c 

J - _ 

- II TF; 


££K OF PRO DUCT I Oil ACTIVITIES AT EACH 


C 

C 

c 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

0 

c 

c 

c 

c 


J = 
J = 


=R CONSIDERATION. 

JN R CD.MS I DCRmT IU'J * 

.CATl-S TYPE uF INFORMATION. 

Jek of THE 1 T. M » 

Is. Trie QUEUE. 

r r.H .b::< op th:- lo t . 

I,j l OFT AILS, 
t ( .,rP, CONSIDERATION. 

NS t CATES TYPO UF INFORMATION. 

TOOK NUMBER uF WORK IN PROGRESS. 

F; WJPK STATION WILL BE THROUGH PROCESSING. 
- GRO'R NU' BE R JF I Tfc ^ L iDER PROCESSING. 


J - 
IF J 


= -PORK 
- . C" 1 


- j R 0 ' 

R NU' 

BE 

- LOT 

s: l: 

OF 

• J .'LU 

ULCC 

V 0 

s at i i 

! *! L ' 

.'.SR 

: T H ; T I NU 

A Cf'. 


J - 

J = 


J - 


TA ( I ) - 
ILT(T) 
usd ) - 


INDICATES TYPE uF INFORMATION. 

- STUCK NU T':R ;jF I T i I HULL BLUCK. 

- i *U \ -:T I I Y I I HOLD BLOCK. 

- u'k.TITY HOLD FLOCK IS TO ACCUMULATE BEFORE 
R_i_<. ASL.-j GOODS TG THE NEXT WORK STATION OR 
I tV- TCRY. 

, - uRG.-R NUNB.R OF 17-M IN HOLD BLOCK. 

XI '■ THE WORK STATION I IS TO WORK IN MINUTES PER WEEK. 
F .0 OVERTIME THE WORK STATION I IS TO WORK. 

U 7- ,r jF SFI PI-MINUTES THE WORK STATION I CAN WORK 

. day. 

i -7 ■» ( I ) * ’ ' ' t 1 II) - T-KP'RARY DUPLICATE PREPARED FOR TIME 

t V AI CABLE AND FORCED OVERTIME VALUES FOR WORK 
ST AT * Uf 4 1. 

qM(K) - UP': I sTA T l ON jUH6i_RS AT WFIICH A COMPONENT ENTERS THfc 
PR , C S S J WK-.Ri- K = COMPONENT SSQUENCe NUMBER IN 
> 4 V 'ING aT WORK STATION. 

SK ( I ) j 1 r II. t ) - 71 MPCR-RY ARRAYS CON T A IN ING THE SH AND QUANTITY 
i if C CMi'ONtr NTS IN SEQUENCE NUMBER K. 

SUTC7 ( I ) - ,;::UP-R OF SETUPS OF WORK STATION 1. 

RMC1(1) - i. L'ST UF RAW MATERIAL I > USiD. 

•jHrP(I ) ,UU l.’tl ) - NAN AND MACHINE COSTS FOR TH-r APPROPRIATE 

SHIFT FOR ITEM I. 

1 L L £ ( I ) - 10LL TI I'U- A f WORK STATION I (TOTAL THIS WEEK). 
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HM { I } , ;i w J m - IDLE YACHINE AND (4 AN COSTS AT WORK STATION I. 
COST ( J ) - CUMULATIVE COST OF ORDERING AND PAYMENTS OF THE ORDER 
l.J RUPEES FCR ITEM J. 


* AUUTI.iES JS~C - 
**.,RIT o ** 


0:?-, 'SION NMS ( >,a£3) ,B{3),C(3),D(3),L{3),F{3),G(3),H(3),U£3), 

1A'<(3) ,PST( '9, 3,4) ,RMR (29, 5,3),OH«»( 8,6) , FTAB(29,I9) 

-'*W< 3) ,S-Xr . ) , -'B; (3 ) ,PPPLY(29) ,FTA8S ( 29 ) , ORDFRt 2„S ,3 5 ,Q( 8,2C,2) , 
*■< 1 L r ‘ ( 0 r A ) , M W ( ' ,'*) , CHRP (29) , OHR M { 29 ) , ALT ( 200 ) , COST ( 29 ) , RMC i £ 29 ) 
■tSUTCTf ) , T A ( 5 ) , I OT ( S ) , DM ( 5 ) , I OLE { 8),OMI( 3),OHJ( 8 ) ,OMJ ( 8), 
51 P I C ( 3 ) , S K ( 5 ) , I NTC T £ S{ 8) ,MIOT{ 8),MTA< 8) 

C'J'SlOH /C lC J R / G,T TP, I X, T I NC , HRT ,OC , SPNLT , US , CMU ,CMC , KTP , I P, VOH, 
I JVTY,SC ,0H, AIDL, CUM , SS u ,TD I , TOM, I FTP, IDL, I TT A, V I NVP , N , DUM 
CL ' 1 IN /PROS I: / R M C‘ , M NS, APPLY, FTAB3, FTA6, PST, RMR»OHR,Q, HOLD, MW, 
t O.-iR P , f F VI , ALT, COST , ORDER, A, B, C» D, 2 , F, G» H, U, AK , SER, ABE, W,SUTCT,TA 
El Cl, j , ,jL , -J ; I ,OH J,UMJ, IDIC,SK, IMTCT,MS,ISP,MIOT,MTA 
I • T -S .R T. , u : , P QT, TM ,Pi I J , HQ , T , HOLD, SN, Q , SUT, TEMP i » T EMP2 »GM, SK» 
J , S ! T C ’ , T 1' 1C, ORDER 
J « ! _ j „ R P ^ ' f ~ ' R 

C - < L* n-\r>'V,f PM,OT,T M ,PNU»HG,T,SN»SUT, TEMPI, TEMP 2, TOC , I SUM , DE 
U..T P\, "7, - ',p'.U, ML,”,S k ,SUT,T^Mp A ,TEMP2,I SUM , TOC, D3/ 11*0, 2*0./ 


>U~,K * K , JbH , K6H ,T , I STRT , I E ’JD , I BY , I DA Y, PNU ,MNU > HQ/ 13* 3/ 


J H, ’ ] S = 


, J5 


, 


, 


t •- ( jh'-’: }- 


IfT£ JUY i ) =■■ 
ILL (JFi )= 
u.'KJ'-w ) = 

T'JtJHM )=>.: 

i 1 1 c r ( j Mj) = j 


KS( JH •.?) = 

1235 SUTC rtJHI )= 

C-u 2-6 J _ ,152 
C OS T { Jh'S ) = u . 

UhP.P( JF *••> = . 

1236 * AIR:- 1 ( JF ’ j ) = . 

DC 4. J~ JH V > , J-rO 

nr.t jh iz )=■ j 

1 Lie ( JF «13 ) = 


1237 S K ( J F ' ! 3 ) = 

0 . 1 U = 0 . 

A I 0 L = J . 

D-~ . 

HAT = ' . 

It Mp = 0 

lc: =io- 

WRITi (6, 

X =• 

C GLNuKAT . P'*.nDIX“IOi'j uRCaRS FUR CURRENT WEEK TO MEET THE DEMAND OF 

c FINAL PliOJUCT. 

0 PN - STUCK NUPFI.-R. 

C QT - uLAUTITY Til Bfc PRCDUCED. 

p 1 ;= 

QT = P( PN ) 



onoo ooo n o o 


Gt Ti_ 9 

lor i f ?■■= j 

I F( P I.GT. I-v ) GC I'.' _ . 8 
wT = H { „ 3 

99 IF { P :.Lt . )) Gri T j ;00 
C CHL Cf' FUP. .fJW-./ISl 4KT P-IODUCT NO. 

IF ( P, ! . LT . . tj ; . . 3 1. GT. I £ ) GO TO i 4 
C PET^RM I it Fir* ST WS In NETWORK FOR THE ORDER. 

DO , C'* 1 'J = . ? J 60 
*1 0 = PS T { P'l » 13, - ) 

IF (mo. u. .) GO TO 5 
i 3 C 1 4 T I * 3 1 . 

C P Ri' DUCT HAS NO iwITIAL WS. 

WRIT, { 5, 65 ) P j >OT 

pA; = . 3 i-! + l 

3 U To 1 < ) 7 

- 4- WRIT, (A, _ ,6) PM »QT 
X = 1 

P< ' = 0^ + ^ 

GO t,i 1C,. 

S. i r ;Ci! At- 1 !; FOR F I PS T AvAILaELE EMPTY POSITION. 

1-5 [•,. , L w 

i r ( <( I _ , - •} _ ) • _ w • U ) GC i 0 1 uT 

loc c* j ‘t 

C V- -nF LC 1 - ecus : S THIS PRODUCTION ORDER TO Be PURGED FROM 

Th_ f Y S 7 l ’ . 

W K I , ( * , . j ) P j » 0 T 

P:' =- 5 ..+ 1 
3 ' i It. 1 - 

P LAC A CKL, I-- PROPER CU -Ur If, FIRST AVAILABLE EMPTY POSITION. 

107 G ( I - ,14?', )=P0 

N = , i+ 

C'i/.riiTY 1 Trif UlJ£IJc = QUANTITY to BE PRODUCFD/LOT SIZE. 

0 ( i ■> » 14 j 3 ) =GT/ 1 FI X ( FTA E ( PN, 2 ) ) 

3 ( 1 . • , 1 4 , 3 ) =N 
WRIT _ (b,L5) °'t , QT 

75 F uP. !■ M (;H , _ Y, BORDER PLACED FOR SN*,I4,* FOR *,I5, * BATCHES* 3 
PA=PM+1 
3 LI T i 1 

108 Cf'NTil'b. 

IF (X. r :Q. . ) WRITE (6,17c) 

ALL P R 0 1 L 3 S A rt L r ORDERS PLACED IN INITIAL QUEUES. 

*** G7?!tR AT T I M E EACH WORK STATION IS TO WORK. 

1' <9 CCNTINL;.- 

I F (K4K.L-. >) GO T 0 111 

1 MK4K5=JBH 
I fiT { KAK )=KCH 

I F (‘GiK.LT. . 0R.K4K.GT .JdG) GO TO 112 
C *** TIT, ALLOU-n is BROKeN INTO ONE PAY SEGMENTS. 



T A { K^K)=T 1 { K- \ ) /l CP 
i C'T (KiK >=! ]T(K- K5/ISP 

C *** ALL 1 S \U. C • i _ C K ; 0 TO I iSURt THAT THEY ARE ACCEPTABLE TIME. 
LHH=: *LS 

if (t:- (« v ) .ot.l-.h) go t., 3 : . 

IF IT A (KrK ) . :;c. L3H) h: S (K K)=.,*US 

IF { T 1 (K-+K) . i Lt H .AND .T ,(K4K) .GE. ( 2*US } ) MSIK4K )=2*US 
IF { TA ( <4-K ) . i \ . ( L*US ) ) :i:(KhK) = 2*US 

IF ( T A ( K4 K ) . L 7 . U*;JS > . AUG .T A ( K4K ) .GE.LS) MS(K4K) =US 
I F { T A ( K>K ) . L - . ;$) MS ( K6X)=US 
J F ( T a ( K4i\ ) .LT. JS) W R I T ( 6 , i 71 5 K4K 
I F (TA(K4K) .LT.'JS) ta ( K-,K)=US 
GC TO 1 ,9 

r *** EXECUTED WH£.fJ TIME AVAIL ‘.BLE VALU" IS GREATER THAN THREE SHIFTS. 
1x0 WRIT:-; C E> , S TO ) K4K 
TA(K4<)=3*US 
M S ( K 4K ) =3 -< t: J S 
G 0 7 J 1 C 9 
ill COM' I MU 5 

DC y . L I K f = i , J 9 
MTA(LIKS ) =T A { LI K_ ) 

MIGT(LIKE) = I07( LIKE ) 

2ui c,'i;timu- ; 

G 0 T C> 1 3 

C *** E R ?. ' R '• S , G7 C J\ EX. ST AMT WORK STATION. 

li 2 WRITE {£}.,:;) K t* 

/ = 

TA(KmO=U. 

’■'S (:<-■:< )='j s 
Go TO I >V 

113 CLNTI-it ' 

C *** CHcCKSNG FOR 2 '-31 • .G 71'. AV mIL AGL z DECISIONS. THGSu MISSING ARE 
C ASSIGiCG "HE V- CUE OF CU_ SHIFT. 

DO . >v •= , ,j L t 

IF ( T A { I ) . Vf . } bO T'j i: 

IF (OHRIIj, ) . L - . .5 G; T- 

r a { i )= is 

M T A ( I )=US 
«S (I ) =L S 

1 14 CONTI vL. 

C PRIJT CUT ALL T i >5 Ai/AIL-BLES AFTER ANY NECESSARY CHANGES. 

00 5 6*6 I=.',J r -J 
I F ( T A ( I ) . L . . i GO TU 66:6 
WRITE to, C& d7 ) I , TA (I) 

6667 FORMAT { ii-iO , . OX , *T I M l AVAILABLE AT WS *,1-,* IS *,110) 

6666 conn. ML. 

C READ PAP.A’i T RS FOR IN TERMED! ATb RESULTS. 

READ (5,i7-) IMflSTRT, IE,,D» I RY, I DaY 
C COMPUTE M'.xriUM TIME AVAILABLE. 

K B H = C 

DO 3166 J BH= i , J 3 

3i66 IF { TA ( JBti) .GT.KBH) KG F=‘l A ( J BH ) 

T']=K8H 
K = 0 

C ** Gr iiRAT HILO QUANTITY. 



o o o o o o o o o ooo ooo oo oooo 


IP ( PHL.l . ■) ji, TU Zt-L 

x= . 

WRI'IE (ot.tipi’ ) ° JU , "i'iU ,1 i -i 

bt(S FOR’ / 1 ( 1 i ! > »_'/.» *HDLD CU7NTITY FOR STOCK NUM3ER *,I5,* AT WORK STA 

3 T I 0 ! Chi . G ^ F j TU *, 15, * BATCHES*) 

H Q-tU / I F I U FT ' B ( P\! U » i ) ) 

IF (HJ.Lc.U) i-.»=_ 

IF (HOLOd.'I’J, l.r'J, PNU ) *-'OLD ( VNUt 3 ) =HC 
P ST ( P ) = 'iu 

3 •' i 7u I..5 
250 CONTINUE 

** START SIMULATION. 

JR5 INC* X STC »S THRU The 6 DAYS OF THE WEEK. 

DO 20 > JRN=. , ISP 

UK; INDEX ST PS THRU 7H WORK STATIONS TO RE-ESTABLISH THE 
INITIAL T I ML VALUES. 

DU 202 UK „ = i, J50 
TA(LIKE)=UTA( LIKU 
I U7 ( L I K •" } = M I OT ( L 1 K Z ) 

202 CONTINUE 

START THE ■ 'AS I, R CLOCK A j* A MINUS *TI NC* VALUE Bj CAUSE A GO TO 
LuOP IS US'D fu'l I sDUI/jG THL MASTER CLOCK. ON THE INITIAL PASS, 

7!U : X T '•-,7 - i ,7 CAUSES THE INITIAL TIME TQ BE ZERO. 

L:.*> T = { -T INS ) 

12 7 T = T+T J NC 

JR i.'iCS/. ST PS Tipui ALL wORK STATIONS TWICE EACH TIME INCREMENT TO 
PR-VUT IjL 71 . -CALS . THc. WORK FLOW NETWORKS ARE NOT IN WORK 

ST AT luN S 2u‘J ‘<Z : . 

Du ’-A JR= , " 

D ' J . * -1 I = , , ,J 
K KK S = 

SUT = 0 

C i 1 _SK S TO o- l i* WORK STATION I HAS REACHeD ITS FINISH TIMe. 

IF ( T.G7. “X( 1 ) ) GO TO 14 > 

IF PR..C .Sbi 1G ; IT FINISH 0, IT GUrS TO CHECK THE NEXT WORK STATION 
I F ( "-W( i, I) .GT.l ) 30 TC 4 -> 

IF TRJ £ , WORK VAT I ON IS IDLE AND IT IS TIME TO LOOK FOR A 
PRODUCT 7) PROS SS. THIS IS DONE AT STATEMENT ^3 
I F { H»< ( I , I).. J. .7) DO TC ,3 . 

REACHING THIS i) .MOTES COMPLETION OF A PROCESSING AT THIS WORK 
S T A i ION. F T HI 3 IS THE FIRST CHECK AT THIS TIME INCREMENT, THE 
NUT 3 2VUUL UaTE.M ENT $ ARE USED TO MOVE THE PRODUCT OFF THE WORK 
STU I UN. IF THIS I S THE sECuND CHECK, THE PRODUCT WAS MOVED 
eARLIER A! u I r IS TIMe TO LOOK FOR A PRODUCT TO PROCESS. 

IF (JR.DO.P) GO TO j 3 - 

IF ( HOLD ( 1,15 . A • ' • AND .MW( I,» 4 ) . EQ.* ■ ) GO TO 130 

I F (HULiX 1,:) i ) GO TO 123 

I F (HOLD (1,1). <..MW(I,1)> GO TO 123 

IF ( HOL'X 1 , 4) .HI .MW( 1 , 3) ) GO TO 123 

I F - (MW ( I, 3 ) ..AND. MM 1,4) .EQ. .0 GO TO i3D 

HU Li) (1,2) =HOLD( 1 ,2 ) +1 

M W ( 1,4) =0 



c 


c 


c 

c 


c 


c 

c 

c 


c 

c 


c 


c 


c 


c 


i 19 


120 


121 

122 


300 


123 


1 24- 


IF P.O-DUCT 
■ !-’. . M 0 L u ( I, 

5 


)) GU TO 130 


■y\i (!,:)=• 

C t ' t L K 3 ■ j ,r i 

IF (HOLLO - ,2) 

I S =HOL C ( 1 , ) 

I rn= pot ( i ■ , i , 

iKn=nc+,. 

DU I 1 J I 6 = . , J 

DIT KM I .i ' ij;; 1 WORK STA] iOH I ft StQUzNC-. 
W,= P3T( IP, 3 o » * } 

IF { MO. .(). J.-3U) On T fl .20 
cnrjriNU. 

A .-ACHING Tillo HII.J 


IN HOLD BLOCK IS KmADY TO BE MOVED . 

\ 1 ^ f \ rr * 


; hi 


0TB $ THAT PROCESSING 0 -i 


THIS STOCK NUFBER 
INVENTORY, 


H^S Bf:£'i CJMPLr f.D AMO THE PRODUCT IS MOVED TO INVENTORY. 

FI Ad ( I 5 , 8 ) = FT mC ( I 5 , 3 ) + F L . AT { HOLD ( 1 ,2) } *FT AB { 1 5 , 2 ) *FT AB ( I 5 , 19 } 
* { Ir ? T )+ FLOAT ( HOLD { I . ?) ) * F T A R f T r. . r )*FTart f T =.. 19) 


FTAB(I5,5)=F 
Z:R05 THE HOLD 
HOLD ( I , M = 

HOLD ( I , ? ) =0 
HOLD! I = 

HOLD ( I , A) = ’ 

GO TO 10 •• 

00 12 j. 17 = ., ,LO 
QUEUES AR_ 

MO AS THE 
IS FOUND. 

IF (C ( Id, i 
IF 

Cf ; ti-il : 
oO TL> 

:*>Mi nU . 


1 / • i l .H i i n ul u \ i f iL } J ^ r I A d i I 

,o{ Ir, f )+ FLOAT ( HOLD ( 1,2) )*FTA6{ I 5 , 2 3 *FTA8 { I z>, 
HLQCK AFT. R THE PRODUCT IS M0V p u TO INVENTORY 


ARCHED UNTIL ONE I S FOUND WITH THE SAME $.N AND ORDER 
I • T H - H CL l > SLOCK OR UNTIL a T EMPTY QUEUE POSITION 


» r ) . 
t- )• 


J.HCLO (I ,^) .AND. Q( 16,17,1). £Q. HOLD! 1,1) ) GO TO 122 
u. ) AC TO 12? 


I F (HOLD! 3. , , ) . . I . ! GO T.l J 

HOtn tL UT'TY f AG FR-VI *JS WORK STATION IS ADDED TO APPROPRIATE 
•QUEUi P',0 'TIO, , ,F THt NL..T WORK STATION'S QUEUE. 

C (IS, I 7 , _ ) =H‘J Lf ( I, . ) 

Q (16, I 7, :.)=!<( Id, 17, 13 + F0LD( 1,2) 

G( I 0,17, 5 ) =t 10 L . II, 'r ) 

1 1- PO HCLD "M.CLX aT TH3 PREVIOUS WORK STATION. 


HOLD (I, l)=<t 
HOLD ( I ,7) = 

H 'IL 0(1, ) ■= 

HOLIK I, '■) =i< 

GO T(J 120 

THIS SECT 10.4 I iRMIiW, S i'Hl NEXT WS IN THE NETWORK, IF THERE IS ONE 

CUi‘|T I ‘iU’i 

I F { HOLD ( 1 ) .. 0.0 ) Sr It ,2S 
I F (HLLC{ I, , ) .. U.O) GO 7 <j .23 
I S=HOLC ( I , ; ) 

I MU=P5T ( I C , I , ) 

I MC= I MC+. 

THE DO LOOP S, ARCH'S FCR THE NEXT WS. 

DO xZh 19=.., JE 
MO=PST( ID, 19, A ) 

IF (MG...O.IMO) GO TO 1.25 
COHTINL . 

SIHC: THL-IU IS GO NEXT WS,ADD THE COMPLETED PRODUCT TO INVENTORY. 



F T A 0 { 1 S , 7 } = F T - 1 ( I ° » ' ) + FLOAT { HOLD ( I,2))*FTA6(IC,2) *FTA3 (15,19) 

FTAF{ I? ,5 ) = FT /■'.( 1 - , 5 ) + FL j AT { HOLD ( 1 , 2) )*FTA3U8, 2 J *FTAB( 15,19) 

3 1 : TO 1..3 

125 nr. .7.6 11-)=;, s i" 

c ** 3U r UxO AR S - Ci-i.-o UP TIL. QN_ IS FOUND WITH THE SAMI $rj AND CRDER 

C ** NO AS TH f UP, I, ' TH.- HCL1- 3 LOCK OR UNTIL AN EMPTY QUEUE POSITION 

c ** I S FHUNP. 

ip (0 (19, i,?u , 5 . L.0 .HOLC( i , 4) . *NQ.Q( 19,120,1 ) .EQ.HGLfM I ,1 ) ) GO TO 

IF U(I9, I* 1 , _)._0. ) GO TO x27 
x 26 CONTINUE 
9 j TC 129 
i.27 CONTI NUL 

C HOW THE PRODUCT F RON TH. HOLD BLOCK TO THE APPROPRIATE QUEUE. IF 

C THl ORD.-R HO i’;\ THE HU LO BLOCK IS ALREADY ZERO, THE PRODUCT IN THE 

C HOLD BLOCK HAS rj R„ VIOUSLY BEEN MOVED . 

IF { HOLD ( 1 , 4 ) . LQ . * ) GO T u xZ3 

C PRODUCT HAS HOT DEEM ?' CV 0 . TH I S IS ACCOMPLISHED BY THE NEXT THREE 

C STATEMENTS . 

0(19, I2l-,x)=H0Lr'(I,!) 

0(19,12 ,25 = 0(19,12 1 , 2 5 + HOLD (1,2) 

Q< 19,12 , J )=Hi.l!?(I 

128 CONTINUE 

C HOLD BLOCK IS L KQ-D. 

HOLTU ,. ) = 

H L L J ( I , 5 = 
hOL J { I » ) = 

hulu ( 1 1 ' ) =;• 

C - D . <'- l " o WH ,,-« IH.R; IS a PRODUCT AT THE WORK STATION. 

IF (Hi. ( i, „ ) . . .l ‘P.'.M --Q. j) GO TO 13 w 

C N -aT SC VE , ; 3 1 ‘ - “ir < T S NO J 2 THE PRODUCT FROM WS TO THE HOLD BLOCK 

C Apr Z f PL T . i 1 ' ■" -■ , r,0 P? I A T _ UrVS CONCERNING THE WORK STATION. 

H UL 1 ' ( I , . )=*.«(!, . ) 

H PL j ( I , ') = M CL . 

I =i* 1 W ( i , . 5 

-iOUH I , 0 =PST ( 1 o I,*) 

H u L U ( I t / )=:-W( it ) 

m w ( i , n = 

M W ( I , .' ) =U 

C 1 HIS CH C :< 3 V i 3 , E IF TH HOLD BLOCK IS READY TO MOVE. 

GO TO 1: H 

C ** IF o QU> U . IS u JT I VAI LAcLE AT THE NEXT WS,ONE TIME INCREMENT IS 

C ** A Drill: TO 111" r IPISH, 'I I M.. UF ws that is associated with the hold 

C ** BLOCK (FULL HOLD BLOCK). 

129 MW (I , ) = •;«( I, , )4TI'.tC 

KKKS= . 

g n f i j. 

C THIS SECTION S' ,RCIT : S TH QtJiUt POSITIONS FOR SOMETHING TO WORK ON. 

130 C OUT INI" 

DO :56 J= ,LC 
I 22=Q ( I,J, , ) 

C I F THE FIRST QU UE POSITION IS EMPTY, GO TO x39, WHICH MOVES ORDERS 

C UP HNF PUS IT I s, . IN TH= QUEUE. 

3 F ( J.EQ., .AN:j. 122.FG.C) GO TO 139 
C CII.CKING FOR ZJIU STOCK NUMBtR. 



o ~ o o o o o o 


II { i . . ' « * . ) 0 I * J . * • J. 

c IF -H "I 'i. : ° : 1 S i "" I r J* i IS EMPTY, INDEX ORD ;RS IN THE QUEUE. 

IF ( ! - ( i ? J ? ) . L • ) 00 TO i js 9 
IF ( u ( I , J , ? ) .Li . ) GO 70 j.3 9 
:■ • »‘ = .iJv 

131 Oi' ! ( < ) =;’ M R ( I "nKf J 
I ,.4 = 0 

ru < i/ iz j< 

:• i.dz - )= i 
0>5 1 1< ! C ( I Z4) =!'■ 

C -."4 THIS SiXTirj?! IS CQ.mC-RNlL; with raw material availability. 

OH 1 / 4 = I 14 + 1 

c ** RAW HAT-RI^L REQUIREMENT 3 ARE DETERMINED. 

IF ! 104.GT.J4G} GO TO 12 
1 F (0/( IS--) .0 .1) GO TC _3L 
SK (T.- i )=k!W( I ..2 f 1 2 A » i ) 

I OIL ( l:?< ) =RHR {!'./, 12- , 2) 

S>'i=S*( I 04 ) 

CHECKS INV' NT CRY TO 1)2 TERMING IF A SUFFICIENT QUANTITY OF ALL RAW 
M t T 2 RIALS IS AVAILABLE .IF A SUFFICIENT QUANTITY DOES NOT EXIST, IT 
GOES TC STTTuMEOT .36 TO CHECK THE PRODUCT IN THE NEXT QUEUE 

posit i c 

IF (FTpMS !,3).LT. ( FLO AT ( IDI C ( 1 24 )) *FTAB { 122, 2 ) ) ) GO TO 136 

C - US : S R-TURiT TO CHECK THE NEXT COMPONENT. 

GO x, xr: 

\ IN 1 ' THIS Pul N i 0 .NOTES THAT SUFFICIENT QUANTITIES OF ALL RAW 
V-*! RiHS HI ST. IT IS NOw POSSIBLE TO PRODUCE THIS PRODUCT. ONE LOT 
O JO-s r < TH NS \ iG PRODUCT I UN STARTS, 
o ' . i i L 
C 0 •" . U ’ T 1 -T H* j 


SHI = - 1 ( I 



j F ( ) ( 1 , J j 

) 1 'T ( T § . ) ) 

SUT = « 

Tii/t uCH 

: th in ip c ;< u 

c SuTUPS. 

tt- ( c ( : , j , 

1 ‘ ( ‘ m % 

) . . . ;iw ( i , . i) 
~ (. 

SUTCT ( I ) = SLJTCT( I )+i 

IF ( ■ 1 f ! * ‘ • 

• • 5 

) .U . i ) GU TL 

.34 


1 .= w . ; < J 

C ** UH« T riV >T 'AY (STOCK ON HAND AND ISSUES). 

FT A - ( ,5)=F T SN,5)-FLJAT( I DIC(I 25) ) *FTA3 ( 122,2) 

FT, A "L,7 ) =FT V ( SM ,7 ) + FLOAT ( I DIG (1 2 5 ) )*FT AB { I 22 ,2 ) 

C ** D 1 R.'M.. COST '’F RAW MATERIALS USED. 

R ' ii. {!>..) = r 'MC . ( i 2 2 ) + F L C A "I ( 1 0 1 C ( I 25 ) ) * FTAB ( 122 , 2 ) *FTA8 ( SN, 6 ) 

IF ( r~ iX ”' ( 0 i , . Q.2. ) Gu TO 134 

C’tur, L+Fl IaT ( I u IC (125 ) ) *"FTA8 ( 122,2) *FTAB( SN, 6) 

1 :-6 CiVlI ;U r 

C TO HU/jIJCT IS 'TjV ::D FRIO THi QUEUE TO TH- WORK STATION. 

M W ( I , . )=I.V 
m w( i , .)=•:)( J , J ,:•. ) 

C THIS CALCULAT. S TH. NEXT FINISH TIME. 

M W { I , ) =T+ SlJT+P ST ( 122, I, . )*IFIX(FTAB( 122,2) ) 

Q ( I , J , 2 ) = i ( I r J , t ) ~ , 

C THIS UHTcR. 1 1 j 4 * • S IF THC FIRST QUEUE POSITION IS EMPTY. 

IF ( o ( I , J , ') . L. . U. ADD. J. „Q. 1 ) GO TO 724 

3 <J T 'J 7 0S 
7:4 rx.TI'iL 



O' ■ ■ J '.*■)= J L<) . 

C ** 1 -.1* 5 ! J Mi’s - FT R FIRM QUEUE HAS BEEN FOUND EnPTY. 

1 + 

' ( i ; I ' 6 J . ) = 0 { i , i I f J, > 

0(1,1 )-')(], K.7, 3) 

x, 4 U ( 1 ? i . j ” ) “ ') ( *. s i. 1 7 , " ) 

Q ( 1 » f . ) — 

(.(!,. Dt 
..(!». » ‘ ) = 

7.5 C J r < ri.il' 

L THIS LUAD.> Uhl LUT QHTC THE WORK STATION. 

M W( l , ' ) = . 

C RHI.-j •. WHAT SHIFT T Hi PLANT IS WORKING. 

IF (T.L (J*US) ) L =3 

IF (T.LT.C*US) ,AND.T.GT.{2*US) ) L=4 

IF {1 .Li.. ( ' J *U5) ) L = 1: 

IF ( r.LT. ( ?*UC) .AND. 7. GT.US ) L=4 
IF { T . L .US) L= < 

I F (T.GT.nSd )) L= f 

c « D- 1 ' THc MAN AND MACHINE COSTS. 

OHP 2 : ) “ 0 H R P ( 122) + FLCAT ( MW( I, 2)-T)*OHR( I,L ) 

’HiR ,( I 22 ) = It 1R -‘( 1 2 2 ) + FL CaT (MW(I»2)-T ) *GHR (1,5) 

C STM -HEHTo FR'jK HERE TC „39 ARt: USED TO PLACE THE PRODUCT IN THc 

C FIRST CUcU- PCS ill UN IF IT WAS NOT ALREADY THERE. 

IF ( W ( I , J , . ) • H L * ) Gu TO X .1 5 

0 { 1 , J * 1 ) = U 
L'(I,J,;-) = 

135 IF ( J . ET • : ) GL TU .44 
T" ■!ld = C( I , J. * ) 

T. MP?=C( 1 , J , ) 

t, r*p: =c ( i t j » * ) 

GO TC 137 
*2 6 C uf\ T I \ L 
Git IQ IF. 

.37 CONTI <L 

DO ■> ■'* L= ,J 

IF ( L . E i < » J ) 3 .j ! D S 

1 HD =J-L 

I Fk = J-L+ . 

0(1 , I FK , . )=U( )y I HD , 1 ) 

Q( I, I F K , 2 ) = 0 ( i? 1 'TD * 2 ) 

0 ( I , I FK,3 ) = 0( I, i HD , 3 ) 

138 CONTINUE 


0 { I 

f i 1 

. 1) = T 

IMP 

* 


0(1 

f - 1 

.2)=T 

F t-1 P 

/' 


on 

i JL ? 

r 3) = T 

i V.P 

3 


I F 

ft! 

[1,1, 

). 

GT 

. .) GO 

D 0 


I = 

’• » L 

q:. 


1 *. f 

= i: 

>6+1 




« ( I 

1 1 * 

6, . S 

= U 

i» 

i — T , . ) 

: ( i 

y 1 . 

. E , ) 

= u( 

I* 

1-7,3) 

4 v ( I 

1 I. 

. 6 , . ) 

=0 ( 

i , 

Il7,E) 


)=•. 

'■-(I, ' i ) ~ ' 

( I , . i ! = 



c 3 F «'h. FIRST DU •; U : POSITION HAPPENS TO BE EMPTY aLL PRODUCTS ARE 

C SHIFT -0 FOIWAM tjfl c POSIT I OH. 

i *’J 9 L i i ‘ ! t I 1 U . 

1 ( Ct I ,2t:. ) .L-; . 3 GO TU 14 i 
Di . 12 = I » L Li 

c Till. 1 1 S *,'(!■; Ijr s AFTER THE FIRST QUEUE POSITION HAS BEEN FOUND 

.. -'P ' Y. 

I 7= I, -6+] 

•;{ I, I 26, a. )=U< L» 127,1) 

'111, 176, j) = 0(i,I.:7,G) 
a A-0 U { i , 1 J f > 1 7 3 = Q ( 1 » I c i » 2.) 
u ( I » 2 ’ » * ) = 

Q U , 2- » •> ) = 

Q ( I » C 0 » v)=0 

3 0 t 0 i 

C ARRIVAL AT THIS POINT C? iOTcS THAT THERE IS NOTHING FOR THE 

C SU3J.CT rf.PK ST ATI ON TC wORK ON. IT IS THEREFORE NECESSARY TO 

C D: T . R ,J I WHETHER THc WORK STATION SHOULD GO IDLE OR SHUT DOWN. 

C THIS IS THE FUNCTION OF THE STATEMENTS HtRE THROUGH STATEMENT 143. 

141 IF ( T.GT. ir,(I). SulD.T.GT. 10T( I ) ) GO TO 143 
IF (KKKS.EQ.: ) GO TO 142 
IF (JR.IQ.j.) GO TO 144 

C THIS CALC. HAT'S IDLE TIME AND IDLE COSTS. 

l l -2 Co 

IF (T.Gl.TMIl.'jR. ( T. GE. V S ( I ) . AND. T .GE . IUT (I ) ) ) GO TO 143 
IF { T . L . {;.*US) ) L =? 

IF ( T.LT. ( -■❖US) .AND.T. GT. I2*US) ) L=4 

IF (T.L .1 - 1 U G ! ) L = 7 

IF (T.LT. { *US) -AND. T.GT. US) L=4 

IP (T.L . IS) L= : 

lr (T. C; . IS ( 1 ) ) L = a 

i dl ( i ) = i jl-( i ) + : 

J 'i( i )= 1 IK I ) +0 HR ( 1,6) 

DHJ( I )=.J 1 J(I ) +QH!’.( I,L ) 

K 

sc r* , -v 
i-rS T.;(I)=T 
1 *t4 ^ GH (A lU - 
K K K < = 

C THIS u-rTcRMIN : S WHETHER I 0 PRINT OUT INTERMEDIATE RESULTS. 

1 F ( ISToj.lt. , ) GO TO V-o 
I F ( 1 !'A Y. L ' . J ) GO TO , :: 23- 
I r ( JR 5.0 . . I DAY ) GO TO _ f i> 

3333 C.jUTI'il,; 

3 1 ■ 3 I 1. i‘. i T = iS TRT , I E N C, i BY 

IF I'-.-l -. i’*.UG D GO TO 44- 4 
i tqum= t 

c THIS PRINTS I NT SR^ EDI ate RESULTS OF PRODUCTION STATUS. 

CALL WRIT. fa (RPC., AMg , APPLY , FT ABE , FTA B, P ST» RMR , OHR , Q , HOLD t MW» OHRP » 
3FR' , ALT, COST, ORD £R , A, B, C, D, t, F, G, H , U» AK, SER , A3E , I Sc , J5 S , J4.. , 

2W , NCO , L w, JT, SUTCT, TA, I CT, OH , I DL E , DM I , OH J , JM J , I DI C , SK , I NTCT , MS , ISP , 

31 TD'JM) 



c Mil. CK'CK’j T L : L T IF A »-;ORK DAY HAS B;=N COMPLETED. 

±46 I 1 . » L . T") GU TO j. i 7 

KK „ 

C TH's LOOP C ALCiJl AT -IS fPlD N"W FINISH TIMES FOR ANY PRODUCTS 

C IGI'iS *7 \\Y WORK STaTIUU, 

U U . ’ 1 = . t J v 1 

IF (.10 ( it'. ) . ,-G. ) GO TC .47 

"N( ) = ’!')( I , : ) -T A ( I ) 

g rn m 1 , 

1*1 Cf'NTI.IU' 

C IK onfhiGG h • OAImS AT A WORK STATION, THE FINISH TIMS IS SET TO - 

C fn GAOS' A SLmRCH FUR WORK TO BE INITIATED AT THE BEGINNING OF TH 

L IriA T WORK DAY. 

MW( 1 , A)=- 

:-.*c Conti ■hi 

L iiT=NT + .1 

DO • o- L = L.gT, i5(J 

C ** UPDATE- I NV: illoRY BY SUBTRACTING BACK-ORDERS. 

F 7 A L ( L , S ) = FT A R ( L » 5 ) -FT A3 ( L » 17 ) 

I F ( FTAB( L, 5 ) .LT.O . ) F7A B{ L , 17 ) =-FTAB ( L , 5 ) 

IF ( f T A B ( L , 5 ) .G-.O. ) F TA 6 ( L , 1 7) =0. 

It- { FTAB ( L , 5 ) . LT . » . ) FTAB( L, 5 ) = ■"■ 

F TAB ( L , ; *J ) =F T -aB (L f a S J + FT aB ( L » 17 J *FTAB CL, A4)* ( i. /FLOATt I SP ) ) 
FlAB(L, ,6 )=FTAB(L, . 6 ) + FT uB ( L , H ) *FT AB ( L , 5 ) * ( 1 . /FLOAT (ISP) ) 

307 C OUT I ”!U - 

Do . „ 1 = , NT 

C ** JH.-T IN 7 NiT'jRY RY SUBTRACTING 1/6 OF WEEKLY DEMAND. 

DM FT 'HU I , 3 ) =F7 AO { 1 , 5 ) -FT AB { I » 17 ) 

= • AM I , j) =b~i M 1 ,5 )-( 1 ./FLOAT ( ISP ) >*W< I) 

C ** CjMNT; CARRY I !G COSTS FDR TH" WEEK. 

FT- A I,7)=FT'-?(i ,T) + (a ./FLOAT (ISP 5 )*W(I ) 

It i != ,r S(F" M { i , 6)- .5) 

If (KMM I ,5) .LT.O. ) WRITE (5, IBS) I , IBH, JR5 

It- ( F T A B ( 1 , 5 ) . L T . 1 • ) FTaMI, j.7 ) =-FTAB ( 1 , 5 ) 

IF ( FT A- i 1, » . ) FTa:(I,a7)=0. 

IF [ FI AM 1 ,5 ) .LT.O . ) FTA? ( I , 5 ) =0 

r i f :;( 1 , C ) = F T - 1 - ( I , . C ) + FT . B ( 1 , 17 )*FTAB ( I » 14) * ( 1. /FLO«T( I SP> ) 

r I A 3 { i , . S ) = F T . G ( I » „ 6 ) + FT B ( I , 3 )*FTAB (1,5 ) * ( 1. / FLOATl I SP ) ) 

3 1 L COM 1 • ' U - 

C ** CM -Ml ORD.R Tor. A Y 10 t) £T -RHIf-lC IF PURCHASED PARTS ORDERS HAVE 

C A R K I V . 0 D'F.INj 1H‘ DAY -UPDATE IMVEuTORY AND ZERO ORuER ( K» J ) . 

J T P - 1 ' R D R ( K , G ) 

IF ( JT°. E ;. ) ) GC TO 432 
1 .'.L = 'L T ( K ) 

ALTD=ALT! <) -FLOAT! I AL ) 
a J TP = FL’J A f ( JT?)+ALFD 

AKT t> = FLiJAT(KT°) +FLUAT ( JR ) / FLOAT ( ISP) 

IF ( A J TP . G T . AKT P ) GO TC - N.*> 

SI ZaiJRD.R ( K, T } 

K K = i R 0 R ( K, 1 ) 

FTArt(K»(,3 )=FTiB(KK,S) + SIZ 
c TmB(KK,5)=FT -- iJ ( KK , 5 ) + SIZ 
FT' -;lKKu 3 ) = F^A M K K , 15 )-GlZ 
,1". • MO) K , KK , 3 I Z , J (5 



UCJil -FTAiKKK, ) 

IF ( FI ii;i( KK , <V ! .Li. . ) GO TO 1513 
IF (GIZ.G, .FT ,«•. (KK,0) ) UCOST=FTAR(KK, 11) 

13.3 CoST ( KK ) = S IZ *UCOST + CCJS T{ KK) +FTAB ( KK ,4) 

HRT = HRT+S IZ*UCf!ST 
OKI). R(K,1 )= . 

3R[J;R(K,2 ) =0 * 

CiRUwRlK, J )= . 

*+ 3 CC JTIwl'f 
2 CHNTINU. 

Src = o 

C THIS CETERMlH. S th> shift CHANGE COST. 

DO ’.49 I = _ , J 3 

IF (OHR( I , : ) .IT . . . ) GO TL 149 
IF (MS(I I ) } SSC = SSC + SC 

NMS( I )=M$ { I ) 
x49 C UH TlHtS 
I TDUM=Tt- 

C THIS WRIT .5 TH_ STATUS OF THE PRODUCTION SYSTEM AT THE END UF 

C THE WFEK. 

CALL WRIT ; C ( K M C„ » , AP^LY , FTA«S , FT* B, PST , RMR , GHR , Q , HOLD, MW , CHRP , 

1 OHRM ,A LT s CuST, ORt’ ER , A, B , C, D, E , F , G ,H ,U , AK , SER , AB E, I 5 0 , J 50 , J4G , 

2M , NOD » LO , U T » SU T CT, TA , I C7 , OM, I DL E, OM I , OHJ , DM J , I DI C , $K , INTCT, MS, ISP, 
3ITDUM) 

WRITF (6, _32) 

C THIS WRITES THE IDLE T PL:S AND COSTS FOR EACH WORK STATION. 

OlJ x 55 I = i , JSO 

I F ( I JL'_ ( I S .EL. ) GG TC _55 

WRITi (5, . *1 ) I , 1 DLr { I ) * OHJ ( I 5 »OMI ( I ) 

155 CONTINU. 

C THIS uCTSiV'T!. S THE TOTAL IDLE TIME AND COST OVER ALL WORK 

C STATION. 

I GU ■!= . 

DO » — > o — . ? J ^ . 

I Sl'0 = I S'OM + I DL .( I ) 

0 = 0+0 11(1) 

A I DL= 4 IPL + ) Jli! 

156 CUuTIUU: 

I DL= I SUM 
BC=A lPL + 0 

WRIT’; (5, 54) 1 SUM, M D L , „ , HC 
WRITE (6,..; 5) 

IT j A — 0 * 

C THIS WRITES TH-. S:T UP C_3T AND SHUT DOWN TIME FuR EaCH WORK 

C STATION. 

DL x*,i6 1 = 1, J 50 

WRITE. (6, . S 7) I , SUTCT ( I ) , T A { I ) 

I TTA= I TTA+T A ( I ) 

186 CONTI ML; 

R t TURN 
C 

158 FORMAT (315) 

159 FORMAT (7 IV.) 

1 60 format (3F‘ o. •, 1 ) 

1 61 =nR'IAT (4h . ) 



o o o o o 


162 Fdki'UT (51 0) 

163 FORMAT ( X, ^PRODUCTION ORDERS, TIME AVAILABLE AMD £RR0RS*/1X, 

3 4 X , ( » H * ) / / / ) 

164 FORMAT (21,. ) 

165 p ; t,< AT { Ji“, ;• , } *3T0CK NUMBER*, 15, EX,* HAS NO FIRST WORK STATION ... 

I ORDvR FOX *, i: ,*BATCHES*) 

166 FORMAT (, H »5X ,* STOCK ,-l LH .. ER* , 1 5 , *DOES NOT EXIST ... ORDER FOR *, 

II : , *BATCH .S* ) 

;t,7 FORMAT (5 a,*I'(I"1AL QMEU. FULL FOR STOCK NUMBER *,I5,* ... ORDER 
„ F U R * , I ’ , * K 4 T C 1 1 . S* ) 

168 FORMAT ( I - ) 

169 FORA. AT (* *'. * 3?1 - . i ) 

170 FORMAT ( l HO , / , - X, 1 b ( 1H * ) , *TI hi AVAILABLE AT WORK STATION *, I3»2X» 
.^GREATER THAN THR-'c SHIFTS - THREF SHIFTS ASSUMED* ) 

171 FORMAT {j.11 ,/ }..X,a ( 1H* ) , *T I ME AVAILABLE AT WORK STATION *,I3,2X» 

1 * L c $ S THAU Of; SHIFT - 0 «E SHIFT ASSUMED* ) 

172 FORMAT (5X »*WUBK ST ATI CM*, 2X , 14, * DOES NOT EXIST*) 

x 73 FUR. i A I { * H ,/, x i,l (- H*),*TI Me AVAILABLE AT WORK STATION *,I3,2X» 


;*Nor sivEu - s. 


SHIFT ASSUMED* 5 


174 Fu‘V aT (110,415) 
j.75 FORM '* “ (31.0) 

176 FORMAT (iH »i?v.»*ALL OPD..RS PROCESSED*//) 

132 FORMAT /. ,*I DLF 7I'E*,/,61X,9(iH*)///33X ,* STATION*, AuX , * TIM 


. r- * , ; OX , *M / ’ !!* , 


♦ MACH IN _*, / , S9X, *NUMBLR*, 1 IX , *1 DL E* , 9 A ,*CQST*» 10X 


2,; X,*CCST*/) 

133 FORMAT { 3 8X» 1 6 , ilX » I 5, 4X, Fi. . 2, 5X » F 1 j .2 5 

*34 FORMAT ( IH ? 27X,*T0TALS*,16X,I1 <»5X,F9.2,6X,F9.2,5X,Fxw.2) 

165 FORMAT (*H. ,56a,*WDRK STATION DATA */57X , 17 ( IH* ) //39X, *WORK STATI 
j 3 N* , .. *, X , * o r T U PS * , 1 o. X , * SHUT - DO WN** / / ) 

1B7 FORMAT ( 40 v , 1 ' » 13X , 14, 2i a, I 4 ) 

i83 FOR. IT ( >. v,*f R JuUCT NC* « 15 , 2X , *OUT OF STOCK *, 17,* UNITS ON DAY * 


500 F I..-P („ !' , * jR'iuR Vj *,15,* FDR SM *,I5»* FOR *,Fi .'.,* UNITS ARR 

llVtO UN D^ v *,i'/) 

5 55 FORMAT (31 ) 

L 


SUBRGUTJ N 
10HRP jfjHRU, 
2W » H G J , L 0 » JT 
3ITDIJM ) 


1/F‘\ IT . O ( P.MC1 ,!'! ■ S, A FPL Y , r I “ c O » r I AO, rol , rvirs, un^nnuuuu, , 
'LT,C'?ST , ORDER , A, B, C, D,L, F, G,H,U,AK,SER,ABc, 150, J5 j , J4-. , 
, SUTCI , TA , I CT , QM , I ULE, OMI , OHJ ,O l 1J,IDlC,SK,INTCT,MS,ISP, 


** THIS IS SUBROUTINE WRITE-.. **-************************************* 
♦♦ PURPOSE - 1 0 WRITE INTERMEDIATE AMD FINAL STATUS OF THo PRODUCTI (DU 

SYS i_v. 


COMMON 
1? 1UVTY 


/CLCJk/ -,ITP, T X,TIUC,HRT,0C,5PMLT,US,CMU,CMC,KTP,IP,V0H, 

SC ,i'H, i . L, C Jfi , SLCjTDI , TDM, I FT P, I DL , ITT A, VI NV P, N ,DUM 



ni •' I'oin.'j ~MC.( 19 ) , ^1S( 3) , APPLY { 29) , FT ABB ( 29 ) ,FT4B{ 29,19) , 
JPSlt ?9, ,4 ) , RMR ( ?. 9 , 5,p) ,OHR( 6,6) ,Q( 8 ,20, 3 ) , HOLD ( 8,4), 

2i^( w,4),«':HR'M 29 ) ,DHRM ( 29 ) , ALT { 200 ) ,COSTC 29 ) , ORDER ( 200, 3 ) , 

3A{ ; ) , 8 ( :• ) , C { . ),0( >),z( 3) ,F ( 3 ) , G ( 3) ,H( 3),U( 3) ,AK( 3), 

R( i ) , m d ( 2 ) ,W( 3) , SIJTCT ( 8 > ,TA< 3),I0T( 3),0M{ 5), 

6IPL.( T ) ,U II ( 0) ,UHJ { 3) ,0MJ( 3), ID IC ( 5) ,SK{ 5>,MS( 3), 

61 i .T CT ( :, ) 

I NT , G -R T A » U S , P vi , U T » TM ,P. !U, HU » T , HOLD , SN,Q, SUT, TEMPI , TEMP2 , DM, SK , 
IT MP3 , SUTCT , T I ilC, ORDER 
I NT lG 1:R P ‘ j T , R m R 

COMMON /PR 3JUu/ PN , OT , TM ,PNJ, HQ, T, SN , SUT, Tl; MP 1 , T i_MP2 , TOC , I SUM, DB 
WRIT.: (6, .9 ) 

WRI rt (6, .77 ) ITOUi'i 
WRITc (6, j76) 

DU 152 1= ; , J 5 ) 

0 0 ;51 J= ? , LQ 

IF { Q ( I * J , . Q. ) GO TO 15 a 

1 A = 0 ( I,J, 1 ) 

I ii = Q ( I , J, 2) 41 F I X ( Fl'ABt IA ,2) ) 

I C = Q ( I , J, j) 

WRITE (6, .79) I , IA, IC, IC , J 

151 C GMT I : ) U " 

152 CONTINUE 

WRIT. (6, ) 

WRIT-" (/,, .7 h ) 

D J .5- I = , J 5 0 
IF (rw(I,-).L ) GO TC .53 
I -j — M W (1,1) 

I 0= I F I X ( FT '3(12,2) )*MW(i ,4) 

IF ! IG.Ll . J ) 3 u TO 153 

v. P, i " ( o , _ ' 9 } 1,MW( 1,1 ),'1W( 1,3), 10 

153' COUTI'lL ' 

WRI It. (-„, . ' . ) 

W R IT z ( A , _ t ! 

DO , 5f I = . » JS 

IF UrJL"( a , ) G1 T : _54 
I +=HCLC ( I - , ! 

I 6= FT A 2 ( J \ , C ) *FLOA f { HOLD ( I , £ ) ) 

IF { Io.L . ) ;n T] '54 

WRIT,. (u, .??) I ,HOLO ( I ,1 ) , HOLD ( I , 4 ) , I 6 
154 CONTINUE 

131 FORMAT ( ) v,*riuLC BLOCK *// ) 

179 FORMAT (Ln.',K, 4/'. , 1 5 » 17 , 1 5 , 16 X , 1 5 , 1 5X , 15) 
a77 FhRM',T (5 v, WAITING LIN_S AT TIMS = 4,15//) 

18. FOR' AT <_!!•), DuX, *19 PRCC.SS 4//) 

178 F0R-1AT (2 ?X,4aOT<K *,*UX,*STOCK NUMBER*, IOX, *QRDER NUMBLR4,1JX, 

1 *OUaUT I TY4 , 1 ^ a, 4RDSI TI CN4/£3X,*STATIQN4/) 

190 F GROAT (1H ,,5:x,*ST*.TUS OF PRODUCTION SYSTEM 
RETURN 
END 


4/53X , 27 ( 1H4 ) /// ) 



o o o o oono o a o o o 


S 1 MRUUT I f4 1_ I N VLuT ( RMC v ,*MS t A PPLY , FT AB 8 , F TAB , PST , RMR , GHR, Q, HOLD , MW, 
.DliRPtOHRM , ALT, COST, ORDER, A, 3, C, D, S, F, G,H,U* AK, S£R, ABc, I5v, J5 0, J4v, 
/W, l4Uf j T LC , i T , 1 SP » SUHRY f IRuW ) 

C 

C 444 THIS IS Stl.-'KOIJTIHP INVE>s' r *44 4# £#£*4# 44444444444444444444444444444 
C 444 PURPCSE- I 3 aCCHMULA i' 2 ALL INFORMATION REGARDING INVENTORY AND ALL 
C CONC.R-.icD COSTS. 

C 

C 4*4 DF.SCR I PT I - 1: ! I IF VAR I AbL £S- 

C APPLY ( L ) - ACCUMULATED COSTS THAT APPLY TO STOCK NO. L. 

C FTABO(L) - A RUNNING TOTAL OF THE QUANTITY OF STOCK NO. LTHAT 

HAS 3ttN PRODUCED. 


444 ROUTINES US*- D- 
4 4 p ^AP 44 


■ /CLCJR/ H, IIP, I x, fINC,HRT, OC, SPfILT,US,CMU,CMC,KTP, IP, VOH, 
x!-: uv 'Y,SC ,0H, aCL, CUM, SSL, TDI »TDM , I FTP, IDL, I TTA, VI NVP ,N» DUM 

"’I .S Z-j RMCx( c9) ,M»S( 3), APPLY! 29),FTA38( 29) ,FTAB( 29,19), 

IP ST ( PS, ’3,4),nnR( 29, 5,3) ,OHR( 5,6) ,Q( 8 , 20 , 3 ) , HOLD ( 8,4) , 

2MW ( 3,4) ,OHRP{ 29 ) , OH RM ( 29) , ALT ( 200 ), COST { 29 ) , ORDER ( 200, 3 ) , 

* A ( ; ) , B { ., ) , C { •>), D( 2),„( 3) ,F{ 2) ,G( 3) ,H( 3) ,U( 3) ,AK{ 3), 

‘+S _ R ( 3 ) , A Rf ( 1 ) ,W( 3) , SUMRY ( 5 3,15), FCST (3 ) 

I IT-GlR T 'i , US ,PN ,QT, TM,PMU ,HQ,T , HOLD, SN, 3 ,SUT,TcMP1*TEMP2 

I I T G 2 R PST, RMR 

I ••t.GER 0 R0_ R , T I NC 

r«st= . 

0?C.lS7 = . 

! ; C S T = . / . 


THIS SU 'Ll ?• 5 , 1 Z S LENT A I i\ .„JSTS, LABOUR AND EQUIPMcNT COSTS 
( OTH 1 R THA. ILL -) AND TUh CARRYING COSTS. 


on _ _ 2 L= _ , I 5 1 
3 SC =C SC+F ~ ' A ( L, ... ) 

0 PCOST = UP COST +0' iRP { L ) 
t WCUST= ’ : 0OiiS T + ! MlKM { L ) 

:CGST=CCUST + FT;,rt(L,26) 

TH= FOLLOWING ST AT CM; NTS ARC USED TO CALCULATE AaiD PRINT UUT THc 
VARIOUS WvLKLY m!,Q CUMLL-TIVE COSTS. 

A PPLY ( L ) = JHRP ( L ) +0) iR’^ { L ) +RPC i. ( L 5 +C0ST ( L ) +APPL Y { L ) + FTA3 ( L , 16 ) 
OHRP(L)=; . 

OHRP ( L ) =0 • 

3HRM( L ) =0 . 

3HRM ( L ) =0 . 

R f C • ( L ) = . 

COST! L)= . 

IF ( APPLY (•_). ■' ,.j. ) GO T-j li2 



A. X tl w L * I 1 i J L ' 

W A 1 1 _ ( j i .. '5 

WUifl, ( A , ,) 

D'i ... .1 !1= - 

J F { FTAiiC 1,1). . ) GC m 1 13 

j J . KK = , : i 
^00 FTVHKKj Jl ‘)) = r . .. (KK t ..r) 

WHI1. {. , . ) I , KU b { ' VU » r T A -3 { M T 7 ) , FTAB ( M, S > » F TA a < M , i S ) , FT AP. ( h , a 

i) , f-TAB F ) , ? 

FTAd< ! , - ) = . 

F f / ,B { i' 1 1 .6 )={»..■ 

171 FORMAT {F .. ), r ,r* . ) 

113 COWT I It' 

DU UK- =. ?i ■ 

FT AC( I JK 5 » ' 5 = 


F T A 3 { IJK 5 , * 

I = 


wri r; { At 

•T ) 


WR1TF {*., . 

r ) 

-COST 

WAIT _ ( 6 ,;.-' 

,) 

. .COST 

WRIT,' (o, . 

7 ) 

cm 

TCST=PPC J3~ 

+ ’ 

,JC .3 T +C«*i , j 

WRIT:. (£>,.-• 

6 ) 

T 3 T 

3; C=Ci ' C + T CC 

T* 

i 


WRIT:. {A». • 

T) 

L 'iO 

WRITF ( 6 » * - 

. ) 


T "C = L','C n ST+ 

< 1 


T i'<C-OPC: jS T + 

,4 I. 

_ ! 

WRIT-. (A,.. - 

A ) 

TOC 

WRIT 3 -■ 

} 

1 :C 

W *\ 1 1 r ( F j . 

7 } 

C”U 

H 1 = . I i + V f H 

• iC 


W IT:. ( 3 1 . 

) 

A j 

WRIT,- .• 

' ) 

«, J w 

TTO = T ‘C+r >'■ 

+ h 

4 > r ,f > 

Yt 1 i : ( " i . 

■ 5 

* L 

T F ', = Ti.y + T- (. 



UUT, j 

i 5 


wp.r* . (a, 

’ ) 


W \ II. ( - ? 

5 

r 

WRIT- 

. ) 

1 , _ G T 

WRII : U, . 

) 

. ‘ 'U 

T IC = 'JC + LC 

r+ 

. 1 - 

WRIT-. (':» • 

b > 

t { 

TDI=TDI + r 11 



W K IT 3 (-, .. 

7 ) 

■ < j. 


V . i 1 =0 . 

to A 7 ;= 5 1 • 

217 V 01 = V 1 1+F T • ,M 
XC=SS=". 

I F (V'jI.GT.p* 
ALU' N = 

50CO COOT Let- 

IF (P-NaLT.L.. 
A L T ? r ‘ J .<• L >' 


f i ! * i- I t- ^ { i ; "> ) 

r> Y) ,-C S:>VUVP*{Vni-Pr!VTY) 


30 TO 5 


) v ■ I c 



o o o 


I 


AH u < - L^-fl-T+p . -lT+XC -SS 
P . i P . ‘ : 1 = 1 ! “’ L 1 i S G 
W \) { f . .... ) 

TK-T7C + TI 
W : \ j T - ( v ; ) ffC 

Cl ;j-C'JF + ~r>r 
H ;ii ( , "■ , ) c i I’l 

W i' ! ™ . { . « . .• i ) H'-' r 

wuir- ( ) m:a 
S UIR - =CT1+ ' , P ■ L7 
WI’.J V - { ' i 7 ) P , i c 
rtnif ( ",) :p,,l t 

W R 17'- ( -> i . V ) £ CJ'~ j 

TTL'iP = J-n- ■. 

wmi {('t: } ] rnp 

1177 F JRPAT { „ H .JrfbcR * > I 5 ) 


CjIjLL 


lap 7 r 


iW'ARiZL vL.KLY RESULTS, 


CALL R_C AP( CC C t , J 1>C , 7 PC i VGI , SCuR* » FC ST f W, .-IT t SU^RY » I ROW t FTA S , I SP , 
.TIC) 


£Cl= J.< 


1 ,. C , 

Fi'k 

- T ( l : ' 

* 

) 

l.H 

F ... 7 

7 (-=, 

i ; « 

0 ) 

1 /, c 

F .. '• 

r < . - 

, 0 

) ) 

j V 

F ' i>‘, 


J 1 

* J . - ^ f = ) 

* c 0 

F. K 

-7 (. ■ 

* ^ 

i Ul «•' D F L I" l CCST HILu HACK*) 

i U 1 

F 

" 7 ' - 

f ^ 

,mI ,V T'. -Y crATUS*»/» 59Xtl6(iH*)///) 

1 J i 

*7 • i ■'“V t 

( * s 


>, y’,*L I IT C;jST*t L v f *isS!J = S*»-*X,*RCC=IPTS*»9X,*0N 

* 

J i 

IJ j • ' i 

' V A 

"7 Y C jCT * ,7 *STCCK C l H AMD* 1 7X , *Sfi* , / 5 

a I'j 

fup 1 

- T (7 t 

I “ 

, J . , c (77, Fi . 1) 1 7X,f:,: .4,BX,F_ .E» 6X, 14) 

* f U 

F i ,• 

: {J, 

7 3 

: , *C S- -Ja7S*,/,60X, 12(*H*) ,*.6(7 ) , 53X, *DIR1CT 



F C ’Ll i 

j 

C S~3*,//> 

\ % 

r j » h 

I. 

1 fc 

• : s *L*,3 r iR C ST (M.'tT I DLc )*, 4X, Fi .2) 

L 2 6 

p : , o 

.7 ! . !• 

? 

' GUST ( JUT IDLE )*,2X» Flo. C) 

i ;7 

F UK 

( - 

’ f 

, ; * . - T r RI^LS ub=n CcfST* j6X ^ F iQ • 2 j / ) 

j, p 

f;f/ 

{.n 

f 

,?T ,T'L Tr-IS p:f.!0C*» ;X,F1 .2) 

U 9 

P u rl r ' 

.! T { ; ‘ 1 

3 * 

- ; . * T , i T / L I _ J4T_*,l?y »F15. 2) 

i 

FUP " 

- V { < 

f ^ 

,.»*!• jTl L <i >• 1UF Af-TUR I \G C'JSTS*///) 

I u 2 

Ffip 

;7 {..j. 1 

; f 

■ ’ t *7 tT/ L L .BCK CCjST*? 9/ , Flb.L 5 


F 

•- 7 ( . M 

? 

’,*T'i7 il ’"CHI j - C-JST*»7>fFi • A ) 

14* 

P L Rf 

, . T ( , 1 1 

? 

' i * .7 V „■ H !•) 1 1~ . * t 7 7 X ? „ * X ) 

*45 

FuR 1 

-.7 < i 

f *« 

-■ y * C ■ . I F 7 Civ’ 1 l'3'„ C 3 S T * i Li X i F * 0 » i. ) 

146 

FOR 

■/' T ( : M 

S f ' 

- ; , *7 •: 7 ■ L 7 -IS PlRIUI)*i1£V» HO.Z) 

1*7 

FORU 

— T { » H 

? - 

, ,*T;iT. L T ,_ DAI .*05/ »F 15.;.//) 

B 

FORu 

O T ( u 

f b 

; un 'V~ n* cv cs ^ ts*// ) 

i 9 

FOR’*’ 

*, r* * 

< *i-i 

t * - 

. %*v LG- l.T • ; A 7 _ R ILLS PcCLWF.C THIS PlRIDD *, 6X» 

£ 

/- r U J m 
= OR»i 

< » ) 

A T ( . r! 

? '* 

, * i-O. ", C ST* s . 5X ? FI . . 2) 



o o o o o o o o 


j 

L 

\ 

M 

( . 1 > 

! t“ ’ 

■ , C 

1 ** * 

F j 



t „ i 

? f 

f "< c j t ) 

' k, . 

f 

»i; 

6 ! 

< .i 

f 

\ f / ( 


F 

*. 

'* j 

Ui 

f ‘l 

| #T u 

, u. r, 

u + 

f * 


u y 

f • i- 

If * 

; i 6 

i. ' *'i 

f 1 * * * S < 

i n. < 

i t 7 

C U 

/ 

f 

. J -1'^ J 

D7 

F 


' 1 

( ;.u 

t * 

, * p 


* r. 


F ) 




*5 >* M 
i* <* 

Fh 

.k n 

r 

f . i: 

'f A 

*tu 


; F * 

««* * 

7 I 




, s C 

F ' 

*K: ' 

; T 

(. 1; 

,////// , 

*. 63 

F>. 

} f 

i 

(A! 

,D 



Ft^Y c- 3 T *, '5Y,FM. 2) 

'VDliJT 'If- .KICK COST*, LX, F*. .D 


i '<■ L PlA'-IT COST, TO 0ATn*?7X,F_5.Z) 

TP h r TOTAL VAL'R OF INVENTORY * , 1 IX , FI . . 2 ) 
,*p DLTY AML LXTcRiRl COSTS THIS PERIOD *,SX, 

. P MALTY AMO EX TcRUAL COSTS TO DATE *, i X, 

;,*iUT4L SCURF TO DAT t *,-lX,F15.2) 


,c ■ v p . c* ° ( ccost ,usc, tpc, voi , score, fc$t,w,ut, sukry, iruw, 

. c ” . . r- id 

c 

C THIS IS S i j ;; . i )'f I • IF R-r,AP ♦ £*##*$ £#**£*$ **$#+********* ************** 

c *** purpose - rn store fi ml results for each period. 

c 

Si.'i "< / / C L " J K / 0,ITP,IX,:iUC,HRT,nC,3P iLT, US, CMU, CMC,KTP , IP, VOH, 

V fU/1 Y , SC , <H, I DL, CU”1, SSS,TIR,TDM,IFTP,IDL,ITTA,VIMVP,M,DUM 
DO , . I < '■( 3),FCS7(D ,5U">RY15 , D ) , F TAB ( ?9 , 19 ) 

ID ' R 71 'CDS 

*** l * i... , ' • , ;-.. l L rX..ClJ'i IU.J OF THE SIMULATOR, A SET OF ib VALUES 

L .. ■ j : . ' US:. IS -,'U i.kY TuRL;. FUR FRH WEEK OF SIMULATED 

DDRATIOU 

*** THI , ■' i ■ J 7U Vi D mj. ( y r‘W [,0. ) FUR STORING The UEw 

.U ’■ ' RY 0 Tt . 

I = ; ^ - 1 

p p i r { - f ) 

D FliR, :.T ( . i ! , ' . j * S 1 1 j L A T J L i UUUDSR *,10} 

SU'RY( ) = :.•„ 

S ! i ?. Y ( I , ^ ) ~CC JS • 

SIHDY(I,„ ) = j'x 
5 H. kY( IrE D = TD 
SO SRY l I , . ) =>_ 'J.‘ 

SUDYU,, . 1 = Vj I 

S U : ' l Y ( I , „ S = S' J i L T 
SU? RY t ) = 3 r D<~ 

S i >UR Y { I, . ) = FRi(. ) 

SUKRY ( I , '.) =Ui ) 

SUHRYS I, A ) — F T -D ( - 5 
SUMRY ( I 3 -TI ' 

S UM RY ( 1 ,7 )=iT *RD-I .,L 



o o o o 


Si H' UJJ T M 


K 


! ■ I ( i A i R * j \ ) 


\ i , 


THIS IS O-i I I I | i fi ' ‘.M'J ' ^^❖❖^❖^^t:^***^^*^^#*#****: 1 !'**#****^*^#*** 

*** HJRt-'l'JS: - <T; ,j_, . RAT, RA .OOM NOS. 

X I = F L 0 A T ( 1; )*, j .n 1-7 

X. = WHDY •(.;!) 

CALL SNOY .<;',.} 

4 '< - s' Y ( ..) 

i ! "s-.L; • ; 1 « t ~ 2 . "" ) b 1 1 7 b 4- 

> — •>.••”> V I l' | 

I '= . ' ! . / 

R ’ 1 >’• 

': f / D 



********* ** 1 :********ZZZZ#%##%*%%#Z*********** 
PROGRAM FOR C3RR .LATIOiM AMD SPECTRAL ANALYSIS 
*********** !■•**********#***#*****#**%%***%%*$% 


C 

^ ❖#* this poor mm is used i m appendix a. 

C *&* THIS IS PR.iSA P.-.UGRAM 

c PURPOSE - TO V..ST DATA F'.R CORRELATION AMD SPECTRAL ANALYSIS. 

C 

D *** DtSCRIPlIJj OF I DPUT PAR / MSTcRS - 
C J X { L } I ) - p;_T DEMAND HISTuRY. 

- L= VARIkul. MO. 

C 1= SAM?L ';0. 

C NS AMP - UU. OF SAMPLES. 

C F jR AUTOCJRREL, Ti M N VAR -- .„ . 

^ :1I -AG - MAX I HUM LAS NJ’-6_R { BFTwEcN M SAMP/zu TO Ns AMP/ E ) . 


C 

C 

r 

w 

c 

c 

c 


c 

c 


c 

f 

u 




r 

L* 


c 

r 

o 

c 


0 1 So RI P TI j'*l OF VARIABLE - 

TP(I) - SUM OF V«RIAB_E JX TILL PERIOD I. 

FPU ) - SUM OF VAR I Ac L c JX TILL PtRIDD I + IP. 

SP(I) - SUM JF S 0 J - R: OF JX TILL PERIOD I. 

GPU! - SUM' : : F SQUARE OF JX TILL PERIOD I + IP. 

SPCn - SUM ,,F PRjDUC" OF JX TILL PERIOD I AiMJ I+IP. 

RP(I> - S.-islAL C'R.R LL - T I ON COEFFICIENT BETWEEN VARlABLn JXU) 

A • 1 »• ja( i +■ I P ) • 

^ H . 5 i p =p = l A G MO. f . , I, L» . . . , MLAG) . 

XP(I) - aUI JC JV -\R i A NC FJ MOTION OF ORDER I. 

FL U ( i > - *- * ESTIMATE OF OF A TRUc POWER SPcCT.RAL DENSITY 
FurCTIOiM OF HARMONIC NUMbER I. 

D p U) - AUMfjTH LSTIMATt OF A TRUt POWER SPECTRE L DENSITY 
FUNCTION OF H RMOMIC NUMBER I. 

-V .G - ]_A . VALUr. JF AP( I ) . 

UGM-, - S i i-n' 'ARD 0u.Vl-.TI3N UF R D ( I ) . 

FR.C - FR_UU cCY. 

C SRNvJ - MU 13 R OF CECIL'S OF FREEDOM FUR SPECTRAL CALCULATION. 
GFRAR - M :ki VALU_ GF SMOOTH SPECTRAL DENSITY FUNCTION. 


01 Mu US I ON FMT(„ }»/(., 1, ),TP( L . )»SP{. - ) , FP ( _ , ) , GP { } t 

iCPI. ),WPU ),« D (i. JtUPUOO) ,FLP{_00) , TITLE (^0) ,JX(^, EOU) 

Ri. AD 60Ut II TL 
R.AD • , MS AMP, ,;VAR 
1 FORMAT (Al; 3 
PRINT 7)0 
3 HI ■;! 6 JO , !' I T L _ 

PRINT 7 .• 
t> ; FORMAT! , 7 ) 



7CD = j.< '\‘\1 i-AA-ZA H*3 )/} 

I r ( WAR) 

4 ' V M V A b = .. 

2* i 3 JNT ulU : 

IAO , , ( ( J). (L, I ) , I=1,NSA IP) , L=i , NVAR 3 

A z DR, MAT (..?!-') 

3i< . i L=d i V k R 
3 3 1 = 4, "ISA MP 

X(L,i ) =J X ( L , I ) 

. 3 ^ JN F I h U . 

M = * J is A F \ P + i 

!WtG=-i./ { FLIWTUSAMP )-l • ) 

S1SMA = SL'R« (FLJmKMGAMP ) / FLOAT ( NS AMP- a ) 

twosig= .yt*SiSf.A 

3 U R l) P = A V j 3 + T W J 3 i. G 
3 R L 0 W — A $/ c s? “ i WJ-> I o 

>RaMT lUi, , AV..G, S i 3 M4 » T W33 1 G » COR UP, C0RL3W 
13 a FORMAT { / /„ X, *PAR4M = T= RS FOR CORRELOGRAM*//5X,*M£ AN =** FI v.6»* iT 
13. utV. = * , F 1 .6,* TWO SIGMA =*tFi,.6 f * JPPER LIMIT =*,Flw.6, 

2* LOWER LIMIT = *, FID. 6//) 

15 3 3NT1NU,-. 

RlwD 3, ML AG 
3 FORMAT (15) 

I FfMLAG.L -.U) v’.i TU AOO 
3u .00 IX- t NVAR 
3D I Y = „ , NVAR 

F P ( a 3 = • J • 

SM . 3=0. 

=P( ) =0. 

3 D ( . ) = *. 

33 . w I = , , !l 1 5 A M P 
Sp( . 3 =5 :> ( . ) + X ( i » I 3**2 
T P ( . )=TP( 3 +/{ I. ,1 3 

GPl 3 =GP ( . 3 + X ( I Y » I 3**4 
F P { . )=FP( „ ) + a ( I Y , I 3 
15 3 3 jTI MJ 
M = ;• L-G+. 

3 j 4 1 = 4, )_ 

J = I - 

< = MSI.’ I P-1 + : 

TP( 1 3 =TP( J) ~.A 1/ ,J) 
sP( 1 ) = SP( i-x )->' (IX * J)**2 
FP ( I )=FP{ J )A ( 1Y,< 3 
3 P ( 1 3 =GP ( I )-Y( 1 Y , K ) * *2 

2 j 3 ON I I *, Ur- 
23 MLAG-=MLAG+4 
3 'J . u I = 1 , M L A : j. 

MMI AP = USA IP-1 + ,. 

GP( I }=•!. 

3<! 15 J = _, 

<4= J+l~ . 

25 3 F ( 1 } =CP ( I 3 + ( .{ 1A, K-:)*X( 1Y, J 3 3 
WP( I )=CP( i 3/ ( FL JAT (MMINP) ) 

RMU.‘.=F Lull { NMi 1 F’ ) *GP ( 1 3 - r- P ( I )*TP{ I 3 
R3_Wi = SGRT( (FLj ,T{ Ai D P) *GP { 1 3 )-FP(I)**2> 



RL - ^. = 5GR T < { FI J. T ( MMINPJ-S'SP ( X ) ) -TP ( I ) **2 ) 

:o *f>( i )=k’]j;j/(Ri- .* *rD'*. n2) 

I F { i IY) - 0 J j fc, , : ; j 
2‘ 3 Pi IT 15 v * I iY 

.iS.i x , <' : CR. 1 j COR 3 ! LAUON X=*, 14,* Y= * , I 4//5X , *P*, 7X, *$UM X*, 

; *0U T Y*, . • *SUM X $ 0* , x x X , *SUM Y SQ*»liX,*C PR0D*,i3X,*R CCcF 

L*/) 

3* 1 lu 3 <t 

2 j 3 Hi !T 9-5,1 I Y 

99 : uA’U.| {//: X,*A.Jj OCORRLL -TIDM X= *» 13,* Y=*, I 3//5X, *P* , 7X, *SUM X*, 

l. T.:, *SU4 Y*, , ,.,*SUM X $.*, _»X,*SUM Y SQ*»1.*X,*C PROD*, 13X,*R COIF 

o*/) 

2 - 4 : jm f inu_ 

0 u »'u I = . , Y L a 9 , 

1 P=I 

5 J 3 R i • . T " 1 IP, IPli) ,Fp{ I) .SPU ),GP(I ), CPU ),RP(I) 

103 FORF XT { 4u',I5j3:’9.'j,F16*8) 

M L f“' ~7 u- ~~ M L *1 o “ _ 

33 DO .9 I = ^,;-;LAu 

FL P { I } = . 

DO - (j 0 -- c. , ^ L *j 

lb Ft_P(i }=F_p(I)+_.u*WP(J )^COS(^.j.4I59^7*FLOAT( { I- » ) *( J-i ) ) /FLOAT ( 

xMLAG) ) 

j9 FLP( i )=FLP { I )+rtP(l )+WP (Mi.A&i)*C3S(3.x416x*FLJAT < I-* ) ) 

JP( * )=0.46*FLP{x J + .S^FL-Mx) 

J P ( . v iL .-^G t ) = 0 • 4 *;*r LP { ML A G) + • 54*FL P ( ML AG I ) 

Dc --3 I = 2, ML AG 

43 J?m = ,.2 J *FLP(I-l ) + . 54*!-LP ( I ) + . 23*FL PC I + i ) 

3 R I .7 9/, IX, I Y 

93 = J R; ) ri7 ( / / _ X , *S PIC TR AL ANALYSIS OF* ,13,* VS*,lJ//* P*,7X,*CuV 

'.AT*,. *.X, *RAW S PcC TRA * ,5 .,*SMU0TH SPECTRA*, 5X,* FREQUENCY*//) 

TOT UP = ■ 

JO To I=;,'!L^, 

I P = I - 

= P. 4=FL DAY ( I? }/{ 1. *FL0AT( MLAS ) ) 

T JTUP = T jTliP + UPt i ) 

75 3 R 1 mT 1 > 0 ,IP,rt p in,FLP(l),UP(i),FR = Q 

CuR -)U= {. . *F LO .7 ( NS-.XIP ) /F-3AT ( MLAG ) )-2./3. 

SFF,~R = T JTUP/FL J- TIMLmGI) 

3 R 1 XT 1 >2, C0R'j-J,GF3AR 

102 SjRMAl ( / / i X , * p . R A M t T E RS FOR SPECTRAL AXALYSIS*//5X,*NU = *,Fx6.4,* 

* M . l, M .//) 

0 yl CO ITLIU- 

00 i 0 1 J 

403 CO. <7 I PJ_ 

STOP 
r 'VJ 



